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ABSTRACT 

This  manual  describes  MICROSCOPE ;  a  portable  FORTRAN  software  system  for 
the  analysis  of  multivariate  functions*  Given  an  interpolation  or 
approximation  scheme,  it  allows  the  following  questions,  among  others,  to  be 
answered:  Does  the  scheme  interpolate?  How  often  is  it  differentiable?  What 
functions  does  it  reproduce  exactly?  If  the  scheme  is  polynomial,  what  is  its 
polynomial  degree?  Where  is  the  smoothness  of  a  function  reduced?  Where  are 
the  bugs  in  a  FORTRAN  implementation? 
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SIGNIFICANCE  AND  EXPLANATION 


In  Computer  Aided  Geometric  Design  and  other  applications  of  Multivariate 
Interpolation  and  Approximation  it  is  often  necessary  to  build  ccxnplicated 
computational  procedures.  These  should  have  certain  properties  according  to 
the  mathematical  theory,  but  the  construction  may  have  failed  due  to 
oversights  or  simple  programming  errors.  This  report  describes  an  interactive 
portable  FORTRAN  software  system,  MICROSCOPE,  that  allows  the  following 
questions,  among  others,  to  be  answered:  Given  an  approximation  or 
interpolation  scheme,  does  the  scheme  interpolate?  How  often  is  the 
approxlmant  differentiable?  What  functions  does  it  reproduce  exactly?  The 
system  also  allows  the  identification  of  bugs  in  a  FORTRAN  implementation  of 
an  approximation  scheme,  often  pinpointing  the  faulty  portion  of  the  code 
within  narrow  limits.  MICROSCOPE  is  also  useful  in  the  corroboration  of 
theoretical  results,  in  guiding  a  theoretical  analysis,  and  for  classroom 
use.  In  addition  to  documenting  the  program  package,  this  manual  also 
contains  a  thorough  discussion  of  numerical  differentiation. 

For  ordering  information  contact  the  first  author. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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e.  About  Thia  Monuol 

This  manual  constitutaa  tha  comprahanaiva  daacription  of 
MICROSCOPE,  its  opsration,  applications,  foundotions,  and  innar  workings. 

Saction  1.1  dascribaa  what  kind  of  problams  con  ba  oddrassad  by 
MICROSCOPE.  Saction  1.2  dascribaa  tha  basic  opprooch  ond  is  fundamantol . 
Saction  1.3  dascribaa  tha  alphanunaricol  acraan  display  and  is  naadad  for 
undaratonding  tha  axomplaa  givan  throughout  thia  manual.  It  also 
introducaa  soma  farms  and  concapts  and  aarvas  as  a  preluda  to  the 
copabilitias  of  MICROSCOPE. 

Saction  2  gives  axomplaa  that  illustrate  four  major  applications 
(interpolation,  amaothnaas,  degree  of  precision,  ond  polynomiol  dagraa). 
Thaaa  do  of  course  not  exhaust  tha  potential  of  MICROSCOPE.  In  saction  3, 
ona  other  application  (discovery  or  i dent i f icot ion  of  o  point  where 
smoothness  is  reduced)  is  used  as  a  vehicle  for  introducing  tha  maehonica 
of  using  MICROSCOPE. 

Section  4  gives  a  detailed  description  of  all  commends  ond 
features,  organized  by  type  of  feature  or  task  to  ba  oceomp I i shad .  In 
ragulor  use,  this  will  probably  ba  tha  most  frequently  consulted  aaction. 

In  reading  sections  2  ond  3,  it  is  dasirobla  thot  o  working 
version  of  MICROSCOPE  is  ovoilobla  so  thot  bonds  on  axparianca  eon  ba 
goinad.  If  this  is  not  tha  cose  tha  instollotion  guide  (saction  5)  should 
ba  consulted  first.  Thot  saction  olso  contains  o  datoilad  discussion  of 
tha  portobiiity  feoturas  ond  possible  vorionts  of  MICROSCOPE. 

Appendices  ora  included  that  describe  soma  technical  features  in 
detail:  tha  differentiation  formulos,  tha  organization  of  tha  program 
pockaga,  and  tha  test  package.  Thera  should  olso  ba  on  oppandix  IV 
contoining  pan  plots  corresponding  to  tha  cruder  olphonumaricol  graphics  in 
tha  body  of  tha  manual.  Hawavar,  thot  appendix  is  not  machine  legible,  and 
is  distributed  separately. 

Tha  manual  itaalf  is  in  machine  readable  form  so  thot  it  con  ba 
distributed  together  with  the  packoga  and  con  ba  printed  locolly.  and  rood 
on  line  using  for  example  the  aaareh  features  of  o  text  editor.  It  wos 
ganaroted  by  o  MICROSCOPE  program,  with  tha  text  of  tha  manual  being 
"notes''  and  the  example  screen  diaploys  being  computed  by  MICROSCOPE  and 
automatically  embedded  into  the  text.  (Actually,  tha  MICROSCOPE  program 
generated  o  data  file  that  was  then  fed  into  the  text  formotting  program 
nnrisiFMT  Beebe,  1980.)  Machine  legibility,  however,  accounts  for  ths 
somewhat  awkward  notation  in  some  plocea.  For  example,  graek  letters  have 
been  given  abbreviations  in  terms  of  roman  letters. 

Definitions  are  denoted  by  writing  the  defined  term  in  single 
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quotation  morka.  Ordinary  quototiona  ora  in  doubla  quototion  morka. 

Tha  axonplaa  in  thia  nonuol  ora  raproducibla  ainea  tha  aat  of 
routinaa  uaad  for  thair  ganarotion  ia  auppliad  with  tha  pockoga  ond  con  ba 
uaad  for  taat  purpoaaa.  To  anaura  unifornity  of  raaulta  obtoinad  on 
diffarant  laochinaa,  o  foeility  hoa  boon  oddad  to  ainulota  rounding,  and  oil 
axomploa  hova  baan  eoiaputad  in  10  digita  orlthiaatic.  Mora  praciaaly.  If  x 
ia  tha  quantity  to  ba  roundad  to  0  digita,  aoy,  than  x  ia  raplocad  by 

—0  —0 
X  :>  (1-hapa  *10  )*x  +  opa  *10 

1  2 

whara  apa  and  apa  ora  random  nunbara  batwaan  -1  and  +1 .  Tha  oddliion  of 
1  2 

tha  apa  tarm  ia  not  atondard  but  oppropriota  In  tha  praaant  contaxt, 
2 

bacouaa  In  invaati gat  Iona  with  MlCBQSCOPE  amoll  numbara  ora  oftan  dua  to 

taking  diffaraneaa  batwaan  vary  cloaa  lorga  nunbara,  laoding  to  a 

coneal lotion  of  aignif leant  digita.  Tha  lorgar  tha  nognituda  of  x,  tha 
nora  inaignlf leant  tha  aaeond  tarn  will  baeona.  If  x  ia  roughly  1  or 
lorgar  (in  nognituda)  our  rounding  ia  aaaant lolly  aqut volant  to  tha 
atondard  rounding  (whara  apa  »  0). 

2 

Thia  nonuol  will  roraly  ba  rood  aaquantlolly  fron  front  to  book. 
Tharafora,  kay  idaoa  and  eoneapta  hova  aonatinaa  baan  raatotad  If  thay  ora 
eruelol  in  tha  givan  eontaxt.  It  ia  hopad  thot  thia  radundoney  will 

eontributa  to  tha  uaafulnaaa  of  thia  doeunant. 

Any  ertielana,  eonnanta  ond  auggaationa  about  MICROSCOPE  ora 

walcona  and  ahould  bo  diraetad  to  Patar  Alfald,  Daportnant  of  Mothanotiea, 
Univaralty  of  Utoh,  Solt  Loka  City,  Utoh  B4112,  801-581-6842.  Wa  ora 
particularly  Intaraatad  in  any  opplicotiona  of  MICROSCOPE  that  ora  not 
daacribad  in  thia  nonuol. 
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1.  Ifitroducl  ian 


1.1  Multivariate  Interoolat ion  and  Aporox imot i on 

This  manual  describes  a  FORTRAN  eoftwore  package.  MICROSCOPE, 
whose  main  applications  are  in  the  analysis  of  functions  of  several 
variobles  (or  ’surfaces')  os  they  occur  e.g.  in  the  opproximat ion  of 
multivariate  data  or  in  the  design  of  geometric  objects  such  os  the  body  of 
an  automobile  or  an  oircroft.  This  is  a  currently  very  active  research 
area  with  many  unsolved  ond  difficult  problems.  For  o  survey  see  Bornhill, 
198S,  ond  the  references  quoted  therein. 

MICROSCOPE  is  portable  to  the  extent  thot  it  has  possed  the  PFORT 
verifier  (Ryder,  1974)  with  the  exception  thot  it  uses  some  non-FORTRAN 
characters  (namely  ?,<,>, I ,:.:) .  Some  non-portable  optional  feotures  are 
supported  by  variants  of  the  code.  These  include  the  use  of  lower  case 
letters  in  Input  and  Output,  ond  the  use  of  modifications  of  on 
olphanumericol  display  on  a  CRT  terminal.  If  the  eoftwore  (<PL0T79>,  Beebe 
1979)  and  hardware  are  ovoiloble  more  sophisticated  ond  pleasing  grophicol 
displays  con  also  be  obtained.  However,  even  if  none  of  the  odditionol 
feotures  ore  provided  by  o  particulor  installation,  o  working,  if  somewhat 
crude,  version  of  MICROSCOPE  con  be  constructed.  See  the  installation 
guide  (section  S)  for  details. 

The  objects  thot  con  be  examined  by  MICROSCOPE  ore  functions  of 
from  one  to  three  independent  variobles.  Allowing  only  up  to  three 
variables  enabled  us  to  contain  all  numerical  information  on  the  CRT  screen 
and  proved  sufficent  for  all  coses  we  hove  encountered.  If  o  function  of 
more  than  three  voriobles  must  be  examined  MICROSCOPE  can  be  opplied  to  o 
eui table  restriction  of  the  function  to  a  subdomain  of  dimension  at  most 
three.  If  there  is  sufficent  interest,  future  vereions  of  MICROSCOPE 
capable  of  handling  more  than  three  variobles  will  be  provided. 

The  functions  of  interest  ore  usually  obtoined  by  opplying  on 
interpolation  or  approximation  scheme  to  discrete  or  transfinite  (i.e. 
infinitely  many)  data.  As  for  os  the  use  of  MICROSCOPE  is  concerned  it 
matters  little  if  a  function  is  obtained  by  interpolation  or  approximation 
(except  that  one  will  not  need  to  test  for  interpolation  in  the  latter 
case),  so  we  will  use  the  two  terms  interchangeably. 

We  will  refer  to  the  function  thot  is  being  examined  as  the 
‘trial  function*.  The  triol  function  may  be  different  from  the 
interpolating  or  opproximat ing  function  in  which  the  finol  user  of  a 
numerical  scheme  is  interested.  For  example,  it  is  often  usefui  to 
generate  data  from  a  'primitive  function'  ond  then  use  the  difference 
between  the  primitive  function  and  its  interpolont  os  the  trial  function. 


1.  Introduction 

1.1  Muitivariote  Interpolot ion  ond  Approximation 
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Thie  makee  diecontinultiee  in  derlvotivee  large  relative  to  the  derivative 
value  and  thereby  amplifiee  diecontinultiee  which  are  then  more  eaelly 
detected. 


Pour  typee  of  gueetlone  occur  moot  frequently  In  the  onolyele  of 
the  trial  function: 

-1-  Doee  the  echeme  Interpolate? 

-2-  How  often  le  the  trial  function  di f ferent loble?  (l.e. 
what  ie  ite  ‘degree  of  emoothneee’7) 

Note:Throughout  olmoet  all  of  thie  monuol  we  will  coneider  only  funetlone 
that  poeeeee  a  etep  diecontinuity  In  eome  derivotive,  for  thie  ie  the 
type  of  emoothneee  limitation  that  oecure  in  multivariate 
interpolation  and  approximation.  There  are  of  couree  other 

VS 

poeeibi 1 1 1 lee.  For  example,  the  function  f(x)  ax  ie  once 
differentiable  but  Ite  eecond  derivotive  exhibite  a  pole  rather  then  a 
etep.  The  teet  package  (eee  oppendix  III)  con  be  ueed  to  gain  eome 
experience  with  euch  funetlone.  In  porticulor,  the  examplee  in 
appendix  III  ehow  how  to  identify  ond  dietinguieh  euepe  and  polee. 

Which  funetlone  ore  reproduced  exoctly  by  the  echeme?  In 
particular,  which  le  the  moximum  degree  up  to  which  ell  polynomlale  ere 
reproduced  exactly?  In  ether  worde,  whet  le  the  'degree  of  (polynomial) 
preclelon'  of  the  echeme? 

-4-  If  the  trial  function  (or  eome  of  ite  derivotivee)  ore 
polynomial  what  ie  their  degree? 

It  ehould  be  obvioue  that  only  in  the  very  eimpleet  of  coeee 
onewere  to  theee  queetione  con  be  obtained  from  o  plain  diepley  of  the 
eurfoce  of  intereet.  Even  given  Infinite  dleploy  occurecy  (although  of 
couree  any  phyeicol  dleploy  ie  ultimotely  piecewiee  llneor  or  even 
diecrete),  it  eeeme  impoeeible  to  dietinguieh  optically  a  once 
differentiable  function  from  one  that  ie  twice  differentiable. 


_ The  Bail 6  AaBrsaeli 

In  the  deeign  of  on  approximation  echeme  one  ueuolly  knowe  the 
‘critical  eete’  where  the  phenomeno  of  intereet  euch  oe  diecontinultiee  in 
derivotivee,  or  interpolation  to  certain  data,  may  occur.  For  example,  in 
0  bivariate  interpolation  echeme  defined  on  a  t riongulot ion,  critical  eete 
ore  ueuolly  edgee  of  trionglee,  poeeibly  internol  edgee,  and  verticee  of 
trionglee.  In  o  tri variate  echeme  defined  on  a  teeeelotion  into 
tetrohedro,  critical  eete  are  focoe,  edgee,  maybe  internal  facee  and  edgee, 
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and  internal  and  external  vertices.  If  the  criticol  set  is  not  known  it 
can  be  discovered  with  MICROSCOPE  (see  section  3),  but  for  the  present  we 
will  assume  it  is  known. 

The  bosic  idea  of  MICROSCOPE  is  very  simple:  Pick  o  random  point 
(the  'point  of  examination')  in  a  specific  critical  set.  We  will  refer  to 
that  set  as  'the  front'.  Then  determine  o  'direction  of  investigation'  and 
consider  a  line  through  the  point  of  examination  in  the  direction  of 
investigation.  This  is  the  'line  of  investigation'.  Depending  upon  the 
application,  the  line  of  investigation  may  be  contoined  in  the  front,  or  it 
may  be  at  an  angle  (not  necessorily  o  right  one)  to  it.  We  will  call  any 
derivative  in  the  direction  of  investigation  o  'tongentiol  derivotive'. 
Any  pure  derivotive  in  a  direction  other  than  the  direction  of 
investigation  is  a  ‘cross  derivative'  ond  its  direction  is  the  ‘cross 
direction'.  MICROSCOPE  numerically  approximates  the  relevont  derivotives 
(pure  or  mixed)  and  displays  the  approx imot ions  on  the  CRT  terminol. 
Depending  on  the  resulting  picture,  one  may  then  toke  further  action,  e.g. 
pick  another  point  of  examination  or  onother  direction  of  investigation, 
consider  other  derivatives,  or  improve  the  existing  plot  by  oltering  the 
underlying  numerical  parameters. 

In  using  MICROSCOPE  it  is  important  to  choose  the  underlying 
domain  and  the  primitive  function  so  as  to  avoid  artifacts  due  to  the 
regulority  of  the  domain  or  pecul iori t ies  of  the  primitive  function,  ond  to 
pick  the  points  of  examination  randomly  or  orbitrorily  within  the  critical 
sets.  Usuolly  different  types  of  critical  sets  will  be  present  and  of 
course  one  will  need  to  consider  at  least  one  representative  of  soch  set. 

Numerical  differentiation  is  a  notoriously  ill-conditioned 
numerical  process.  However,  in  MICROSCOPE,  the  appeoronce  of  round-off 
errors  is  immediately  apparent  because  the  points  on  the  graph  of  the 
relevant  derivative  ore  scattered  across  the  screen  in  a  random  fashion. 
Moreover,  by  plotting  derivatives  of  o  high  degree  with  o  large 
discretization  parameter  the  limitations  imposed  by  round-off  errors  can 
usuolly  be  overcome.  For  example,  if  the  expected  discontinuity  in  o 
specific  derivative  is  too  unpronounced  to  be  visible  at  the  lorge  value  af 
the  discretization  porameter  necessitated  by  round-off  errors,  then  higher 
order  derivatives  con  be  computed  (on  on  even  coarser  discretization). 
These  will  be  approximations  of  the  Diroc  Del  to  function  and  its 
derivatives,  and  usually  will  pinpoint  even  weak  discontinuities. 

Numerically,  we  proceed  os  follows.  Let  F  denote  the  trial 
function,  P  the  point  of  examination,  and  d  the  direction  of  investigation. 
MICROSCOPE  displays  derivatives  of  F  in  the  direction  of  d  along  the  line 
of  investigation  with  P  being  the  center  point.  A  built-in  facility  allows 
the  replacement  of  F  with  a  derivative  of  F  in  some  other,  independent, 
direction.  More  compiicated  mixed  partial  derivatives  cart  be  built  by 
modifying  the  trial  function,  and  MICROSCOPE  offers  some  help  with  that  os 
well.  However,  for  the  present  we  will  restrict  our  attention  to  a 
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tongantlol  darlvatlv*  of  F. 

Th*  diraetlonol  darlvotiva*  ora  approx imotad  at  pointt 

P(l)  -  P  -f  i•••d,  whora  i  -  -N,  -Mf-1 . N-1,N 

ond  tho  quantity  ■  >  6  datarmlnaa  tha  spoclng  of  tho  pointa.  It  will  ba 
rofarad  to  a«  tha  'display  Intsrvol'  or  'intsrvol*  for  short.  Ths  nunbsr  N 
dstsrmlnss  ths  rsaolution  of  ths  display  (as  wall  os  tha  numsrieal  sffort). 
At  soch  point  P(l).  ws  approximato  ths  k-th  dsrivotivs  of  F  In  ths 

(k) 

dirsetlon  d  ot  ths  l-th  point  by  a  quontity  0 

I 

Ths  rangs  of  k  In  ths  prsssnt  vsrslon  of  MICROSCOPE  Is  from  zsro 
to  six.  Ws  hovs  found  this  sufficisnt  for  oil  of  our  applications.  Ths 
(k) 

quontitiss  0  ors  dstsrminsd  by  numsrieal  di f fsrsnt lot  ion  (i.s.  by 

I 

dl f fsrsnt lot ing  on  Intsrpolot ing  polynomial).  Ths  formulos,  givsn  bslow, 
wsrs  dstsrminsd  by  ths  following  critsrlo: 

-1-  Control  Olffsrsness  should  bs  smploysd  in  ordsr  to  ovoid 
art  I  foot  dus  to  osymmstry. 

-2-  Ths  diserstixotton  poramstsr  (which  doss  not  nsesssorily 
squol  s)  should  bs  such  that  ths  numbsr  of  points  offsetsd  by  ths  suspsctsd 
discontinuity  (in  soms  dsrivotivs)  is  indspsndsnt  of  ths  ordsr  of  ths 
dsrivotivs. 

-3-  Othsr  than  ths  abovs,  ths  dsgrss  of  ths  intsrpolotlng 
polynomial  should  bs  as  low  os  possibis. 

Furthsrmors,  wo  pieksd  discrst Izot Ion  formulas  thot  uss  points  at 
squally  spaced  Intsrvols.  This  may  bs  chongsd  in  future  vsrsions  of 
MICROSCOPE- 

DsnotIng  ths  'diserstizotion  poramstsr'  by  h  >6,  ths  above 
criteria  give  rise  to  ths  following  formulas: 

(0) 

0  -  F5 

I 

(1) 

Q  -  (-F1  +  Fg)/(2»h) 


1.  Introduction  page  11 

1.2  Ths  Basic  Approach 


V 


MICROSCOPE:  A  Softwor*  SyttM  for  Multivoriot*  Anolyaio  pog*  12 

Uaor'a  Manual  Varaion  of  April  1984 


(2)  2 

0  -  (n  -  2F5  +  F9)/h 


(3)  3 

0  .  4*(-F1  +  2F3  -  2F7  +  F9)/h 


(4)  4 

0  -  16»(F1  -  4F3  +  6F5  -  4F7  ♦  F9)/h 


(5)  5 

0  -  243«{-F1  4  4F2  -  5F4  4  5F6  -  4F8  4  F9)/(2h  ) 


(6)  6 

0  -  729«(F1  -  6F2  4  15F4  -  -2eF5  4  15F6  -  6FB  4  F9)/h 

i 

whara 

FI  -  F(P(I)  -  hd) 

2hd 

r2  -  F(P(1) - ) 

3 

hd 

F3  -  F(P(i)  -  — ) 

2 

hd 

F4  .  F(P(i) - ) 

3 

F5  -  F(P(i)) 

hd 

F6  -  F(P(i)  4  — ) 

3 

hd 

F7  .  F<P(i)  4  — ) 

2 

2hd 

F8  -  F(P(i)  4  - ) 

3 
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F9  -  r(P(i)  hd) 

The  firit  four  of  tht  obov*  formulo*  can  bo  found  in  Abromowitz 
and  Stogun,  1960,  p.  914,  and  tho  loot  two  con  bo  dorivod  by  tho 
tochniguoo  doocribod  thoro.  Highor  ordor  formutao  con  of  couroo  aloo  bo 
dor i vod. 


(k) 

It  io  intorooting  to  noto  thot  0  conoidorod  ao  a  function  of  x 

i 

(or  a  function  of  a  continuouoly  varying  paromotor  i)  It  d  function  of  tho 
oonio  dogroo  of  omoothnooo  ao  tho  trial  function  itoolf.  Just  tho  numbor  of 
points  of  roducod  oflioothnooo  it  incrooood,  thoro  boing  ono  corrooponding  to 
oach  point  in  tho  difforont lotion  fornulo  potting  through  tho  point  of 
roducod  onoothnooo.  Thus  numoricoi  difforontiot Ion  con  bo  thought  of  oo  o 
onoothi ng  procooo.  (Thio  controoto  otorkly  with  ito  numoricoi  property  of 
o  roughing  proeoduro  duo  to  round-off  offocto.)  Tho  figuroo  in  aoctlon  4.2 
illuotroto  tho  omoothlng  proportioo  of  numoricot  di f foront lot  ion. 

For  a  moro  dotoi ltd  diocuooion  of  tho  truncotion  ond  round-off 
proportioo  of  thooo  formuloo  too  Appendix  1. 


1.3  Tho  Alohonumtrieol  Dioolov 

In  thio  oubooction,  tho  alphanumerical  tcroon  dioploy  gonorotod 
by  MICRQSCQPr  It  doocribod  in  dotoi I .  Tho  purpose  of  tho  dooeription  io 
twofold:  First,  it  io  noeoooary  for  on  understanding  of  tho  oxomploo 
tprinklod  throughout  this  manual.  Second,  however,  it  also  servos  os  on 
illustration  of  some  of  tho  capabilities  of  MICROSCOPE  and  os  a  first 
Introduction  to  tho  power  and  versatility  of  tho  pockogo. 

During  a  MICROSCOPE  oosoion  tho  CRT  screen  will  usuolly  bo 
occupied  by  tho  'alphanumerical  display'  like  ths  ono  illustrated  in 
figure  1.  (Moro  oophisticatod  graphical  displays  can  also  bo  obtoinod,  see 
sections  4.14  and  5.) 
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CD:  deg  -  1  dir  -  (  8.0808680408.  1.0008000480)  ch  -  6.60800-84 

Point  •  (  0.0008800400.  1.8080000480)  e  -  5.00800-83 

Direction  >  (  1.0808000480.  0.0808000480)  h  -  3.08000-02 

FO  (  9.390-09,  1.27D-02)  F2  (  1.200-01.  2.220400)  F3  (-1.200481,  1.210481) 
I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  current  CALLS  -  174 

Figure  1  Oemonst rat  ion  of  the  Alphonumericol  Oisploy 


2 

In  this  example,  the  trial  function  is  f(x,y)  *  x  .y.obs(x*y). 
We  analyze  a  crass  derivative  of  it  in  the  direction  (0,1).  The  direction 
of  investigation  is  (1,0).  and  the  point  of  examination  is  (0,1).  Because 
af  the  absolute  value  term,  the  cross  derivative  is  twice  but  not  three 
times  tangentially  differentiable  at  the  point  of  examination. 

We  explain  the  display  starting  at  the  top  and  referring  to  lines 
either  by  their  position  or  by  the  first  word  contained  in  them. 

The  display  is  naturally  divided  into  two  ports:  The  first  15 
lines  constitute  the  'graphical  display’,  and  lines  17  through  21  form  the 
'numerical  display'.  Line  16  separates  the  two  disploys  and  is  present 
only  if  there  is  sufficient  space  available. 

-1-  The  Grophicol  Display 

Every  MICROSCOPE  display  shows  the  graphs  of  o  set  of  tongentiol 
derivatives  of  degrees  between  0  ond  6.  The  function  that  is  being 
differentiated  tangentially  is  called  the  'display  function'.  It  may  be 
the  trial  function  itself  (the  most  frequently  occuring  cose),  or,  os  in 
this  example,  o  pure  cross  derivative  of  o  degree  between  1  and  6. 
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Distance  along  the  horizontol  oxis  corresponds  to  displacement  in 
the  direction  of  investigation  from  the  point  of  exominotion  (which 
corresponds  to  the  center  of  the  horizontal  axis).  The  ecale  on  thot  axis 
has  the  eome  meaning  for  all  graphs  shown.  The  vertical  axis  corresponds 
to  function  values,  but  the  verticol  scale  is  different  for  each 
derivative,  and  is  chosen  such  thot  the  graph  fills  the  entire  verticol 
extent  of  the  graphical  display.  This  normalization  is  corried  out  and 
modified  outomot iccl ly  os  the  invest igot ion  proceeds.  A  period  denotes  ths 
graph  of  the  display  function,  ond  the  digit  k,  soy,  where  k  is  between  1 
and  6  ond  denotes  the  k-th  tangential  derivotive.  In  this  exomple,  the 
dieploy  function  and  ite  eecond  and  third  tangent iai  derivatives  ore  being 
shown. 


Sometimes  some  derivotive  is  more  interesting  than  othsrs,  ond 
then  it  may  be  useful  to  accentuate  it  with  asterisks.  In  the  above 
example,  the  third  derivative  hos  been  eo  emphasized.  Usuoiiy,  the  graph 
of  0  given  derivative  overwrites  the  groph  of  ony  lower  order  derivative. 
However,  the  graph  of  an  occentuated  derivative  overwrites  all  other 
graphs. 


Both  the  number  of  lines  ond  the  number  of  columns  in  ths  display 
con  be  set  according  to  specific  purposes.  Throughout  this  manual  we  will 
use  7S  columns  and  IS  rows  for  the  grophicol  dieploy  os  this  is  appropriate 
for  the  width  of  the  printed  output.  A  width  of  79  might  be  preferable  for 
a  standard  CRT  terminal  disploying  8B  columns:  135  columns  ond  a  lorger 
number  of  rows  might  be  used  for  printed  output;  ond  a  smoller  number  is 
sometimes  useful  for  preliminary  investigations  when  the  triol  function  is 
expensive  to  evaluate.  An  even  number  of  columns  is  normally  undesiroble 
because  then  the  point  of  examination  does  not  lie  in  the  precise  center  of 
the  display. 

Line  6  in  the  example  contoine  a  horizontal  scole  counting  from 
the  center  to  the  left  and  right.  This  will  usually  be  omitted  so  as  not 
to  clutter  the  display,  but  it  is  sometimes  useful  for  the  identification 
of  points  of  interest.  The  distance  from  the  center  can  be  used  to  inquire 
about  numerical  values  of  certoin  derivatives,  or  to  shift  the  graph  left 
or  right  os  required  by  the  c i rcumstonces.  The  scole  con  be  turned  on  or 
off,  or  replaced  with  a  line  of  hyphens.  The  groph  of  any  function 
overwrites  the  horizontal  scale.  The  center  of  the  dieploy  con  optionally 
be  morked  with  o  "•f"  sign.  This  is  sometimes  useful  to  pinpoint  symmetry 
properties.  In  the  above  example  the  groph  of  the  third  derivative  posses 
through  the  center  of  the  display.  The  center  mork  overwrites  the  graph  of 
ony  function  passing  through  it. 

The  central  column  of  colons  morks  the  column  corresponding  to 
ths  point  of  examination.  The  two  columns  of  I's  outline  'the  window', 
i.e.  that  set  of  points  at  which  the  numerical  differentiation  is  offected 
by  the  behavior  of  the  trial  function  ot  the  point  of  exominotion.  (This 
is  not  always  quite  true  when  cross  derivatives  are  being  examined,  see  the 
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discussion  in  section  4.11).  for  every  point  in  the  window,  the  interval 
about  the  point  spanned  by  the  differentiation  stencil  contains  the  point 
of  examination.  The  differentiation  formulas  are  chosen  such  that  the 
window  is  identical  for  all  tangential  derivotives.  Obviously,  for  the 
display  function  itself  the  window  contains  only  the  point  of  examination. 
For  its  derivatives,  the  window  contains  all  points  thot  con  be  written  os 
P  -f  td,  where  abs(t)  does  not  exceed  the  value  of  h/s,  (and,  os  before,  P 
is  the  point  of  examinotion,  d  is  the  direction  of  investigation,  h  is  the 
discretization  parameter,  and  s  is  the  interval).  We  define  w:»2*h/s  to  be 
the  ‘window  width'.  Obviously,  by  picking  the  window  width  <  1,  it  is 
possible  to  generate  displays  in  which  the  only  point  contained  in  the 
window  is  P  itself.  However,  for  several  reosons  it  is  usually  preferable 
to  use  a  larger  window  width.  This  ospect  is  discussed  in  more  detail  in 
section  4.1.  The  default  value  of  the  window  width  is  12,  os  in  the 
example,  but  this  con  be  changed  os  needed. 

Without  examining  the  numerical  display,  it  is  apporent  from  the 
graphical  display  that  some  function  is  being  investigated  that  is  twice 
but  not  three  times  differentiable  since  the  (occentuoted)  third  derivotive 
shows  a  clear  step  discontinuity.  Note  how  thot  step  is  smeored  out  over 
the  window. 


-2-  The  numerical  display 

The  first  line  of  the  numerical  disploy  starting  with  "CO:" 
indicates  that  the  derivatives  are  mixed  portiols,  i.e.  they  ore 
tangential  derivatives  of  a  cross  derivative  (of  the  trial  function).  This 
line  is  absent  if  the  trial  function  itself  is  being  examined,  and  the  line 
carries  only  information  that  is  pertinent  to  the  cross  derivative.  In  the 
present  example,  the  degree  of  the  cross  derivative  ("deg")  is  1,  the  cross 
direction  is  (0.1),  and  the  discretizotion  parameter  used  for  computing  the 
derivative  (denoted  by  "ch"  corresponding  to  "h"  for  the  tongential 
derivative)  is  60-4.  The  formulas  for  computing  cross  derivotives  ore 
equivalent  to  those  for  computing  tongentiol  derivatives. 

The  line  starting  with  "Point"  describes  the  point  of  examination 
(which  is  (0,1))  and  the  intervol  (which  is  SD-3) .  This  display  implies 
that  the  trial  function  is  bivoriote.  Appropriote  modifications  toke  place 
for  functions  of  1  or  3  variables.  The  next  tine  starting  with  "direction" 
describes  the  direction  of  investigation  (which  is  (1,0))  and  the 
discretizotion  parameter  (which  is  30-2).  Notice  that  the  window  width  is 
2h/s  •  12  which  is  consistent  with  the  grophicol  disploy 

The  "Point"  and  "Direction"  descriptions  ore  always  present. 

The  next  line,  starting  with  ''F0'',  gives  the  ranges  of  the 
displayed  derivatives.  Ranges  are  labeled  Fk  where  k  is  the  order  of  the 
tangential  derivative.  For  example,  F2  corresponds  to  the  range  of  the 
second  (tangential)  derivative  (of  a  first  order  cross  derivative)  which 
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h«r«  i*  (0.12.2.1).  If  mor*  function*  (up  to  ••van)  or*  di*play*d,  than 
thair  ronga*  or*  alao  givan,  adding  up  to  too  lino*  to  th*  nunarlcol 
dlaploy  a*  nacaaaory. 

Th*  naxt  and  loat  I  in*,  atorting  with  "1/0",  contain*  information 
about  th*  praaant  configuration  of  MICROSCOPE. 

Th*  10  Intagar*  following  "I/O"  or*  th*  FORTRAN  chonnal  numbara 
0*  thay  or*  amployad  by  MICROSCOPE  at  th*  dlaploy  tima.  Thay  corraapond  to 
th*  following: 


-1-  ‘Input  Oavlea'  (har*  25),  I.*.  th*  davic*  from  which 
command*  or*  raod.  In  on  Intaroctiva  aaaaion,  thi*  will  uauolly  b*  th* 
tarminol.  But  It  moy  olao  b*  o  fll*  of  MICROSCOPE  commond*  thot  ganarot*  o 
apacific  aatup,  or  o  documantot Ion  Ilk*  thi*  manual. 

-2-  'Output  Oavlea*  (har*  6),  i.a.  th*  davic*  raeaiving  th* 
command  prompt*  ond  any  halp  information.  Uauolly,  thi*  will  olao  b*  th* 
tarminol,  or,  if  UICROSCOPC  la  baing  dr  I  van  by  command*  in  o  fila,  it  moy 
b*  0  Null  or  Dummy  davic*. 

-5-  ‘Grophle*  Oavlea'  (har*  olao  6),  I.a.  th*  davic*  raeaiving 
th*  alphonumaricol  dlaploy  lltuatrotad  in  figur*  1.  Thi*  will  ogoin 
uauolly  b*  th*  tarminol,  but  it  moy  olao  b*  o  aoporot*  dlaploy  davic*  or  o 
fll*  raeaiving  informotlon  for  aubaaquant  printing. 

-4-  ‘Halp  Oavlea'  (har*  1),  i.a.  th*  davic*  containing  th* 
MICROSCOPE  halp  flla,  which  I*  auppllod  with  th*  pockoga.  MICROSCOPE  ho* 
aavarol  intaroctiva  halp  focilltia*  ond  th*  information  naadad  for  running 
thorn  ia  raod  out  of  thi*  fila  whan  MICROSCOPE  ia  firat  collad.  Uauolly  you 
will  not  b*  concarnad  with  praciaaly  whar*  MICROSCOPE  racaivaa  it*  halp 
Information.  Howavar,  th*  chonnal  numbar  t*  printed  in  th*  dlaploy  to  halp 
pravant  on  occldantol  ovarwrit*  of  th*  halp  fll*. 

-5-  'Recording  Davic*'  (har*  2),  o  davic*  thot  con  racaiv* 
aalactad  information  from  th*  MICROSCOPE  aaaaion,  a.g.  copia*  of  th* 
acraan  dlaploy  or  'not**'  that  explain  th*  dlaploy*.  Thi*  monuol  i*  th* 
raault  of  reading  o  oat  of  inatructiona,  i.a.  o  MICROSCOPE  'program'  from 
th*  Input  Oavlea  (a  data  file),  procaaaing  It  and  writing  th*  raault*  onto 
th*  Recording  Davic*. 

-6-  'Raatort  Oavlea’  (har*  3).  At  any  time  during  o  aaaaion. 
th*  currant  configuration  of  MICROSCOPE  con  b*  aovad,  allowing  o  lotar 
raatort  from  th*  currant  atotua  of  th*  invaatigotion.  Thi*  i*  occompliahad 
by  writing  th*  volua*  of  all  ralavont  voriabla*  into  o  file  pointed  at  by 
th*  Raatort  davic*. 
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-7-  'Logging  Davict'  (hara  27).  It  ia  poasibla  to  kaap  o  diary 
of  th  MICROSCOPE  aasaion  listing  tha  conMonds  ond  date  antarad  os  wall  os 
th  rasults  affactad  by  than.  This  feature  is  activated  by  setting  tha 
Logging  Device  number  to  a  non-zero  value. 

-8-  'Point  Loading  Device'  (hare  28).  A  device  allowing  you  to 
rood  points  from  o  file,  rather  than  having  to  type  them  in. 

-9-  'Direction  Loading  Device'  (here  29).  Similar  to  -8-  for 
loading  the  direction  of  investigation. 

-10-  'Cross  Direction  Loading  Device'  (here  30).  Similar  to  -8- 
for  loading  the  cross  direction. 

The  firat  six  I/O  channel  numbers  are  parameters  af  the  master 
MICROSCOPE  routine  (see  section  3.1)  and  con  each  be  changed  during  a 
session.  The  last  four  channels  are  zero  (i.e.  unused)  by  default,  and 
con  be  set  to  non-zero  (i.e.  active)  values  during  a  MICROSCOPE  session. 
It  is  of  course  unnecessary  to  hove  all  channel  numbers  distinct.  For 
example,  all  three  loading  devices  might  be  identicol,  allowing  for  the 
possibility  of  having  just  one  list  of  points,  directions,  and  cross 
di rect ions. 

Following  the  device  numbers  is  o  flog  ”NRML  on"  or  "NRML  off" 
(hare  on),  indicating  whether  the  direction  of  investigation  entered  into 
MICROSCOPE  has  been  normalized  or  not.  The  "Direction"  line  always  gives 
the  direction  as  entered  by  the  user.  If  NRML  is  off,  the  direction  of 
investigation,  d,  is  identical  to  that  printed  in  the  display,  otherwise  it 
is  that  printed  in  the  display  normalized  to  hove  unit  Euclidean  length. 
With  tha  normalization  on,  tangential  derivotives  ore  the  standard 
directional  derivatives  as  they  ore  defined  in  Calculus  Textbooks.  With 
the  normalization  off,  they  ore  the  more  versatile  Cotaoux  derivotives  as 
they  are  used  frequently  in  multivoriote  interpolot ion  ond  approximation. 
The  NRML  flog  opplies  similarly  to  the  cross  direction. 

Following  the  NRML  flog  is  the  word  "current”.  This  indicotes 
that  the  numerical  display  corresponds  to  the  grophics  display.  To  explain 
why  thia  might  not  be  so  requires  o  brief  excursion  into  the  Organization 
of  MICROSCOPE.  When  drastic  parameter  chonges  are  taking  ploce  it  is  often 
more  efficient  to  accumulate  them  before  generoting  the  new  disploy. 
Usually  in  that  case  the  al phanumericol  display  is  scrolled  while  commands 
are  being  given,  and  vanishes  from  the  screen.  However,  sometimes  it  is 
preferable  to  view  the  old  display  rather  then  the  commonds  given  so  far. 
In  that  case,  the  display  can  be  recalled  on  the  screen,  and  the  word 
"currant"  is  replaced  by  "GO  pndg"  (i.e.  the  GO  command  needs  to  be  given 
to  incorporate  into  the  graphical  display  changes  of  porameters  that  may 
already  be  shown  in  tha  numerical  display). 
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Tha  loat  itan  in  tha  loat  lino  ia  "CALLS  >  174“  whora  174  la  tha 
nuiabar  of  avoluotiona  of  tha  trial  function  axpandod  up  to  thia  point  in 
tha  yiCRQSCOPE  aaaaion.  Thia  olloaa  Monitoring  tha  computational  effort 
apant.  A  typical  Multivariate  aurfoca  will  be  axpanaiva  to  avoluata,  and 
ita  coat  will  uaually  doMinota  tha  coat  of  running  UlCROSCOPt:  and  tharaby 
datarMina  haw  long  you  hava  to  woit  for  a  roponaa  froM  tha  tarMinal. 
Tharafora.  a  aignlf leant  effort  hoa  boon  Moda  to  kaap  tha  nuMbor  of 
nacaaaary  function  avoluotiona  aMOII.  Thia  feature  eontributaa  in  no  email 
Maoaura  to  the  complexity  of  tha  tilCBOSCOPE  aourca  coda.  Tha  example 
dlaploy  hoa  75  eolumna  and  wa  ora  axainining  a  firat  order  croaa  darivotiva 
of  tha  trial  function.  Each  avoluotion  of  tha  croaa  darivotiva  raquiraa 
two  avoluotiona  of  tha  trial  function.  Diaploying  tha  groph  of  tha  croaa 
darivativa  alone  would  tharafora  require  o  value  of  CALLS  equal  to 
2*75  •  158.  Note  that  computiig  ond  diaploying  the  aacond  and  third 
tongantiol  darivotlvaa  of  the  croaa  darivotiva  raquiraa  on  additional 
effort  of  only  18X.  Thia  ia  not  o  trivial  occompi iahmant  oa  moy  be 
auggoatad  by  tha  fact  that  tha  aoma  dlaploy  with  tha  elightly  different 
window  width  of  16  would  require  530  avoluotiona.  For  o  detailed 
diaeuaaion  of  tha  nufflarical  effort  oaaoelotad  with  vorioua  window  widtha 
aaa  aaction  4.3. 

Notice  that  in  tha  above  example  wa  idantifiad  a  atap 
diecontinuity  in  a  mixed  partial  darivotiva  of  fourth  order  (o  third  order 
tongantiol  darivativa  of  o  firat  order  croaa  darivotiva)  without  fuaa  or 
doubt.  If  thia  woa  o  real  I nvaat Igot ion  wa  could  atata  now  that  tha 
underlying  achama  In  any  coaa  ia  at  moat  three  timaa  different  iobta. 


With  MICROSCOPE,  one  anawara  quaationa  Ilka  thoaa  in  aaction  1.1 
for  apacific  axomplea.  Thua  only  nagotiva  raaulta  (a.g.  thot  o  achama  ia 
BOl  di f farant i able)  can  be  proved  atrictly.  However,  poaitiva  raaulta  can 
alao  be  obtained.  If  tha  axomplea  ora  choaan  corafully  to  exclude 
ortifocta  and  make  tham  rapraaantotiva,  our  progrom  allowa  for  anawara  to 
tha  above  quaationa  with  a  dagraa  of  confidence  thot  bordara  on  cartaintyl 
Naturally,  a  computational  opprooch  doaa  not  raploca  a  rigoroua  proof. 
However,  in  tha  outhora  axperianca  it  ia  often  poaaibla  to  confirm  (or 
ahattar)  fuzzy  notiona  with  MICROSCOPE,  and  to  generate  further  inaight 
into  tha  phenomena  under  conaidarat ion.  Eneouroging  MICROSCOPE  raaulta  may 
help  in  building  up  tha  atomina  needed  for  attempting  o  rigoroua  proof,  ond 
datollad  inveat igot iona  moy  help  in  providing  idaoa  for  carrying  out  tha 
onalyaia. 


Often,  of  couraa,  MICROSCOPE  wilt  be  employed  in  hindaight,  whan 
tha  onalyat  already  knowa  tha  onawara  to  quaationa  like  tha  obove  from 
theoretical  raaaoning.  Tha  axaminotion  with  MICROSCOPE  than  aarvea  to 


corroborate  tha  raaaoning,  and  olao,  more  mundanely  but  no  laaa 
importantly,  to  eliminate  buga  in  tha  implamentot ion. 
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In  tho  outhora  axparianca.  MICROSCOPE  provad  Invaluabla  oa  a 
dabugging  tool,  oftan  pinpointing  within  narrow  limita  thot  port  of  a  code 
thot  woa  at  fault,  but  alao  aonotimaa  uncovering  flowa  in  tha  theoretical 
approoch,  and,  on  a  few  exhilarating  occaaiona.  lending  to  the  diacovery  of 
new  and  unlookad  for  raaulte  that  could  then  be  proved  theoretically.  Thua 
MICROSCOPE  ia  uaeful  in  the  identification  oa  wall  oa  in  the  verification 
of  reaulta. 

In  the  davalopnent  of  MICROSCOPE  we  were  occaaionolly  confronted 
with  tha  objection  that  our  alphonunerical  diaplaya  ore  rother  cruda.  Our 
primary  roaponao  to  thia  ia  thot  the  very  purpoae  of  MICROSCOPE  conaiata  of 
tronaloting  vary  aubtle  propertiea  into  diaploya  that  moke  them  glaringly 
apparent.  A  large  part  of  our  intoreat  centera  around  atop  functiona,  and 
it  aeama  fair  to  diaplay  them  aa  auch.  Indeed,  any  curve  fitting  might 
diagulae  the  vary  faoturea  in  which  we  are  intereated.  However,  we  do 
provide  on  interface  to  a  aophiat icoted  grophica  pockoge.  See  aectione 
4.14  and  S  for  more  information. 
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2.  AoDlicotlons 

In  this  section,  we  illustrate  the  main  applicotlons  of 
MICROSCOPE.  Consider  the  simple  example  of  Interpolating  to  a  primitive 
function  p  by  o  cubic  spline  S,  say,  occording  to  the  conditions: 


2  3 

S(x)  -  00  +  o1»x  +  a2»x  +  a3«x 

and 

2  3 

S(x}  -  be  4  b1*x  4  b2«x  4  b3*x 
where  the  coefficients  are  defined  by 


If  K  <  1 

otherwi se 


S(e)  -  p(e).  S(1)  -  p(i),  S(2)  -  p(2),  S'(2)  -  p'(2),  S"(e)  -  e 


and  S  is  twice  differentiable  at  x  >  1.  Thus  we  require  interpolation  to 
'position'  (i.e.  function  values)  ot  0,1,2,  and  interpolot ion  to  the  first 
derivative  at  2  by  a  piecewise  cubic  function  thot  is  twice  continuously 
differentiable.  At  0  we  impose  odditionolly  a  "free  end"  condition. 
(There  is  an  extensive  literature  on  Splines,  see  e.g.  ds  Boor,  1976,  or 
Schumaker,  1981) 

For  our  present  purposes  we  assume  that  p(x)  is  the  exponential, 

i  .e. 


X 

p(x)  -  e 

where  e  is  the  base  of  the  natural  logarithm.  This  function  oppeors  to  be 
sufficiently  non-polynomial  to  ovoid  the  introduction  of  artifacts. 

A  simple  calculation  shows  thot 

00  -  1 

2 

of  -  (-2ee  4l2*e-9)/7 
a2  -  0 

2 

o3  -  (2*e  -5*e42)/7 
and 

2 

b0  -  (See  -16*e4l2)/7 
2 

bi  -  (-17»e  460«e-24)/7 
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2 

b2  -  (15»e  -48«e+15)/7 
2 

b3  ■  (-3»e  +11«e-3)/7 

In  the  following  exoieplee  we  will  verify  thot  S  with  the  above 
coefficients  doee  indeed  poeeeee  the  required  propertiee. 


2.1  Interpolation 

The  boeic  idea  of  verifying  interpolot ion  ie  to  plot  the 
difference  between  the  priieitive  function  and  ite  interpolont.  and  to  check 
that  the  error  and  ite  appropriate  derivotivee  ore  indeed  zero.  So  we 
incorporote  the  exponential  and  S  oe  defined  obove  into  MICROSCOPE  (see 
section  3.1  for  detoile).  ond  start  on  interoctive  eeeeion.  After  giving 
the  appropriate  coemonde  we  obtoin  the  following  picture  illuetroting  the 
si tuot ion  at  x  ■  0: 


+  I 


I 

I 

I 


Point  •  0. 9909000+90  e  -  5.00000-93 
Direction  •  1.0090000+00  h  -  3.00000-02 
FO  (-3.330-02,  6.760-02) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRMl  on  current  CALLS  -  249 


Figure  2  Interpolot ion  to  Position  at  Zero 


We  notice  that  the  ronge  of  the  error  is  from  -0.0333  to  0.0676. 
Since  the  graph  of  the  error  function  passes  through  the  vertical  axis 
slightly  below  the  center,  this  ronge  ie  consistent  with  the  error  being 
zero,  i.e.  with  interpolot ion.  A  more  relioble  confirmation  con  be 
obtained  by  asking  for  the  volue  of  the  displayed  function  at  the  center 
(i.e.  0)  of  the  display.  This  procedure,  which  is  not  shown  here,  yields 
o  value  of  -9.40-11,  which  is  zero  up  to  round-off  errors.  (If  you  try  to 
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raproduca  thia  axampla  with  your  own  varalon  of  MICROSCOPE  you  may  gat  a 
alight ly  diffarant  raault  bacauaa  of  tha  random  noturo  of  round-off 
arrora.)  An  altarnativa  approach  conaiata  of  halving  tha  intarval,  and 
tharaby  dacraaaing  tha  ronga  eovarad  by  tha  graph.  In  tha  numarical 
diaploy  tha  ronga  ahould  ahrink  eorraapondingly ,  maintaining  tha  value  of 
zaro  in  ita  intarior.  Halving  tha  Intarval  ylalda  tha  diaploy: 


+ 


Point  ■  e.eeeeeeof^e  a  -  2.seeei>-03 

Olraetlon  -  1 . OeoeoeOf^O  h  «  1.S90e[>-«2 

Fe  (-2.830-82.  2.880-82) 

1/0:  29  8  8  1  2  3  27  28  29  38  NRML  on  currant  CALLS  -  287 

FIgura  3 


Tha  ronga  of  tha  arror  hoo  baan  raducad  to  from  -8.8283  to  8.8288 
and  0  zaro  volua  of  tha  arror  at  tha  point  of  axamlnotion  la  atill 
conalatant  with  tha  graphical  diaploy.  ThIa  procadura  con  ba  continuad. 
Procaading  aimilarly  at  tha  point  x  «  1  yialda  a  eimilar  raault.  and  you 
might  try  thia  for  youraalf.  At  tha  point  x  »  2.  howavar,  wa  obtain  tha 
diaploy; 
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Point  -  2.eeeeeeof«e  a  -  s.eeeeo-e^ 

Diraction  -  i.eeeeaeo+M  h  •  3.eeeeo-e2 

re  (  2.640-11.  1.420-62) 

I/O;  2S  6  6  1  2  3  27  28  29  36  NRML  on  currant  CALLS  -  362 

Figura  4  Intarpolot ion  ot  x  «  2 


Tha  graph  of  tha  arror  function  touchaa,  but  doaa  not  croaa,  tha 
X  oxia,  which  ia  eonaiatant  with  our  intarpoloting  to  tha  firat  darivotiva 
oa  wall  oa  poaition.  Tha  volua  of  tha  arror  ot  tha  cantar  eon  ba  raod 
diraetly  from  tha  numaricol  diaploy  ond  ia  -2.40-11  (you  nay  gat  o  alightly 
diffarant  volua  bacouaa  of  t*'*  rondom  round-off  arrora).  Thia  confirna 
intarpolot  ion  to  poait.wi..  Intarpolot  ion  to  tha  darivotiva  ia  likaly 
bacouaa  tha  graph  oppaora  to  hova  o  horizontol  tongant.  Howavar, 
confidanea  con  ba  roiaad  by  computing  tha  firat  darivotiva  of  tha  arror  oa 
wall,  yialding  tha  diaploy: 
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111 
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1111111111111111 


Point  -  z.eaeeeeiHeo  e  -  s.eeeeo-e^ 

Direction  -  1.0808000400  h  «  3.0OOOD-02 

F0  (  2.64D-11.  1.42D-02)  FI  (-6. 110-02,  1. 770-01) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  current  CALLS  -  374 

Figure  5  Interpolot ion  to  the  First  Oerivotive 


Inquiring  about  the  value  of  the  firet  derivotive  ot  the  center 
yields  the  number  50-4.  This  is  much  larger  then  the  volue  for  the  error 
itself,  and  hence  might  couse  some  doubts.  However,  the  derivotive  is 
calculated  numerically  and  hence  not  exoctly.  The  volue  reported  by 
MICROSCOPE  is  oetually  the  truncation  error.  This  can  be  confirmed  by 
halving  the  discretization  parameter  yielding  the  display: 
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11 
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111 


Point  -  e.eeeeeeo+ee  •  -  2.5eeeD-«3 

Oiroetion  -  1.8808800488  h  -  1.S080D-82 

F0  (  8.830-81,  1.120480)  FI  (  1.260408,  1.270408)  F2  (-2.S30-81,  2.530-81) 
I/O:  25  6  6  1  2  3  27  28  28  38  NRML  on  currant  CALLS  -  585 

Figure  7  Froa  End  Condition 


That  tha  aaeond  darivativa  la  indaad  zero  eon  new  bo  varifiad 
uaing  tha  aama  taehniquaa  ea  obova  (i.a.  inquiring  directly  or  ueing  o 
decreaead  diaeretization  paromotar). 


2.2  .SaQ8lhn888 

In  the  praeent  oxaiapla,  the  only  critical  act  ia  the  one 

containing  the  point  x  ■  1  (aaauming  we  are  confidant  that  S  hoe  boon 
progronoaod  to  be  polynomiol  olaewharo).  When  Inveet  igot  ing  oiaoothnaea  it 
ia  uaually  prafarobie  to  axonina  on  error  (if  o  primitiva  function  ia 
available)  rother  than  the  interpolont  itaelf.  In  the  preaent  exonpla,  the 
error  function  at  the  aingle  criticol  point  yielda  the  diaploy: 
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Point  ■ 

1.000009D400  a 

2.50000-03 

Direction  ■ 

1.00000OD400  h 

• 

1 . 50000-02 

Fe  (  2.490480.  2.970400)  FI  (  2.390400,  2.900400)  F2  (  2.480400,  3.110400) 
F3  (  2.730400,  4.060400) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  currant  CALLS  *  592 


Figure  8  Soothnaaa  ot  *  ■  1 


Thia  diaplay  can  ba  intarpratad  ainilorly  aa  figure  1.  and 
clearly  ahowa  a  function  that  ia  twice  but  not  three  tinea  di f ferent ioble. 
Thia  ia  what  wa  expect  fron  aur  conatruct ion. 


2.3  Degree  of  Praciaion 

The  ‘praciaion  cloaa'  of  an  approxinot ion  achene  ia  the  aet  of 
functiona  that  are  reproduced  exoctly  by  the  achene.  Of  particular 
intereat  are  the  polynoniala  in  the  praciaion  cloaa.  The  largaat  nunber  g, 
aoy,  auch  that  oil  polynoniala  of  degree  up  to  q  ora  in  the  praciaion  cloaa 
ia  the  degree  of  polynoniol  praciaion. 

To  decide  if  any  particulor  function  ia  in  the  praciaion  cloaa 
one  approxinataa  that  function  and  inveatigotea  if  the  error  ia  zero  within 
round-off  of  facta.  Thia  oppliea  oven  to  non-polynomial  functiona. 

To  onawer  the  more  narrow  queation  of  polynomial  praciaion  it  ia 
poaaibla  to  proceed  more  ayatemot ical ly  and  to  examine  polynoniala  of 
increaaing  degree.  The  polynomiola  ahould  be  oa  ganarol  oa  poaaibla,  i.e. 
they  ahould  be  unlikely  to  introduce  mialeoding  ortifocta.  One  poaaibility 
ia  to  chooae  the  coafficianta  randomly.  If  really  in  doubt,  one  might 
examine  each  polynomial  in  a  auitable  boaia  of  the  oppropriote  polynomiol 
apoca.  The  pointa  of  examination  ahould  ba  choaen  auch  that  eoch 
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appropriots  port  of  ths  domain  is  covsrsd.  For  instoncs,  in  ths  present 
sxompls  (which,  Ineidsntol ly ,  hos  linsor  prscision)  ens  should  chooss  ons 
point  >  1.  and  ons  <  1.  covsring  both  cosss  in  ths  dsflnition  of  S. 

Ths  following  figure  illustroiss  ths  typicol  oppsorones  of  o  zero 

error . 
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.  .  1 

.1  . 
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I  . 
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.  .  1 

I 

I 

.  I 

.  .  I. 

I 

1 

Point  -  2.3423880-81  s  -  2.58880-83 

Direction  «  1.8888880+88  h  >  1.58880-82 

F8  (-9. MO-11.  9.780^11) 

I/O:  25  6  6  1  2  3  27  28  29  38  NRML  or)  current  CALLS  -  667 

Figure  9  A  Zero  Error 


oL 


Notice  how  the  values  of  the  function  ore  spreod  randomly  over 
the  entire  display,  occurotely  pinpointing  the  round-off  nature  of  the 
error.  Sometimes,  however,  one  hos  o  situation  in  which  ths  coefficients 
of  the  interpoiont  ore  evaluated  inoccurotely,  leoding  to  on  error  thot  has 
0  scattered  distribution  with  a  bios  like  thot  in  the  following  display: 
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Point  -  2.3423080-91  e  -  2.50000-03 
Direction  ■  1 . 0000000400  h  •  1.50000-02 
F0  (  2.160-10.  7.310-10) 

I/O:  25  6  6  1  2  3  27  28  20  30  NRML  on  current  CALLS  «  742 
Figure  10  Round-off  With  o  Biae 


2.4  Degree  of  a  Poivnoeiiol 

Verifying  that  a  certoin  function  ie  o  polynomial  of  o  certoin 
degree  is  often  useful  when  examining  cross  derivotives.  For  example.  In 
the  bivariate  Clough-Tocher  scheme  (Atfeld,  1984)  which  ie  piecewise  cubic 
on  a  triangle,  differentiability  between  triangles  is  forced  by  the 
requirement  that  the  normal  derivatives  across  edges  be  linear  (rother  thon 
quodrotic)  along  the  edges.  In  this  illustration,  let  us  verify  that  our 
interpolont  is  indeed  cubic.  This  con  be  occomplished  by  showing  that  the 
third  derivative  is  constant  or  that  the  fourth  derivative  is  zero  (within 
round-off).  One  obtains  for  example  the  following  display: 
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Point  -  2.3423880-81  •  -  2.38880-83 

01  root  ion  -  1.8888880488  h  -  1.58880-82 

F8  (  1.180488,  1.430488)  FI  (  1.290488,  1.410488)  F2  (  3.870-81,  8.920-81) 
F3  (  2.730488,  2.730488) 

I/O;  25  8  6  1  2  3  27  28  29  38  NRML  on  currant  CALLS  «  841 

Fiquro  11  A  Cubic  Function 


Notico  thot  tho  third  dorivotivo  ic  cloorly  offoetod  by 
round-off,  without  ony  opporont  blot.  On  tho  other  bond,  tine*  tho 
ditployod  limitc  of  it«  rang*  or*  idonticol  tho  tizo  of  tho  rongo  io  iooo 
thon  1/273th  of  tho  dorivotivo  voluo.  A  ciooor  oxominotion  rovooio  thot 
tho  third  dorivotivo  rongoo  from  2.7386  to  2.7323.  A  tolltolo  oign  of 
round-off  orrora  ia  thot  thoy  con  bo  decraoood  ovon  furthor  by  Ine rooming 
tho  diacrot i zot ion  poromotor.  Doubling  h  yiolda  tho  diop  loy: 
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Oi  root  ion 

Fe  (  i.eeo+ee. 

F3  (  2.730486, 
1/0:  25  6  6 


-  2.3423000-81  s  -  S. 00060-03 

-  1.0000000400  h  «  3.00000-02 

1.560400)  FI  (  1.270400,  1.500400)  F2  (  1.340-01, 
2.730400) 

1  2  3  27  28  29  30  NRMl  on  curront  CALLS  -> 


1.150400) 

885 


Figoro  12 


Ths  range  is  now  from  2.731361  to  2.731569  which  is  e  reduction 
by  a  factor  8.  This  is  consistent  with  the  third  derivative  being  constant 
since  the  round-off  errors  in  the  third  derivative  increase  i i ke  the  third 
power  of  h,  see  Appendix  1.  Further  increases  do  not  olter  the  basic 
scattered  structure  of  the  third  derivative  (until  the  range  covered  in  the 
display  includes  the  point  x  ~  1,  where  the  third  derivative  has  o  step 
discontinuity).  Increasing  the  volue  of  h  should  eventually  decrease  the 
round-off  errors  below  a  level  at  which  the  true  change  in  the  derivotive 
becomes  visible  in  the  display  if  the  third  derivative  wos  not  constant. 
The  fact  that  this  does  not  happen  indicates  that  there  is  no  such  change. 
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3.  Uoino  MICROSCOPE 

This  toctlort  contain*  o  d**criptIon  of  th*  master  MICROSCOPE 
routine  MCRSCP  ond  o  first  Introduction  to  th*  command  longuog*  of 
MICROSCOPE.  Tor  a  more  detailed  description  consult  section  4. 

During  an  interactive  session  th*  terminal  screen  moy  be  either 
’active'  or  'Inactive'.  On  on  octiv*  screen,  chonges  In  parameters  or* 
Incorporated  immediatsly  into  the  olphonumericol  disploy  which  is 
maintoinsd  on  th*  screen.  On  an  Inoctiv*  screen,  conmond*  ere  occumuloted 
and  carried  out  only  when  th*  screen  i*  octivoted  by  either  th*  connond  GO 
or  FORCE.  An  inoctiv*  screen  shows  th*  lost  few  command*  and  th* 
corrsaponding  data  whereas  on  active  screen  show*  at  most  th*  currant 
command  and  Its  data,  in  addition  to  th*  olphonumericol  display.  (However, 
an  incomplet*  alphanumsrical  display  moy  b*  colled  on  on  inoctiv*  screen, 
see  th*  description  of  th*  command  RSCREEN  In  section  4.5  and  th* 
discussion  of  th*  alphanumerical  disploy  in  section  1.3.) 

A  second  distinction  applies  to  th*  woy  in  which  screen  chonges 
or*  implemented  on  on  octiv*  screen.  In  on*  version  of  MICROSCOPE,  which 
w*  coll  th*  'screen  version',  th*  olphonumericol  disploy  is  kept  stotionory 
ond  only  modified  rather  than  reploced.  This  version  requirss  you  or  your 
system  manager  to  provide  four  simple  screen  editing  routine*  thot  or* 
specific  to  each  computer/term  I  no  I  combinotlon  (see  chapter  S.)  If  this 
focility  is  not  avollobl*.  or  if  you  us*  o  hordcopy  tsrminol,  then  th* 
'scroll  version'  of  MICROSCOPE  mov  be  employed,  which  scroll*  th*  old 
screen  ond  reploces  It  by  a  new  version  of  ter  eoch  chong*. 

Th*  screen  version  is  much  preferobls,  ond  in  this  monuol  w* 
usually  ossum*  It  is  oval  lob  I*. 

Another  distinction  pertains  to  th*  utilization  of  lower  cos* 
letters.  Th*  us*  of  lower  cos*  letters  in  a  FORTRAN  program  is  not 
portable,  but  now  most  installation*  do  provide  this  focility.  In  this 
monuol  w*  ossum*  lower  cos*  letters  or*  ovoilabl*.  If  so,  any  lower  cos* 
letters  occurlng  in  th*  Input  to  th*  program  or*  treated  os  if  they  were 
upper  cos*.  Thu*  commands  eon  be  given  In  lower  cos*  letters.  Output  from 
MICROSCOPE  consist*  Of  0  mixture  of  lower  ond  upper  cos*  letters.  However, 
to  maintain  portability  w*  provide  o  version  of  MICROSCOPE  that  works 
solely  with  upper  cos*  letters.  Section  5  contains  more  information. 
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To  put  MICROSCOPE  into  operation,  you  will  hove  hove  to  write  o 
FORTRAN  moin  progrom  that  calls  the  routine 

MCRSCP(F.  INPUT.  OUTPUT.  GRAPHC,  HELP.  RECORD,  RESTART, 

X  LINES.  WIDTH.  PLOT.  NUMRCL.  PROMPT) 

The  first  parameter  is  the  name  of  a  DOUBLE  PRECISION  function 
(the  triol  function).  It  has  the  form 

DOUBLE  PRECISION  FUNCTION  F(X) 

DOUBLE  PRECISION  X 
DIMENSION  X(1) 


F  -  .  .  . 

RETURN 

END 


and  must  be  declared  EXTERNAL  in  the  colling  program.  MICROSCOPE  will  set 
at  most  the  first  three  components  of  X  before  colling  F  (i.e.  it  can 
handle  only  functions  of  up  to  three  variobles).  Notice  thot  F  does  not 
hove  any  parameters  other  than  the  values  of  the  independent  voriobles. 
Usually,  however,  it  will  also  depend  on  parameters  describing  e.g.  the 
domain  data  structure.  These  are  most  conveniently  possed  to  F  in  o  common 
block. 

All  of  the  remaining  parameters  of  MCRSCP  ore  input  integers. 
MCRSCP  will  not  change  them,  and  they  moy  therefore  be  integer  volues 
rather  than  variables.  The  first  six  of  them  (i.e.  INPUT.  OUTPUT,  CRAPHC. 
HELP,  RECORD,  RESTART)  are  FORTRAN  device  numbers  corresponding  to  the 
first  six  I/O  numbers  listed  in  the  numerical  display  (see  section  1.3). 
It  is  your  responsibility  to  set  up  the  necessary  correspondences  between 
device  numbers  and  physical  devices  or  disk  files. 

The  last  five  porometers  describe  the  physical  dimsnsions  of  the 
screen  display.  Their  meaning  is  os  follows:  (Numbers  in  porenthesss  give 
the  value  used  for  the  examples  in  this  monuol.  squore  brackets  contain  the 
allowable  ranges  for  the  parameters). 


LINES  (24): 

screen  [10  <  LINES  <  S8]. 


The  totol  number  of  lines  used  on  the  terminal 


WIDTH  (75):  The  number  of  columns  used  for  the  grophical  display 
(0  <  WIDTH  <  136].  Usually  one  will  choose  as  large  o  value  os  possible 
for  moximum  resolution.  However,  small  values  moy  sometimes  be  preferable 
for  increased  speed.  For  a  symmetricol ly  ploced  vertical  axis,  the  value 
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of  width  ahould  ba  odd.  Tha  nunaricol  dlaploy  occupies  up  to  75  colunns. 
This  ia  therefore  tha  nlnisiun  phyaicol  width  of  the  output  device.  WIDTH 
is  also  the  nuisber  of  points  ot  which  derivatives  ore  soispled  in  tha 
<PL0T79>  display. 

PLOT  (15):  The  nuisber  of  lines  el  located  to  the  graphical 
display  [8  <  PLOT  <  49]  Again,  a  value  of  PLOT  os  large  os  possible 
ahould  be  supplied.  The  nuoMrieol  effort  ie  independent  of  PLOT. 

NIMtCL  (8):  The  nusiber  of  lines  obove  the  screen  bottois  at  which 
the  nusMrIeal  disploy  begina.  [5  4  PROMPT  <  MAMCL  <  LINES  -  PLOT]. 

PROMPT  (2):  The  niesber  of  lines  obove  tha  scrs*n  botteei  et  which 
the  proMpts  ere  given.  {1  <  PRObTT  <  LINES  -  PLOT  -  8]. 

Ueuelly.  MCRSCP  will  be  colled  only  once  during  o  session. 
However,  it  is  peesible  to  coll  your  own  routines  fro*  MCRSCP.  You  hove  to 
supply  0  subroutine  SU8USR  without  any  poroweters  which  con  be  colled  by 
MCRSCP  A  detailed  discussion  of  this  procedure  is  given  in  section  4.15  on 
user  intervention.  For  new  we  way  oseuMe  that  SU8USR  is  o  dusssy  routine 

subroutine  SUBUSR 

RETURN 

END 

(which  is  supplied  with  the  MICWQ  '  pockoge  ond  which  you  eon  of  course 
write  yourself) 

When  called,  MCRSCP  first  .  vs  an  identifying  stencil  on  the 
OUTPUT  device,  containing  its  version  n  r  and  the  lost  SMdificetion 
dote.  Nest,  it  checks  to  see  if  v  •  screen  parameters  satisfy  the 
reetrlctions  listed  above.  If  not  then  MCRSCP  exits,  rather  then 
terminates,  in  erder  to  give  you  a  chance  to  correct  your  mistake. 

On  the  first  cell  to  MCRSCP  It  computes  the  swchine  round-off 
unit  (see  section  4.1)  end  prints  its  common  logarithm  to  give  you  en  idea 
of  how  many  decimal  digits  are  corrled  en  the  computer.  The  actual  number 
relevant  to  your  investigation  moy  of  course  be  less  depending  en  the 
detoils  of  the  arithmetic  involved  in  evaluating  the  trial  function. 

Also  on  the  first  coll  only.  MCRSCP  then  prints  the  first  line  of 
any  available  news  (see  section  4.18)  to  help  you  decide  if  you  wish  to 
reed  them. 

After  setting  some  defaults  MCRSCP  then  enters  the  command  mode, 
in  which  it  accepts,  interprets  end  carries  out  the  coessonds  entered  en  the 
Input  Device. 
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3.2  The  Command  Mode 

The  MICROSCOPE  prompt  ia  of  the  form; 

n  » 

where  n  is  the  number  of  the  command  to  be  given.  It  starts  at  1  and  is 
incremented  by  1  after  each  command. 

A  MICROSCOPE  command  ia  a  string  of  ot  leost  two  letters  ond 
digits.  However,  only  the  first  two  charocters  are  significant  and 

recognized  by  the  parser.  They  must  be  entered  in  the  first  two  columns 
immediotely  following  the  prompt.  On  on  active  screen,  this  will  be  in  the 
some  line  as  the  prompt,  on  an  inoctive  screen,  or  in  the  scrolling 

version,  on  the  next  line.  Blank  spoces  are  significont.  f MICROSCOPE 

simply  reads  the  command  using  a  2A1  formot.  The  remaining  columns  con  be 

filled  with  arbitrary  characters.  They  are  not  read  or  processed  by 

MICROSCOPE.  This  ia  sometimes  useful  for  including  comments  in  o  program. 
A  command  line  ia  terminated  by  pressing  the  RETURN  button. 

After  a  command  you  will  be  prompted  for  the  oppropriote  data  (if 
any  are  needed).  The  prompt  describes  what  dota  are  needed  and  ia  specific 
to  the  command.  There  ore  seven  types  of  doto  requests; 

-1-  No  Doto 

-2-  A  command  name  (to  obtain  help  on  the  requested  command) 

-3-  1  Integer 

-4-  2  Integers 

-5-  1  Double  Precision  variable 

-6-  1  Double  Precision  Vector 

-7-  2  Double  Precision  Vectors 

The  dimension  of  the  vectors  is  identical  to  the  number  of 
variables  of  the  trial  function.  (The  default  ia  2,  see  section  4.18  for  a 
listing  of  oil  defaults).  Numbers  ore  processed  by  on  internol  porser  and 
are  essentiolly  format  free.  The  following  distinctions  and  restrictions 
apply,  however: 

-  Blanks  are  ignored.  For  example  "1  2"  is  considered  the  single 
integer  12. 

-  Numbers  on  the  some  line  must  be  separated  by  commas. 

-  Numbers  can  also  be  seporoted  by  writing  them  on  seporate  lines 
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(omiting  tha  coma). 

-  If  two  vaetora  ora  raquaatad  tha  aacond  vactor  muat  atort  on  o 

naw  I i na  . 


-  An  aflipty  atring  (obtoinad  by  juat  puahing  tha  RETURN  button)  ia 
conaldarad  to  maan  zaro  whonavar  a  aingla  numbar  (Intagar  or  Raol)  la 
raquaatad.  Thia  ia  particularly  uaaful  whan  o  zaro  antry  impliaa  that  o 
aubffloda  (aea  tha  daacriptiona  of  tha  PLOT  ond  tha  USER  comanda  in  aactiona 
4.14  and  4.15  raapact i valy)  ia  to  ba  axitad. 

-  Whan  two  Intagara  ora  raquaatad  aithar  ona  or  both  con  ba 
rapraaantad  by  an  anpty  atring.  (Tha  only  comand  for  which  thIa  ia  uaaful 
ia  TYPE.) 

-  Floating  Point  numbara  (lika  1.0E-12  or  ID-4)  ora  procaaaad 

corract ly . 


-  Floating  Point  numbora  may  olso  ba  rapraaantad  oa  intagara. 

Aftar  racaiving  tha  data,  MICROSCOPE  proeaaaaa  tha  comand  and 
tha  data.  If  tha  acraan  ia  activa  tha  comond  ia  corriad  out  and  tha 
alphanufflorieol  diaploy  ia  updatad,  othorwiaa  tha  naxt  comand  ia  promptad 
for. 


If  a  comand  la  not  racognizad  o  itaaaoga  la  printed  to  that 
affact  and  a  naw  comand  ia  promptad  for.  If  data  ora  raquaatad,  ond  o 
numbar  ia  not  racognizad,  than  naw  dato  will  ba  raquaatad.  Howavar,  if 
data  ara  ayntoctical ly  correct,  but  othorwiaa  maaninglaaa  (a.g.  you 
raquaat  tha  groph  of  tha  117th  darlvotiva),  than  the  comond  will  ba 
aborted  and  a  now  comand  will  bo  prompted  for. 

There  are  currently  63  comanda  which  divide  into  11  major  groupa 
corraaponding  to  aactiona  4.5  through  4.15.  Rather  than  deacribing  tha 
groupa  in  thia  aection,  we  illuatrote  aome  of  them  in  tha  couraa  of  o 
aearch  for  a  point  of  reduced  amoothneaa  (the  'active  point'),  giving  tha 
nomea  of  the  comanda  uaed,  and  further  illuatroting  the  capabilitiaa  of 
MICROSCOPE. 
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3.3  A  Diaeonf  inui  tv  in  the  second  defiyfltiita 
The  trial  function  underlying  this  exomple  ia 

2  X  X 

f(x)  -  (eta*x  )/2  4  e  if  x  >  8  and  f(x)  >  e  otherwise 

(where  eta  •  -8.805).  Clearly  f  is  once  but  not  twice  differentiable  at 
the  origin.  The  exponentiol  tern  eervee  to  conceol  the  discontinuity  in 
the  second  derivative.  (Similar  functions  that  are  0  through  5  tisMS 
di f f erent iobi e  con  also  be  exomined  for  orbitrory  voluea  of  eto.)  The 
parameters  ore  chonged  by  giving  the  user  intervention  command  USER 
(assuming  the  test  package  has  been  looded.  see  oppendix  III  and  section 

5). 


Let  us  pretend  we  do  not  know  the  the  degree  of  smoothness  but  we 
wish  to  examine  it.  We  first  inform  MICROSCOPE  that  the  trial  function  has 
one  independent  variable,  using  the  DMNSN  commond.  Then  we  enter  a  seorch 
interval  (here  [-1.5],  say)  using  the  IINTVL  command.  At  this  stage,  the 
screen  is  still  inoctive.  This  is  changed,  ond  the  first  display  is 
obtained,  by  giving  the  GO  command.  The  display  will  remain  active 
throughout  the  remainder  of  this  section.  We  obtoin: 


I 


Point 
Oi  recti  on 
F0  (  3.680-81, 
I/O:  25  8  6 


-  1.8888880488  s  «  5.40540-82 

-  4.0800000480  h  -  3.24320-81 

2.010401) 

1  2  3  27  28  29  30  NRML  on  current 


CALLS  > 


960 


Figure  13 
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Thia  display  shows  no  lock  of  sffloothnsss  at  oil.  So,  not  knowing 
about  ths  natura  of  ths  suspsctsd  discontinuity  at  all  ws  turn  on  ths 
horiiontol  seals  (using  DSCALE)  to  locots  ony  phsnonsno  of  intsrsst  and 
activats  (using  OCRAPH)  ths  groph  of  ths  highsst  possibis  (i.s.  sikth) 
dsrivativs.  Thsss  ehangss  ors  incorporotsd  ons  by  ons  (not  shown  hsrs) 
into  an  octivs  serssn.  Ths  final  rssult  is  this: 


. 666668866 
666666666666666666 


66666666666666321 .987654666666666666666666666666123466789.123. .6789.1234567 


Point  -  1.0066000466  s  -  5.40540-62 

Oirsetion  -  4.0600000400  h  -  3.24320-61 

FO  (  3.680-01,  2.010401)  F6  (-6.070401.  6.270401) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  currant  CALLS  -  972 

Figurs  14 


Now  It  is  cisar  thot  thsrs  is  o  phsnonsnon  of  sons  kind  securing 
at  ths  point  -17  nsasursd  on  ths  horizontal  seals.  To  sxomins  this 
furthsr,  ws  shift  ths  graph  by  17  (using  SHIFT).  This  csntsrs  ths 
intsrssting  parts  of  ths  graph,  yislding: 
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16:1 
I  :  I 

16:1 
I  :  I 

1:61 
I  :  I 

16:61  ■ ■ 6666666666 
6666666666666666666666666666666654321 . 12366666666666666666666666689.1234567 

6:6  1 
I  1 

16  1  .... 

1  I  . 

I  6:  . 


I  6 


Point  -  8.1081080-02  e  w  5.46540-02 

Direction  -  4.0000000+00  h  *  3.24320-01 

F0  (  1.470-01.  8.000+00)  F6  (-6.070+01,  6.270+01) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  current  CALLS  -  989 

Figure  15 


Already  at  tfiie  etage  it  is  poeeible  to  tell  the  natur*  of  the 
reduced  emoothnees.  The  eixth  derivative  clearly  ehowe  four  distinct  local 
extrema.  This  qualifies  it  os  the  third  (numericol)  derivative  of  o  delta 
function.  The  second  derivotive  of  the  triol  function  must  therefore  be  a 
step  function.  The  trial  function  is  hence  once  but  not  twice 
di f ferent iobi e.  Derivatives  of  step  functions  ore  discussed  in  detoil  in 
section  4.2  and  in  appendix  I. 

To  confirm  our  diagnosis  directly  we  examine  the  second 
derivotive.  Turning  on  its  groph  ond  erosing  the  distrocting  graph  of  the 
trial  function  (using  EGRAPH) ,  however,  yields  the  following  disploy  which 
shows  no  lock  of  smoothness  ot  all  in  the  second  derivative: 
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226666666666 
12366666666666666666666666689 . 1234567 
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1  2222 
1  22222 
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Point  -  8.1661880-62  s  -  5.46540-62 

Oirection  -  4.6666660+66  h  -  3.24320-61 

F2  (  1.480-61,  8.680+66)  F6  (-6.670+61,  6.270+61) 

I/O:  25  6  6  1  2  3  27  26  29  36  NRMl  on  current  CALLS  -  989 

Figure  16 


The  second  derivative  appaara  smooth  becouse  the  discretizot ion 
is  too  crude  to  show  its  discontinuity.  So  we  keep  decreasing  the  volue  of 
the  di scret i zot ion  parameter  (and  simultaneously  the  intervol,  keeping  the 
window  width  fixed).  Since  the  active  point  is  only  approximately  in  the 
center  of  the  display,  this  will  cause  it  to  drift  to  one  side  which  hos  to 
be  counteracted  by  more  shift  commonds.  Because  of  the  orgonizotion  of 
MICROSCOPE  it  is  particularly  efficient  and  convenient  to  divide  (or 
multiply)  with  powers  of  two.  Halving  h  four  times  (using  HALVE)  followed 
by  a  shift  of  24  and  four  more  halvings  yields,  after  erosing  the  scole  and 
the  graph  of  the  sixth  derivotive: 
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Point  ■  -1.0408340-17  a  -  2.11150-04 

Oirection  -  4.0000000400  h  -  1.26690-03 

F2  (  9.920-01,  1.000400) 

I/O;  25  6  6  1  2  3  27  28  29  30  NRML  on  current  CALLS  -  1353 

Figura  17 


Now  it  ie  cleor  thot  there  2e  o  jump  diBeont Inuity  In  the  aecond 
derivative  of  the  trial  function.  The  location  of  the  octive  point  con  be 
aeen  to  be  0  4/-  about  a  ■  0.0002.  Even  the  eize  of  the  Jump  can  be 
eatimoted.  Inquiring  (uaing  TYPE)  for  the  voluea  of  the  derivotive  at  the 
left  and  right  end  pointa  of  the  window  (where  in  thia  crude  diaplay  it 
aaaumea  ita  extrema)  yielda  voluea  of  9.988460-1  and  9.965220-1 
reapectively,  correapondi ng  to  a  jump  of  -2.32D-3.  The  reol  jump  la  of 
courae  eto  •  -5.00-3.  The  diacreponcy  ia  duo  to  a  truncotion  error,  not 
juat  in  approximating  the  derivotive,  but  olao  in  picking  the  pointa  of 
evaluation.  The  eatimote  can  be  improved  by  ext ropo I ot i ng  (lineoriy)  the 
left  and  right  porta  of  the  aecond  derivative  to  the  point  of  exomination, 
uaing  on  eatimote  of  the  third  derivative  which  con  be  obtained  from  one  of 
the  earlier  diaploya.  Nith  the  third  derivotive  aet  equal  to  1,  thia 
procedure  yielda  on  eatimote  of  -2.32D-3  -  12a  ■  -4.860-3  which  ia  o  very 
good  reault  conaidering  the  difficultiea  of  numerical  differentiation. 

In  the  preaent  aetting,  it  ia  not  poaaibla  to  dacreoae  tha 
d i ac rot i zot i on  porometer  further  without  auccumbing  to  round-off  errora  in 
the  aecond  derivative  (you  moy  wiah  to  try  thia  youraelf).  However,  the 
obove  example  illuatrotea  how  to  alleviate  the  limitotiona  of  round-off 
errora  We  oireody  recognized  the  true  noture  of  the  reduced  amoothneaa 
when  axomining  tha  aixth  rather  than  the  aecond  derivative.  With  the 
accuracy  uaad  for  the  exomplea  in  thia  laonual  (10  digita),  atapa  in  the 
aecond  derivative  300  timea  aa  amol I  oa  the  one  above  con  be  detected.  The 
following  diaploy  ahowa  tha  groph  of  the  aixth  derivative  for  the  aoma 
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function  oa  obova  with  ato  -1.70-5  ond  h  ■  0.24.  Thaaa  voluaa  hava  baan 
tokan  frofli  table  4  In  oppandix  t  (tou  >  1. 00-10). 
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1  666 

I 

I  66 

I 

I  666 

I 
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1  8 

6  1  6666 

1 

6  666 

666  6  I 

666666666666  1  66 

6  I 

666666666666666666  1 

6  1 

Point  -  -1.0408340-17  a  -  4.00000-02 
OIraction  >  4.0000000400  h  •  2. 40600-01 
PS  (  2.230-01,  4.400400) 

I/O;  25  6  6  1  2  3  27  28  29  30  NRML  on  currant  CALLS  «  1440 
Figure  18  At  tha  limlta  of  raaolution 


Tha  aixth  darivotiva  ahowa  4  axtramo,  idantifying  tha  atap  in  tha 
aacond  darivotiva. 
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4.  A  Detailed  Description  of  Features  and  Commands 


4.1  Round-Off  Effects 

The  foct  thot  oil  orithmetic  operotions  ore  carried  out  with 
finite  accuracy  is  fundamental  and  oil  pervasive  in  working  with 
MICROSCOPE.  Numerical  Differentiation  is  porticularly  sensitive  to 
rounding  since  it  involves  the  subtraction  of  nearly  equol  quantities, 
leading  to  a  cancellation  of  eignificont  digits.  Consequentially,  the 
error  in  numerical  differentiation  does  not  decrease  indefinitely  with  the 
d i sc  ret i 2at i on  parameter,  but  rather  storts  to  behave  erratically  beyond  o 
certain  minimum  value  and  tends  to  i ncreose  as  the  discretization  porometer 
is  dee reosed  further.  However,  it  is  its  very  rondomneee  thot  mokes 
round-off  errors  readily  apparent  in  o  MICROSCOPE  disploy  so  that  you  will 
not  be  misled  into  occepting  meaningless  results  ot  their  foce  value. 

The  precise  effect  of  rounded  orithmetic  depends  on  several 

factors: 


-1-  The  number  of  significant  (binory)  digits  carried  by  a 
particular  computer.  Most  users  will  prefer  on  approximation  of  the 
equivalent  number  of  decimal  digits,  and  so  MICROSCOPE  computes  thot  number 
and  prints  It  out  at  the  beginning  of  o  session.  More  precisely, 
MICROSCOPE  computes  the  smallest  power  of  2,  tou,  soy,  such  that  within  the 
given  computer  arithmetic  the  quantity  1  -f  tou  satisfies  a  test  for  greater 
than  1,  ond  then  prints  the  rounded  value  of  the  negative  common  logarithm 
of  tou. 


-2-  The  accuracy  with  which  the  triol  function  is  being 
evaluated.  That  accuracy  will  always  be  less  thon  or  ot  most  equoi  to  the 
accuracy  indicated  in  -1-.  Usually  it  will  be  much  less,  for  any  of  the 
fol lowing  reasons: 

-2.1-  Typical  multivariate  interpolants  are  very  complicated 
algebraically,  involving  many  operations  which  tend  to  increase  round-off 
errors. 


-2.2-  If,  as  is  frequently  the  cose,  the  trial  function  is  the 
difference  between  a  primitive  function  and  its  interpolant,  then  a 
cancellation  of  significant  digits  occurs. 

-2.3-  If  the  disploy  function  is  a  cross  derivative  of  the  trial 
function  it  is  already  contaminated  by  round-off  errors  due  to  numericol 
differentiation,  thus  exhibiting  effects  that  one  would  normolly  expect  to 
occur  only  for  higher  order  tangentiol  derivatives. 
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-3-  The  precise  nature  of  round-off  effects  is  also  very 
sensitive  to  the  relative  magnitudes  of  vorious  tongential  derivotives 
involved.  For  exomple,  a  constont  function  will  exhibit  round-off  errors 
at  any  value  of  the  discretization  porometer. 

Of  the  above  three  sources  of  round-off  errors,  you  will  usually 
have  a  firm  knowledge  only  of  the  first.  Knowledge  about  the  other  two  may 
be  gained  by  using  MICROSCOPE  itself.  This  is  usuol ly  more  effective  than 
attempting  a  theoretical  analysis. 

The  precise  nature  of  round-off  errors  is  discussed  in  Appendix 
I.  Detailed  formulas  and  tables  ore  given.  However,  the  usual  approach 
will  be  to  experiment  with  various  values  of  the  discretization  porometer, 
choosing  it  as  small  os  possible  without  generating  the  telltale  rondom 
pattern  that  indicotes  round-off  errors. 


4.2  Derivatives  of  Step  Punetions 

One  perhops  surprising  feature  of  numerical  di f ferent lot  ion  ie 
that  doubtful  discontinuities  in  some  derivatives  con  often  be  displo/ed 
very  clearly  by  plotting  higher  order  derivotives,  even  though  o  lorge 
di ecret i zot ion  porometer  is  necessory  to  keep  them  free  of  round-off 
errors.  This  fact  becomes  quickly  apparent  empiricolly,  and  is  onolyzed  in 
Appendix  I. 


The  conceptual  basis  for  the  approach  is  provided  by  a  Diroc 
Delta  function  (see  Couront  ond  Hilbert,  1962,  p.456).  A  step  function  con 
be  thought  of  os  rising  at  the  discontinuity  with  infinite  slope  from  one 
discrete  value  to  another.  If  it  is  piecewise  constont  its  derivative  con 
be  thought  of  as  a  function  thot  is  zero  everywhere,  except  ot  the 
discontinuity,  where  it  is  infinite  such  thot  the  integral  over  the 
derivotive  equals  the  jump  in  the  step  function.  Thot  derivative  is  o 
delta  function. 

More  precisely,  o  delto  function  obout  o  point  P,  soy,  is  defined 
os  the  limit  of  a  sequence  of  non-negotive  functions  thot  ore  different 
from  zero  only  in  neighborhoods  of  P.  The  diometer  of  those  neighborhoods 
tends  to  zero,  and  the  value  of  the  integrol  over  the  functions  remains 
constant  (at  1  in  the  classical  definition). 

Derivatives  of  a  delta  function  can  be  visualized  by  considering 
o  function  in  the  obove  sequence.  Suppose  it  is  bell  shaped.  Then  ite 
derivative  will  have  two  local  extrema  of  opposite  signs.  The  second 
derivative  will  have  three  extrema,  etc.  Taken  to  the  limit,  the  k-th 
derivative  of  o  step  function  hoe  k  local  extremo  of  alternating  signs,  all 
located  at  the  some  point. 
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Numerical  derivatives  of  a  step  function  exhibit  the  some 
phenomena,  if  the  discontinuity  is  smeored  out  over  a  window  containing 
sufficiently  many  points.  This  is  illustroted  in  the  pictures  below.  The 

step  discontinuities  occur  in  the  m-th  derivotives,  where  m  w  0,1 . S. 

The  discretization  parameter  is  h  •  1,  ond  the  size  of  the  jump  is  1.  More 
precisely,  the  m-th  function  is  defined  by 


m 

s(x)  ai  6  if  X  <  0  and  x  /m!  otherwise, 
where  the  axclomotion  mork  denotes  the  foctoriol. 


Notice  that  os  the  degree  of  the  discontinuous  derivative 
increases,  the  smoothness  of  the  graphs  of  higher  d*rivativ0t  increases. 
This  is  consistent  with  the  interpretation  of  numerical  differentiation  as 
a  smoothing  process.  (The  numerical  derivatives  have  the  some  degree  of 
smoothness  as  that  of  the  display  function.  Just  the  number  of  the  points 
of  reduced  smoothness  is  increased.) 
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Direction  •  1.6660660466  h  -  1.66660466 
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F6  (-7.290463,  7.290403) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  current  CALLS  «  1662 


Figure  19  Derivatives  of  a  Step  Function  (m  m  0) 
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01  root  ion  ■  1.6666660466  h  «  1.66660466 

n  (-1.660466,  4.280-11)  F2  (-1.660466,  2.850-16)  F3  (-2.660466,  2.660466) 
r4  (-8.660466,  1.660461)  r5  (-8.160461,  8.160461)  r6  (-1.460463,  9.720462) 
I/O:  25  6  6  1  2  3  27  28  29  36  NRML  on  currant  CALLS  >  1797 


Figure  26  Darivativaa  of  a  C  Function  (m  •  1) 
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3333333333333 
544  44  33  I 

45  4  33  I 

4  5  44  3  I 

45  644  I 
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4  63  6666  I 

62  563  666666666866 

22  563  555555555 
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222  63  55  I 

6  2625  55  I 

33  222  5  I 

6  33  6  55222225  I 
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55555  4444444  33333333366666 


55555222222222222222222 


Point  -  e.eeeeeeotee  a  -  3.33330-62 

Direction  -  1.6666690406  h  -  1.06060466 

F2  (-1.660+66,  1.990-16)  F3  (-1.660+66,  1.740-69)  F4  (-2.670+60,  2.670+06) 
F5  (-8.980+66,  1.350+61)  F6  (-1.460+62,  1.460+62) 

I/O:  25  6  6  1  2  3  27  28  29  36  NRML  on  currant  CALLS  -  1932 


1 

Figure  21  Darivativaa  of  a  C  Function  (m  >  2) 
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I  55  465 

6666666666666666  5  65 

I  66  5  45 

I  66  55  6  45 

I  6  5  6  55  4 


33 


33  6  4 
336 
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44  6333 
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4444444 


55  I 

5  I 

555  I 

55555555555555 
I 
I 

6666666666666666 
66  I 

66  I 

6  33336  I 

666666333333333333333333333 


Point  -  0.6606860406  a  -  3.33330-02 

Direction  m  1.0006600400  h  •  1.06660400 

r3  (-1.660466,  1.600-69)  F4  (-1.330406,  1.140-08)  F5  (-3.000400,  3.000400) 
F6  (-1.280401,  2.700401) 

1/0:  25  6  6  1  2  3  27  28  29  36  NRML  on  current  CALLS  *  2067 


Figure  22  Derivotivea  of  a  C  Function  (m  -  3) 
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6  6  55 
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65 
566 
5  6 

5  666 

5  66666666666666666 


44  5 

45 

55  444 
5  444 

444444444444444444444444 


Point  •  0.0000000400  a  -  3.33330-62 

Direction  m  1.6000000400  h  •  1.06000400 

F4  (-1.000406,  9.600-69)  F5  (-1.380460,  6.610-68)  F6  (-4.130400,  4.130406) 
I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  current  CALLS  -  2262 
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Figure  23  Derivotivea  of  a  C  Function  (at  ■  4) 
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6666666666666666666668555 

66  555 
66  55 
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6  5 
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666666666666666666666 


66 


66 


55  6 

5  6 
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655 

6  55 

6  555 


66666  5555555555555555555555555 


Point  -  e.eeeeeeo^e  *  >  3.33330-82 
Direction  •  1.00060804^8  h  1.00000400 
P5  (-1.000400,  5.090-06)  F6  (-1.650400,  2.670-06) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  current  CALLS  -  2337 


4 

Figure  24  Derivatives  of  a  C  Function  (m  >  5) 


Ui _ Ihr.  IttJjuIaK 

Sometime*  it  is  desirable  to  use  very  smoll  window  widths  or  even 
completely  separate  differentiation  stencils,  for  exompls  when  seorching  o 
large  domain  for  a  point  of  reduced  smoothness.  Thsn  steps  in  o  dsrivotive 
will  appear  as  such  rather  than  os  groduol  slopes  that  might  be  obscured  by 
other  features  of  the  relevant  graphs.  Also,  cusps  and  pole*  ore  best 
disploysd  using  small  window*  (so*  appendix  III  for  some  exomples). 
Usually,  however,  it  will  be  preferable  to  use  lorger  window*. 

Thsre  or*  two  main  reason*  to  us*  overlapping  differentiation 
stencils:  first,  as  illustrated  in  the  preceding  subsection,  this  is 
necessary  to  resolve  all  the  extrema  of  o  derivotiv*  of  a  delta  function, 
and,  second,  it  ollow*  the  us*  of  on*  svoluotion  of  a  trial  function  for 
sevsral  values  of  a  dsrivative,  thereby  increasing  sfficisncy  and 
aecsisrating  your  intsraction  with  MICROSCOPE.  Using  overlapping  stencils 
is  equivolsnt  to  having  a  window  width  (w  >  2h/s)  that  is  greotsr  than  1. 

Du*  to  the  internal  organization  of  MICROSCOPE,  for  any  given 
window  width,  the  number  of  evaluotion*  (of  the  display  function)  required 
is  determined  solely  by  the  largest  degree  of  the  plotted  derivatives.  The 
following  table  list*  formula*  for  computing  the  number  N  of  evaluotion* 


A  Detailsd  Description  of  Feature*  and  Commands 
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required  given  the  window  width  (w).  the  highest  degree  of  a  plotted 
derivative  (k)  and  the  number  of  columns  (C)  in  the  grophical  display. 


mox  k: 

8 

1  or 

2 

3  or 

4 

5  or 

6 

w-2h/s 

<  1 

C 

3C 

5C 

9C 

2 

C 

c 

+ 

2 

2C  + 

3 

4C  4 

5 

4 

c 

c 

+ 

4 

C  + 

4 

3C  4 

8 

6 

c 

c 

+ 

6 

2C  + 

9 

2C  4 

9 

6 

c 

c 

+ 

8 

C  + 

9 

3C  4 

16 

18 

c 

c 

+ 

18 

2C  + 

15 

4C  4 

25 

12 

c 

c 

-f 

12 

C  4 

12 

C  4 

12 

14 

c 

c 

+ 

14 

2C  4 

21 

4C  4 

35 

16 

c 

c 

+ 

16 

C  4 

16 

3C  4 

32 

18 

c 

c 

+ 

18 

2C  4 

27 

2C  4 

27 

28 

c 

c 

+ 

28 

C  4 

28 

3C  4 

48 

22 

c 

c 

22 

2C  4 

33 

4C  4 

55 

24 

c 

c 

+ 

24 

C  4 

24 

C  4 

24 

Table;  Efficiency  of  Window  Widths 


The  special  role  played  by  w  •  12  is  obvious;  the  numericol 
effort  is  independent  of  the  degree  of  the  plotted  derivatives,  and  exceeds 
the  number  of  columns  in  the  disploy  only  by  o  constant.  This  is  true 
whenever  w  is  a  multiple  of  12.  The  reoson  for  this  is,  of  course,  that 
the  display  function  is  sampled  at  multiples  of  1/2  and  1/3  of  the 
discretization  porameter.  For  efficiency,  the  defoult  value  of  w  is  12. 
This  is  also  sufficient  for  most  exominotions  involving  delta  functions. 
Sometimes,  however,  o  lorger  multiple,  like  24  or  36  may  be  desirable. 
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LA _ Roforoneo  Tablo  of  Coamondo 


This  •ubaoction  contain*  a  quick  rof*r*nc*  tobl*  of  oil  ovollobl* 
conmond*.  For  ooch  cofnnond,  o  on*-lino  doocription  I*  givtn.  Ofton,  th* 
information  eontolnod  in  tho**  short  doacription*  will  b*  oil  thot  I* 
n**d*d.  Th*r*for*,  th*  d**cription*  or*  u*od  In  ••v*ral  oth*r  ploc**  An 
oppropriot*  ■•loction  of  tham  occur*  ot  th*  baginning  of  aoch  of  th* 
following  aubaaction*.  Th*  HSUMRY  commond  lists  tham  ordarad  by  topics. 
Th*  mor*  datollad  HELP  commond  givas  o  documantot ion  of  on  individuol 
command  haodad  by  th*  on*- 1  in*  dascription. 


Noma 


Dascription 


Data  Error  Scrasn 


ACCENT  mark  dor  I  vat iv*  with  ostarisk*  TOGGLE 
C1CROSS  raod  1st  componant  of  cross  direct  Ion 

C2CROSS  rood  2nd  componant  of  cross  direct  ion 

C3CROSS  raod  3rd  componant  of  cross  direction 
CCHANNL  Chong*  chonnal  number 
CDIRCTN  raod  direction  of  cross  direction 
CHVALUE  rood  ch  for  cross  di f farant lot  ion 
CORDER  read  order  of  croas  darivotlv* 

CWINOON  Chong*  window  width  by  factor  SIGN 
DICHOIR  raod  1st  componant  of  dir.  of  Inv. 

02CHDIR  raod  2nd  component  of  dir.  of  inv. 

D3CHDIR  raod  3rd  componant  of  dir.  of  inv. 

DCENTER  draw  canter  of  groph.  display  TOGGLE 
OGRAPH  draw  graph 

DIVIDE  divid*  h  by  intagar  m  SIGN 

OkMSN  read  dimension  v  of  domoin 
DOUBLE  doubi*  h 

DSCALE  draw  scol*  In  graphical  display  TOGGLE 

DXAXIS  draw  horizontal  axis  TOGGLE 

EGRAPH  eras*  graph 

EXIT  return  to  colling  program 

FLIP  rap  lac*  d  by  -d  TOGGLE 

FORCE  racofflput*  everything 

GO  compute  whatever  is  naadad 

HALVE  halva  h 

HELP  obtain  datollad  help  on  commond  cmd 
HSUMRY  obtain  a  summary  of  oil  commands 
IDIRCTN  rood  th*  direction  of  investigation 
IHVALUE  read  value  of  h 

lINTVL  rood  endpoints  of  line  of  invest ig. 
IPOINT  read  th*  new  point  of  examination 
LCROSS  lood  0  now  cross  direction 
LDRCTN  load  a  new  direction  of  Investigation 
LIST  print  all  currently  avoilobl*  command* 
LOG  give  channel  number  for  logging 
LPOINT  load  o  new  point  of  examination 


k 

-1<k<7 

immadiot* 

Cl 

cxe 

daactiv. 

C2 

cxe:i 

daactiv, 

C3 

Cxe;l 

daactiv. 

ic.nc 

4 

— 

C 

cxe 

daactiv. 

ch 

5 

daactiv. 

k 

-1<k<7 

daoctiv. 

m 

immadiat* 

DI 

Dxe 

daoctiv. 

D2 

Dxe:i 

daactiv. 

03 

Dxe;1 

daoctiv. 

— 

— 

Immadiot* 

k 

-1<k<7 

immadiat* 

m 

dxe 

immadiot* 

V 

6<v<4 

daoctiv. 

— 

— 

immadiat* 

— 

— 

immadiat* 

_ 

— 

Immadiot* 

k 

-1<k<7 

immadiot* 

— 

— 

daactiv. 

— 

— 

immadiot* 

— 

— 

act  i  V 

— 

— 

oct  i  V 

— 

— 

immadiot* 

cmd 

2 

— 

D  Dxe  daactiv. 


h 

h>6 

daoc  t i V 

A.B 

AXB 

daoct i V 

P 

— 

daoct i V 

kc 

3 

daactiv 

kd 

3 

daoct  iv 

nc 

kp 

3 

daoct i V 

A  Data! lad  Dascription  of  Feoturas  and  Commands 

4.3  Rafaranc*  Table  of  Command* 
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MULTPLY  multiply  h  by  factor  m  SIGN 

NEWS  print  the  news 

NORMAL  turn  normalization  of  D  on/off  TOGGLE 

OUTPUT  write  ecreen  immoge  on  recording  device 

P1CHPNT  read  first  component  of  P 

P2CHPNT  rood  second  component  of  P 

P3CHPNT  rood  third  component  of  P 

PAUSE  pause  until  RETURN  from  channel  nc 

PLOT  enter  plotting  (<PL0T79>}  mode 

QUIT  terminate  MICROSCOPE  session 

RCROSS  read  channel  number  for  crose  I  let 

ROIRCTN  read  channel  number  for  direction  iist 

RESTART  reestablish  earlier  parameter  eettings 

ROTATE  change  di recti  an  of  i nvest igot ion 

RPOINT  read  channel  number  for  point  list 

RSCREEN  refresh  screen 

RWIND  rewind  device  nc 

SETDF  reset  defaults 

SHIFT  shift  point  of  examination 

STORE  store  current  parameter  settings 

TCENTCR  type  value  of  k-th  derivative  ot  center 

TNOTE  take  a  note  on  recording  device 

TYPE  Type  value  of  derivative  at  some  point 

UNDO  undo  previous  parameter  changes 

USER  enter  user  mode 

WAIT  make  program  wait  for  GO  or  FORCE  com. 
ZOOM  change  window  width  by  foetor  SIGN 


mxe  immediote 

-  deactiv. 

-  daoctiv. 

1  deoctiv. 

1  deactiv. 


Dx0  deactiv. 

-  pending 

-  deactiv. 

-  immediate 


-1<k<7  immediote 


-  deoctiv. 

mxe  immediate 


—  Description  — 

TOGGLE  the  command  reverses  itself  if  called  again 

SIGN  means  that  reversing  the  sign  of  the  input  porometer 

generates  the  opposite  effect  (e.g.  multiply  insteod  of 
divide). 

—  Data  — 

A  one  endpoint  of  the  line  of  investigation 

B  the  other  endpoint  of  the  line  of  investigation  (P  •  (A  -f 

B)/2) 

C  cross  di rect ion 

Cl  first  component  of  cross  direction 

C2  second  component  of  crose  direction 


A  Detailed  Description  of  Features  and  Commands 
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third  coniponont  of  croo*  di  root  ion 
diicrat iiot ion  poramotor  for  croat  darivotivea 
confflond  noma 

diraction  of  I nvaat I  got  I  on 

firat  coffiponant  of  diraction  of  i nvaat i got i on 
aacond  componant  of  diraction  of  i nvaat i got i on 
third  coffiponant  of  diraction  of  invaatigot ion 
corractlon  to  ba  oddad  to  0 

typa  of  davica  (1  •  input,  2  ■  output,  3  •  graphic,  4  > 

racord,  5  ■  raatort) 

Intagar  ahift  in  numbar  of  columna  on  diaploy,  -fva  ahifta 
graph  to  right,  point  to  laft,  -va  la  ravaraa 
dagraa  of  o  darivotiva  (to  ba  plottad,  araaad,  ate.) 
poaition  of  tha  croaa  diraction  in  a  liat 
poaltion  of  tha  diraction  of  i nvaat igot ion  in  o  liat 
poaition  of  tha  point  of  axominotion  in  a  liat 
Intagar  laul  t  Ipl  iar. 

Intagar  dividar 
channal  (davica)  nuffibar 
point  of  axoffiinotion 

firat  coffiponant  of  tha  point  of  axomination 
taeond  componant  of  tha  point  of  axominotion 
third  componant  of  tha  point  of  axominotion 
no  doto  ora  raguirad 

—  Error  — 

Raatrletiona  on  tha  doto  ora  liatad  in  thia  column.  Tha  antry 

" - "  maona  thot  no  chack  ia  moda  for  arrora.  Tha  aymbo'  "X"  maona  "not 

aquol".  If  tha  arror  conditiona  cannot  ba  givan  axplisitly  than  tha 
following  rafaranca  numbara  hova  baan  uaad: 

1.  Componanto  to  ba  chongad  muat  not  hova  on  indax  that  ia 
graotar  thon  tha  currant  dimanaion. 

2.  Tha  command  muat  ba  o  volid  MICROSCOPE  conxaond.  If  it  ia  not 
raeognizad  you  will  ba  prompted  for  another  coiamand  to  ba  helped  on. 
Typing  any  valid  command  wilt  gat  you  aoma  information  on  it  ond  will  than 
return  you  to  command  moda. 

3.  Thera  muat  ba  one  vector  in  lagol  formot  in  aoch  line  of  tha 
lift.  If  illegal  charoctara  are  encountered,  on  error  termination  of  tha 
conaaond  will  occur.  Alao,  tha  indax  number  muat  net  exceed  tha  number  of 
linaa  in  tha  liat.  If  It  doaa,  tha  program  will  raoch  tha  end  of  tha  file 
containing  tha  liat,  and  tha  praciaa  action  after  that  dapanda  upon  tha 
computing  eyatam. 


C3 

ch 

cmd 

D 

01 

02 

03 

00 

i  c 

ia 

k 

kc 

kd 

kp 

ffl 

n 

ne 

P 

PI 

P2 

P3 


Screen  — 
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the  coflMtond  octivotas  the  acraen 

the  cocmand  daactivatea  the  acraen  -  i.a.  the  GO  or  FORCE 
conmand  muat  be  given  before  the  conputotion  proceada.  For 
aoma  conaitonda,  there  ia  a  check  if  the  newly  input  datum  ie 
identical  with  the  old  one.  in  which  coae  the  acraen  remaina 
oct i va . 

If  the  acraen  ia  active  then  the  affacta  of  the  command  are 
incorporated  before  the  next  command  ia  requeated 

If  the  acraen  ia  inactive,  thia  command  will  put  the 
graphical  diaplay  on  the  acraen,  even  if  computotiona  are 
pending. 

The  graphical  and  numerical  dieplay,  and  the  acraen  atatua 
are  not  changed.  Notice,  however,  thot  the  octuol  acraen 
diaplay  may  be  deatroyed  by  error  meaaogea  or  help 
information.  The  next  diaplay  changing  command  (or  RSCREEN) 
will  then  refreah  the  acraen. 


4.S  Commanda  to  Control  the  Dieelav 

In  thia  and  the  following  aubaeetiona,  input  to  MICROSCOPE  ia  vio 
the  aaaigned  input  channel,  unleaa  otherwiae  atated.  Thia  will  uaually  be 
the  terminol.  but  may  also  be  a  file  for  o  botch  type  aeaaion. 


ACCENT 

mark  derivative  with  aeteriaka 

TOGGLE 

k 

-1<k<7 

immediate 

OCENTER 

draw  center  of  graph,  display 

TOGGLE 

— 

— 

immediate 

DGRAPH 

draw  graph 

k 

-1<k<7 

immediate 

OMNSN 

read  dimenaion  v  of  domain 

V 

e<v<4 

deoct iv. 

DSCALE 

draw  acole  in  graphical  diaplay 

TOGGLE 

— 

— 

immediate 

DXAXIS 

draw  horizontal  axia 

TOGGLE 

— 

— 

immediate 

EGRAPH 

eraae  graph 

k 

-1<k<7 

immediote 

FLIP 

replace  d  by  -d 

TOGGLE 

— 

— 

immediate 

RSCREEN 

refreah  acraen 

— 

— 

pending 

TCENTER 

type  value  of  k-th  derivative  at 

center 

k 

-1<k<7 

— 

TYPE 

type  value  of  deri votive  at  aoma 

1  point 

k,  ia 

-1<k<7 

immediate 

We  now  turn  to  the  deacriptiona  of  individuol  commonda.  Screen 
diaplay  commanda  control  whot  oppeara  on  the  terminal  acraen.  (For  the 
affacta  of  theae  commanda  on  the  <PL0T79>  drawing  aee  aection  4.14.)  There 
ore  several  aubgroupa  of  commands; 

-1-  Deciding  on  the  Tangential  Derivatives  to  be  Displayed. 


oct  IV 
deoct i V . 

immediate 

pending 
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By  default,  only  the  diaploy  function  ia  diaployad.  The  graph  of 
the  diaploy  function  itaalf  ia  morkad  by  parioda,  the  graph  of  ita  k-th 
tongantiol  derivative  (k  ~  1.2,  ....  ft)  by  the  digit  k.  Ueing  DGRAPH  or 
EGRAPH  raapact I va ly ,  you  can  turn  on  or  off  the  grapha  of  tha  zeroth 
through  aixth  darivotivaa.  It  ia  alao  poaaibla  to  diaploy  no  functiona  ot 
all.  In  any  coaa,  tha  diaploy  function  ia  avoluotad  only  ot  the  pointa 
where  it  ia  needed. 


-2-  Accentuating  a  Croph. 

If  you  wioh  to  anphaaiza  tha  graph  of  o  porticulor  derivative 
than  you  eon  plot  Ita  graph  ueing  oatariaka  rather  thon  o  period  or  tha 
digit  k.  Tha  coiMond  ACCENT  turne  thia  option  on  or  off.  Note  that  tha 
numerical  diaploy  contoina  no  direct  information  on  which  darivotivaa  hove 
bean  oceantuotad.  Alao,  if  you  occantuata  aavarol  darivotivaa  at  once 
their  grapha  cannot  be  diet inguiahad  from  each  other.  Tha  groph  of  on 
accented  derivative  ovarwritaa  onything  alee  in  tha  diaploy,  except  for  tha 
canter  mork.  Tha  ACCENT  command  alao  hoa  tha  affect  of  turning  tha  graph 
of  o  darlvotiva  on,  i.a.  it  doaa  not  need  to  be  preceded  by  o  DCRAM 
command.  By  default  no  derivative  ia  occantad. 

-5-  Drawing  a  Horizontal  Axia 

A  horizontal  oxia  con  be  drown  through  tha  canter  of  tha 
graphical  diaploy,  uaing  tha  commond  OXAXIS.  If  daairad,  o  ecola  con  be 
obtoinad  by  uaing  DSCALE.  Both  commonda  ora  on/off  awitehaa.  Tha  point  at 
which  tha  oxia  or  acola  Intaraacta  tha  central  column  of  tha  graphical 
diaploy  doaa  not  in  general  corraapond  to  the  value  zero.  (That  point  ia 
different  for  each  of  tha  tangential  darivotivaa  plotted,  and  moy  not  even 
be  contained  in  tha  graphicoi  diaploy.)  Tha  purpoaa  of  tha  oxia  ia  to  give 
help  in  examining  aymmatry  propertiaa,  ond  tha  acola  ia  uaaful  in  locoting 
pointa  of  intaraat  for  more  datoilad  axominotion.  Both  tha  acola  and  tha 
oxia  ora  off  by  default. 


-4-  Marking  the  Canter 

Tha  on/off  command  OCENTER  moy  ba  uaad  to  mark  tha  canter  of  tha 
graphical  diaploy  with  o  plua  aign.  Thia  ia  eomatimaa  uaaful  for 
confirming  aymmatry  propertiaa.  Tha  canter  mork  ia  off  by  default. 

-5-  Flipping  tha  Direction  of  Invaat igat ion 

The  FLIP  command  hoa  the  affect  of  raplocing  tha  direction  of 
invaatigotion  d  by  -d.  Thia  turna  tha  groph  of  any  darivotiva  about  tha 
vertical  axia,  and  it  alao  turna  tha  graph  of  any  odd  degree  derivative 
about  the  horizontal  oxia,  thereby  poaaibly  diaantongl ing  grapha  that  hove 
bean  overwritten  by  othara. 
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Coneidor  for  oxofflplo  tho  oxponontiol  function  and  ite  first 
dorlvotiv*.  Ordinorlly  tho  two  gropht  would  bo  idonticol,  yielding  o 
dieploy  in  which  tho  display  function  is  conplotoly  covoroo  by  tho  groph  of 
its  first  derivative: 


I 

I  111 

I 

I  111 

I 

I  1111 

1 

I  111 

I 

I  1111 

I 

I  1111 

I 

I  1111 

1 
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I 

111111  I 
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1111111  I 

I 

11111111  I 

I 

111111111  I 

I 

Point  -  6.6666660+66 

t 

1 . 66670-62 

Direction  -  1 . 6666660+66 

h 

m 

1.66660-61 

F«  (  5.460-61.  1.850+66)  FI  (  5.410-61,  1.860+66) 

I/O:  25  6  6  1  2  3  27  28  29  36  NRML  on  current  CALLS  -  2424 

Figure  25  An  Invisible  Croph 

After  using  the  FLIP  conimand  one  obtoins: 
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111  1  . 

1111  I  I 

111  1:1 
111  I  :  I 

111  I  :  I 


Point  -  e.eeeeeeo+ee  a  -  i.66670-e2 

Diraction  -  -1.0000000+00  h  -  1 .00000-01 


F0  (  5.400-01,  1.850+00)  FI  (-1.860+00.-5.410-01) 

I/O:  25  6  6  1  2  3  27  28  29  30  NRML  on  currant  CALLS  -  2424 

Figure  26  Tha  Diaontanglad  Gropha 


whara  tha  two  grapha  ora  now  claorfy  vi'albla. 

-6-  Obtaining  Numarieol  Valuaa 

Tha  olphanufflarieal  diaploy  doaa  not  provida  praciaa  numarieol 
voluaa.  Valuaa  of  any  computad  darivotiva  to  21  digita  may  ba  obtainad  by 
uaing  tha  TCENTER  or  tha  TYPE  commond.  TCENTER  providaa  valuaa  at  tha 
cantar  of  tha  diaploy,  TYPE  doaa  ao  for  any  point  in  tha  diaploy.  For 
TYPE,  maaaura  tha  diatonca  in  columna  form  tha  cantar  of  tha  graphical 
diaploy,  poaitiva  to  tha  right  and  nagotiva  to  tha  laft.  For  clarity,  if 
tha  acraan  varaion  of  MICROSCOPE  ia  amployad,  tha  column  in  which  tha  volua 
la  raod  ia  briefly  marked  with  axclomotion  morka.  Returning  to  command 
mode  raquiraa  typing  a  corrioga  return.  Thia  ia  naceaaary  ao  that  tha 
prompt  for  tha  next  command  doaa  not  ovarwrita  tha  numarieol  information 
before  it  con  ba  raod. 

Only  darivativa  valuaa  that  hove  bean  previoualy  computad  con  ba 
read.  Thia  includaa  oil  darivotivaa  through  [k]  where  k  ia  tha  highaat 
degree  of  o  darivativa  that  hoa  been  previoualy  diaployad  on  tha  currant 
active  acraan,  and  [k]  la  k+1  if  k  ia  odd  or  k  i f  k  ia  even. 
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-7-  Refreshing  the  screen 

The  conmand  RSCRCEN  hos  two  distinct  applications: 

On  an  active  screen,  it  eeryes  to  refresh  the  screen  image  after 
it  has  been  destroyed  e.g.  by  a  system  messoge  or  by  HELP  information. 

On  on  inactive  screen,  the  lost  few  commands  instead  of  the 
alphanumerical  display  are  shown.  In  that  situotion,  RSCREEN  calls  the 

display  onto  the  screen  without _ corrvinp  out  ony  eomputat ione .  Thus 

changed  porometers  may  be  shown  in  the  numerical  display  without  having 
been  incorporated  into  the  graphical  disploy.  This  condition  is  indicated 
by  the  flag  "GO  pndg"  in  the  numerical  dieploy.  (If  the  screen  is  octive 
then  the  flag  "current"  is  displayed  ineteod).  The  second  application  is 
useful  if  the  alphanumerical  disploy  has  to  be  viewed  in  order  to  decide  on 
the  next  command. 

When  entering  MCRSCP  the  screen  is  inoctive. 

-8-  Setting  the  Domain  Dimension 

The  OMNSN  command  informs  MICROSCOPE  of  the  number  of  independent 
voriables  (1,2  or  3).  It  is  listed  here  becouse  it  does  affect  the 
numerical  display.  However,  this  commond  is  usuolly  given  ot  most  once  at 
the  beginning  of  the  MICROSCOPE  eeseion.  The  default  dimension  is  2. 
DMNSN  also  resets  the  point  of  examinotion,  the  direction  of  investigation, 
the  cross  direction,  and  the  order  of  the  cross  derivative  to  their  defoult 
values. 


_ Commonds  to  Control  the  Point  of  Examinotion 


IINTVL 

reod 

endpoints  of  line  of  investig. 

A, 8 

AXB 

deoct iv. 

IPOINT 

read 

the  new  point  of  examination 

P 

— 

deoct i V. 

LPOINT 

load 

0  new  point  of  examination 

Irp 

3 

deoct iv. 

P1CHPNT 

read 

first  component  of  P 

PI 

— 

deoct i v . 

P2CHPNT 

read 

second  component  of  P 

P2 

1 

deoct i V. 

P3CHPNT 

reod 

third  component  of  P 

P3 

1 

deoct iv. 

RPOINT 

read 

channel  number  for  point  list 

ne 

— 

— 

SHIFT 

shift 

point  of  examination 

i  s 

— 

immediate 

The  simplest  way  of  defining  the  point  of  examinotion  is  to  read 
it  from  the  input  device  by  ueing  the  IPOINT  commond.  Individuol 
components  may  be  defined  similarly  by  using  the  commands  P1CHPNT,  P2CHPNT 
or  P3CHPNT. 

Sometimes  it  is  useful  to  generote  o  list  of  possible  points  of 
exominotion  before  the  MICROSCOPE  session,  ond  then  refer  to  them  by  their 
position  (i.e.  line  number)  in  the  list,  rather  then  by  typing  the 
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individual  components.  This  is  focilitoted  by  the  LPOINT  command,  which 
prompts  for  the  position  of  the  point  and  then  reods  it  from  o  device 
(usuolly  0  file)  with  o  device  number  that  has  been  set  by  the  RPOINT 
command  (the  default  is  zero).  All  point  loading  devices  must  contain  one 
vector  per  line. 

The  SHIFT  command  can  be  used  to  shift  the  graph  by  a  specified 
number  of  columns  to  the  right  (if  the  number  is  positive)  or  left  (it  is 
negative).  There  is  no  restriction  on  the  number  of  columns  by  which  to 
shift  the  grophs.  In  particular,  the  number  moy  exceed  the  number  of 
columns  in  the  display. 

It  is  olso  possible  to  define  the  port  of  the  domain  thot  is  to 
be  shown  in  the  display  by  defining  its  two  endpoints.  This  is 
occomplished  by  the  IINTVL  command,  which  will  prompt  for  the  left  and 
right  endpoint  (in  that  order)  of  the  grophicol  disploy.  These  hove  to  be 
seporated  by  at  least  one  line  feed.  Notice  that  IINTVL  defines 
simultaneously  the  point  of  examination,  the  direction  of  i nvest igot ion, 
the  discretization  parameter,  and  the  display  interval. 


4.7  Commonds  to  Control  the  Direction  of  Invest igot i on 


01CHDIR 

read 

1st  component 

of 

dir.  of  inv. 

DI 

DX9 

deoct i V. 

02CHDIR 

read 

2nd  component 

of 

dir.  of  inv. 

D2 

Dxe;1 

deact i V. 

03CH0IR 

reod 

3rd  component 

of 

dir.  of  inv. 

DJ 

0X9;  1 

deoct iv. 

IDIRCTN 

read 

the  di rect ion 

of 

invest igot ion 

D 

Dxe 

deoct iv. 

LORCTN 

lood 

0  new  di rect ion  < 

of  investigation 

kd 

3 

deoct  iv 

NORMAL 

turn 

normal izotion 

of 

0  on/off  TOGGLE 

— 

— 

deoct iv. 

ROIRCTN 

read 

channel  number 

for  direct  ion  list 

nc 

— 

— 

rotate 

change  di rect ion  of 

invest igot ion 

DD 

Dxe 

deoct i V. 

The  direction  of  invest igot ion  con  be  read  from  the  input  device 
by  using  the  commend  IDIRCTN.  Individual  components  can  be  chonged  by  the 
commonds  D1CHDIR,  D2CH0IR  or  D3CH0IR. 

Similarly  to  the  point  of  exominotion,  the  direction  of 
investigation  may  also  be  defined  by  giving  ite  position  in  o  list.  The 
command  LORCTN  loads  the  direction  from  the  list  (usually  a  file)  whose 
channel  number  has  been  defined  by  the  ROIRCTN  commond. 

A  new  direction  may  also  be  defined  by  rotating  the  old  one.  The 
command  ROTATE  will  read  the  difference  between  the  old  ond  the  new 
direction,  and  updote  the  direction.  This  is  sometimes  useful  for 
exomining  the  effects  of  deviating  slightly  from  the  previous  direction  of 
investigation. 

The  numericol  display  will  olwoys  contain  the  direction  of 
investigation  os  it  has  been  input  to  MICROSCOPE.  Internally  the  direction 
of  the  derivative  (i.e.  the  vector  d  in  the  di f ferent lot  ion  formulas)  is 
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usually  the  given  direction  divided  by  it*  Euclidean  length.  This  leads  to 
the  standard  definition  of  a  directional  derivative.  However,  sometime*  it 
is  more  useful  to  use  the  unnormalized  direction  for  the  differentiation 
(i.e.  to  hove  d  identical  to  the  displayed  numerical  value).  The 
normalization  can  be  turned  on  or  off  using  the  commond  NORMAL,  and  is  on 
by  default.  Normalization  change*  the  values  of  the  directional 
derivatives,  ond  the  domain  that  is  covered  by  the  graphical  display.  An 
identical  graphical  display  (but  with  different  numerical  values)  con  be 
obtained  by  dividing  the  display  interval  inetead  of  the  specified 
direction  of  investigation  by  the  some  foctor. 


4.8  Commands  to  Control  Aceurocv.  the  Window,  ond  Round-Off  Effects 


CWINOOW  change  window  width  by  factor 

SIGN 

m 

(*>0 

immodiot* 

DIVIDE 

divide  h  by  integer  m 

SIGN 

m 

dx0 

immediate 

DOUBLE 

double  h 

— 

— 

immediot* 

HALVE 

halve  h 

— 

— 

immediate 

IHVALUE 

read  value  of  h 

h 

h>0 

deoct i V . 

MULTPLY 

multiply  h  by  factor  m 

SIGN 

m 

mX0 

immediate 

ZOOM 

enlarge  central  portion  of  displ. 

SIGN 

m 

mx0 

immediate 

This  group  of  commands  controls  the  three  parameter*  h 
(discretization  parameter),  s  (the  disploy  intervol),  ond  w  (the  window 
width)  which  or*  related  by  w  ■  2h/s  (w  must  lie  in  the  interval  [1/16,576] 
Olid  can  ossum*  only  certain  rational  values).  Each  of  the  commond*  chonge* 
two  of  the  three  parameters,  leaving  the  third  constant,  occording  to  the 
following  toble: 


Parameter  left  conetant 

w  s  h 

DIVIDE  CWINDOW  ZOOM 

DOUBLE 

HALVE 

IHVALUE 

MULTPLY 

Table  of  Discretization  Control ing  Commands 

Most  frequently,  the  window  width  is  left  unchanged,  often 
throughout  the  entire  MICROSCOPE  session.  This  occounts  for  the  obundance 
of  commands  covering  that  case. 

The  command  IHVALUE  serves  to  define  on  entirely  new  value  of  the 
discretization  parameter  (and  the  display  interval). 

Often  you  will  want  to  modify  the  current  volu*  of  h,  e.g.  to 
get  rid  of  the  influence  of  round-off  errors,  or  perhaps  to  increase  the 
accuracy  of  numerical  values.  Due  to  the  orgonizotion  of  MICROSCOPE  it  is 
particularly  efficient  (i.e.  require*  few  evaluations  of  the  trial 
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function)  if  h  is  multiplied  or  divided  by  powers  of  2.  The  factor  '  is 
used  so  frequently  that  two  specific  commands  (HALVE  and  DOUBLE)  hove  been 
provided  for  this  case,  which  require  no  orguments. 

Multiplications  or  divisions  with  other  integer  foctors  can  be 
occomplished  by  using  the  commands  MULTPLY  and  DIVIDE  respectively.  It  is 
possible  to  reverse  the  meaning  of  the  commends  by  entering  negative 
factors  (e.g.  MULTPLY  followed  by  -3  is  equivalent  to  DIVIDE  followed  by 
■f3).  This  is  due  more  to  the  internol  workings  of  MICROSCOPE  than  to  ony 
attempt  to  provide  a  truely  meaningful  convenience. 

The  command  CWINOOW  chonge*  the  current  window  width  by  the 
specified  factor.  A  positive  factor  corresponds  to  on  increase  of  the 
window  width,  a  negative  one  to  o  decreose.  The  change  is  obtained  by 
altering  the  value  of  the  discretization  parameter  while  leaving  the 
display  intervol  fixed. 

The  command  ZOOM  ha*  a  similar  effect.  However,  it  olters  the 
disploy  interval  and  leaves  the  discretization  parameter  unchanged.  This 
is  useful  when  operating  close  to  the  round-off  level  with  o  value  of  h 
that  must  not  be  decreased.  On  the  other  hand,  if  round-off  error*  or* 
irrelevant  then  the  CWINDOW  command  should  be  used  because  it  is  more 
*f f icent . 

The  primary  purpose  of  the  command*  described  her*  (other  thon 
IHVALUC  la  to  alter  the  display  while  the  screen  is  octiv*.  If  severol 
command*  are  given  while  the  screen  is  inactive  only  the  last  on*  will 
actually  be  implemented.  If  severol  chonges  ore  to  be  effected  before  the 
triol  function  is  evaluated,  then  this  con  be  occomplished  by  erosing  the 
graphs  of  all  derivotives,  and  then  octivating  the  screen  (using  GO  or 
FORCE),  and  moking  the  desired  change*. 

For  the  following  illustrotion  it  is  ossumed  that  the  screen  is 
octiv  derivatives  through  second  order  hove  just  been  disployed,  79  columns 
are  being  carried  in  the  graphical  display,  end  the  current  window  width  is 
12.  The  table  below  give*  the  additional  number  of  evoluotions  of  the 

display  function  needed  to  accomodot*  changes  by  factors  2,3 . B.  An 

entry  " — "  indicates  that  the  change  is  impossible  (e.g.  repeoted 
opplicotions  of  the  HALVE  commond  can  only  occomplish  chonge*  by  powers  of 
2).  Entries  in  parentheses  indicate  thot  the  change  has  led  to  o 
non-integer  half  window  width  and  therefor*  is  unlikely  to  be  employed. 
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2  3  4 

5 

6 

7 

8 

DIVIDE 

46  91  68 

91 

76 

91 

80 

MULTPLY 

46  66  68 

91 

76 

91 

80 

DOUBLE 

46  —  92 

— 

— 

— 

138 

HALVE 

46  —  92 

— 

— 

— 

138 

IHVALUE 

91  91  91 

91 

91 

91 

19 

CWINDOW 

12  24  36 

48 

66 

163 

84 

(increosing) 

CWINOOW 

6  6  (82) 

(237) 

6 

(237) 

(237) 

(decreasing) 

ZOOM 

S2  76  96 

139 

151 

163 

154 

( i ncreasi ng) 

ZOOM 

40  52  (116) 

(182) 

66 

(198) 

(204) 

(decreosing) 

Table:  Relative  Efforts 

4.9  Commands  to 

Control 

Executi 

,911 

EXIT 

return  to  calling  program 

deact  i  V 

FORCE 

recompute  everything 

— 

— 

act  i  V 

GO 

compute  whatever  is  needed 

— 

— 

act  i  V 

QUIT 

terminate  MICROSCOPE  session 

— 

— 

— 

RESTART 

reestablish  eorlier  parameter 

sett ings 

— 

— 

— 

SETDF 

reset  defaults 

— 

— 

deact  i  V 

STORE 

store  current  parameter  settings 

— 

— 

— 

UNDO 

undo  previous  parameter  changes 

— 

— 

— 

WAIT 

make  program  wait  for  GO  or  FORCE  com. 

— 

— 

deac  t i v 

To  activate  an  inactive  screen  the  GO  or  the  FORCE  command  must 
be  used.  GO  assumes  that  previous  function  evoluotions  are  still  valid  and 
FORCE  reevaluates  the  trial  function  at  all  needed  points  regardless  of 
prior  evaluations.  Ordinarily  GO  wi I  I  be  ueed  for  the  sake  of  efficiency, 
but  the  FORCE  command  is  necessary  if  changes  hove  been  applied  to  the 
trial  function  (either  after  leaving  MCRSCP  or  by  using  the  USER  command) 
that  invalidate  the  earlier  evaluations. 
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Uauolly,  th«  conmond  QUIT  !•  usad  to  tarminota  a  MICROSCOPE 
aaaaion.  That  eoniniand  laada  to  a  aimpla  FORTRAN  STOP  atotamant.  Howavar, 
it  la  poaaibla  to  laova  tha  routina  MCRSCP  vio  tha  EXIT  command  (loading  to 
a  RETURN  atotamant),  and  than  coll  MCRSCP  ogoin,  parhopa  of  tor  changing 
coma  aeraan  poromatara.  Aftar  raontaring  MCRSCP,  thaaa  will  again  bo 
chaekad  for  conaiatoncy,  and  tha  aeraan  will  ba  morkod  inoctiva.  Howavar, 
tha  dafoulta  will  not  ba  aat,  ao  thot  it  ia  poaaibla  to  raturn  to  MCRSCP 
and  find  avarything  aa  it  haa  boon  laft  with  tha  axcaption  of  what  you  may 
hova  raaat  youraalf  batwaan  MCRSCP  colla.  Aftar  raontaring  MCRSCP  it  may 
ba  noeaaaary  to  uao  FORCE  inatood  of  GO  to  avoid  paculior  diaployo  duo  to 
oltarad  aeraan  poromatara. 

Tha  raeommandad  way  of  ehonging  problam  dafining  poromatara  (a.g. 
aalaeting  o  now  trial  funetion)  ia  to  uaa  tha  eommond  USER  (aaa  aaction 
4.13). 


Uauolly  tha  aeraan  ia  aat  inoetiva  by  aoma  droatie  chonga  of  kay 
poromatara  aueh  oo  tha  point  of  axomination  or  tha  diraction  of 
invaat igat ion.  It  eon  olao  ba  aat  inoctiva  by  uaing  tha  eommond  WAIT.  Thia 
option  ia  porticulorly  uaoful  in  tha  ocrolling  varaion  whan  it  may  ba 
eonvaniant  to  affoet  aavarol  ehongao  bafora  ocrolling  tha  diaploy  Juat 
onea,  rothar  than  hoving  it  aerollad  oftar  aoch  individual  command. 

Tha  atota  of  o  MICROSCOPE  Invaatigot ion  eon  ba  atorad  uaing  tha 
eommond  STORE,  and  lotar  raotorad  uaing  tha  eommond  RESTART.  STORE  writaa 
tha  voluaa  of  oli  ralavont  poromatara  onto  tha  Raatort  Daviea  ond  than 
rawinda  thot  davieo.  Thua,  if  tha  command  STORE  ia  givan  aavarol  timaa 
only  tha  atotua  during  tha  loot  iaaua  of  tha  command  con  ba  raatorad.  If 
it  ia  nacaaaory  to  atora  aavarol  doto  aattinga,  thon  tha  raatort  chonnal 
numbar  hoa  to  ba  chongad  by  uaing  tha  CCHANNL  eommond  (aaa  aaction  4.13), 
or  o  diffarant  fila  hoa  to  ba  oaaignad  to  It  by  tha  command  USER. 

If  daairod,  aoma  of  tha  dafoult  aattinga  con  ba  raaatobl iahad  by 
uaing  tha  SETDF  command.  Mora  praciaaly,  tha  affacta  of  SETDF  oro; 

-1-  Tha  aeraan  ia  daoctivotad 

-2-  Ail  pravioua  function  avoluotiona  ora  diacordad 

-3-  Tha  window  width  w  ia  aat  to  12 

-4-  Only  tha  diaploy  function  itaalf  ia  floggad  to  ba  plottad 

-3-  Any  croaa  darivotivaa  ora  turnad  off 

-6-  Tha  diaeratizotion  poromatar  h  ia  aat  to  ita  dafoult  volua  (10-4) 

-7-  Tha  x-oxia,  cantar  mark,  and  horizontal  acolo  ora  oil  turnad  off 
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-8-  If  the  test  package  hoa  been  loaded  (see  Appendix  III)  then 
the  test  function  is  reset  to  its  defaults  (N  •  1 ,  D  >  10. 
no  exponential  is  added). 

While  the  screen  is  inactive,  it  may  be  possible  to  recover  from 
on  erroneous  command  using  the  command  UNDO.  This  hos  on  effect  only  on 
quantities  sat  oftsr  the  lost  GO  or  FORCE  command,  ond  resets  them  to  their 
values  prior  to  that  command.  The  following  porometers  are  reset  (i.e. 
ony  change  is  undone); 

-1-  The  point  of  examination 

-2-  The  direction  of  investigation 

-3-  The  left  and  right  endpoints  of  the  investigation  interval 
-4—  The  value  of  the  di acrat i zot ion  porometer  h 
-5-  The  window  width  w 
-6-  The  number  of  independent  variables 
-7-  Marking  the  center  point 
-8-  Drawing  a  horizontal  scale 
-9-  Drawing  a  horizontol  ox  is 
-10-  Selection  of  graphs  to  be  plotted 

Note  that  any  changes  regarding  o  cross  deri votive  or  the 
norma H zot ion  of  the  direction  of  investigation  ore  not  undone. 


HELP  obtain  detailed  help  on  command  cmd  cmd  2 

HSUMRY  obtain  o  summary  of  all  commands  — 

LIST  print  oil  currently  available  commonds  -  — 

NEWS  print  the  news  — 


The  quickest  way  of  obtoining  information  on  possible  commonds  is 
the  command  LIST  which  will  print  o  liat  of  oil  currently  ovoilable 
commonds  on  the  screen.  More  detailed  informotion  on  oil  commands  can  be 
obtoined  through  the  HSUMRY  commend  which  prints  the  table  in  section  4.3 
reorganized  by  topics. 

Information  on  a  specific  command  con  be  obtained  by  typing  HELP, 
foi  lowed  by  the  name  of  the  command  you  wont  help  on.  MICRC^QPE  will 
print  not  only  a  short  description  of  the  commond,  but  also  the  currant 
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values  of  those  parameters  that  are  affected  by  the  command.  Thus  if.  for 
example,  you  wish  to  find  out  what  the  current  number  of  independent 
variables  is  you  can  request  HELP  on  DMNSN. 

The  command  NEWS  will  print  o  text  thot  may  chonge  from  time  to 
time  and  will  also  depend  on  your  porticulor  installation.  At  the  time  of 
this  writing  the  "news"  ore  a  very  brief  primer  on  MICROSCOPE. 

Finally,  If  you  type  a  question  mork  to  MICROSCOPE.  it  will 
invite  you  to  enter  a  command,  ond  will  suggest  HSUMRY  ond  LIST  as 
possibi I ities. 


4.11  Commands  to  Keep  o  Record 

LOG  give  channel  number  for  logging  nc  -  - 

OUTPUT  write  screen  immage  on  recording  device  -  -  - 

TNOTE  toke  o  note  on  recording  device  text  -  - 

Results  obtained  with  MICROSCOPE  con  be  recorded  permanently  on 
the  Record  Device.  The  command  TNOTE  will  copy  any  text  you  type  on  the 
Input  Device  until  it  encounters  (he  string  'EC*  (End  Comments,  both  E  ond 
C  must  be  copitol)  in  columns  1  and  2.  It  then  returns  to  command  mode. 
The  length  of  text  lines  may  be  up  to  86  columns  or  the  number  of  columns 
specified  for  the  graphical  disploy.  whichever  is  larger.  Trailing  blanks 
are  stripped,  to  moke  the  record  file  more  compoct. 

The  command  OUTPUT  causes  MICROSCOPE  to  write  o  copy  of  the 
grophic  and  the  numerical  display  onto  the  recording  device.  This  monuol 
is  generated  by  using  TNOTE  and  OUTPUT  throughout  a  lengthy  string  of 
MICROSCOPE  commands. 


There  is  also  a  logging  device  which  is  active  if  its  chonnel 
number  is  different  from  taro  (It  Is  inoctive  by  defoult).  The  commond  LOG 
sets  the  logging  device  number.  If  the  device  is  active  all  commands, 
data,  disploye,  and  notes  ore  recorded  on  it. 
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Cl CROSS 

C2CR0SS 

C3CR0SS 

CDIRCTN 

CHVALUE 

CORDER 

LCROSS 

RCROSS 


read  1st  component  of  erase  direction 
read  2nd  component  of  cross  direction 
rood  3rd  component  of  cross  direction 
read  direction  of  cross  direction 
read  ch  for  cross  differentiation 
reod  order  of  cross  derivative 
load  0  new  cross  direction 
reod  channel  number  for  cross  list 


Cl 

cxe 

deoct  i  V 

C2 

Cxe;1 

deoct i V 

C3 

CX6;1 

deoct  i  V 

C 

CxB 

deoct  i  V 

ch 

5 

deoct  i  V 

k 

-1<k<7 

deoct  i  V 

kc 

3 

deoct  i  V 

nc 

— 

— 
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Cross  derivatives,  like  all  derivotives  in  MICROSCOPE,  are 
calculated  using  numerical  di f f erent iot ion.  The  formuloa  ora  onalogoua  to 
those  given  for  tongential  derivatives.  In  generol,  the  cross  stsncils 
will  hove  no  points  in  common.  Therefore,  one  evoluotion  of  o  k-th  order 
cross  derivative  requires  k-fl  evaluations  of  the  trial  function. 

By  choosing  the  cross  direction  to  be  porollel  to  the  direction 
of  investigation  the  range  of  tangential  derivotives  con  bs  shiftsd  from 
[0,6]  to  [k.S^k].  This  situation,  however,  is  not  treated  separately  by 
MICROSCOPE  and  does  not  lead  to  sovings  in  function  evoluotions. 

The  commands  controlling  the  cross  deri votive  (if  ony)  ore 
simllor  to  those  controling  the  direction  of  investigation.  By  default,  no 
cross  derivotivs  is  computed.  It  is  activated  by  reoding  o  nort-zero  order 
of  the  derivative,  using  COROER,  or  by  defining  a  cross  direction  using 
COIRCTN  (which  reads  the  cross  direction).  C1CR0SS.  C2CROSS,  or  C3CROSS 
(which  read  individual  components),  or  LCROSS  which  loods  o  cross  direction 
from  0  channel  defined  by  RCROSS  (the  default  channel  number  is  zero).  The 
dafoult  order  is  1  if  the  cross  derivative  is  activated  by  setting  a  cross 
direction.  If  it  is  activated  by  setting  the  order  its  default  direction 
is  the  vector  with  oil  of  its  components  aquol  to  1. 

The  cross  discretization  poromater  ch  is  either  set  explicitly 
using  CHVALUE,  or  by  default,  in  which  cose  it  ossumes  the  some  volue  os 
the  tangential  parameter  h  at  the  time  when  the  cross  derivative  is 
octivoted.  For  efficiency,  later  changes  of  h  hovs  no  effect  on  the  volue 
of  the  cross  parameter.  Note  that  ony  change  of  eh  raquiree  o  reevoluotion 
of  the  display  function  at  all  needed  points. 

Cross  derivatives  are  turned  off  by  setting  their  order  to  zero 
(using  COROER).  The  NORMAL  commend  hos  the  some  effect  on  cross 
derivatives  as  it  does  on  tangential  derivatives. 

When  plotting  cross  derivatives  it  is  importont  to  keep  in  mind 
that  the  significance  of  the  window  moy  be  blurred  if  o  stencil  used  for 
the  cross  different iot ion,  and  centered  outside  of  the  window,  nonetheless 
crosses  the  front.  This  con  be  avoided  either  by  keeping  the  cross 
discretization  parameter  much  smaller  than  the  display  interval,  or  by 
choosing  the  cross  direction  parallel  to  the  front. 

In  the  following  example,  the  trial  function  is 

2 

f(><.y)  “  *  •ob8(x)»y»obe(y) 

the  point  of  examinotion  is  [0,1],  the  direction  of  investigation  is  [l.B], 
ond  the  cross  direction  is  [1.-1].  The  value  of  h  •  3t>-3  and  ch  ■  6t>-3. 
Obviously  the  affects  of  the  discontinuity  in  the  second  tongential 
derivative  ore  carried  beyond  the  window. 
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Figure  27  An  Cxtandad  Oiacont inui ty 
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CCHANNL  Change  channel  number 

lc,nc  4 

— 

PAUSE 

pouae  until  RETURN  from  channel  nc 

nc  - 

— 

RWIND 

rewind  device  nc 

nc  — 

Soma  rudimantary  progronviiing  of  MICROSCOPE  ia  poeaibla.  Tha 
baaic  opprooch  conaiata  of  writing  o  aaquanca  of  comaonda  into  a  file  ond 
oaaigning  it  to  tha  input  device.  The  vorioua  davicaa  can  be  awitchad  by 
uaing  tha  CCHANNL  conmand.  It  will  prompt  for  two  intagera,  tha  firat 
defining  tha  type  of  channel  to  be  awitchad,  tha  aacond  giving  tha  new 
channel  number.  Tha  typaa  are 

1 :  Input  Device 
2:  Output  Device 
3:  Graph ica  Device 
4:  Recording  Device 
5;  Reatart  Device 
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For  on  axplonotion  of  theaa  davicaa  aea  aaction  1.3.  Thair 
channel  nuiabara  ora  defined  originally  by  poromatera  poaaad  on  to  the 
routine  MCRSCP  (aea  aaction  3.1).  Upon  entering  MCRSCP  theae  ora  copied 
into  vorioblaa  that  may  be  changed  by  CCHANNL.  Thua  the  input  poromatera  of 
MCRSCP  ore  not  changed  by  that  commend  ond  oaaume  their  old  voluaa  ahould 
MCRSCP  be  exited  (uaing  EXIT)  and  later  reentered  without  changing  ita 
poromatera 


The  loading  devicea  for  the  point  of  examination,  the  direction 
of  i nvaat i got i on,  and  the  croaa  direction  con  be  awitchad  by  uaing  the 
commonda  RPOINT,  RDIRCTN  and  RCROSS,  reapect i ve I y .  The  Help  Device  cannot 
be  switched  aa  the  informotion  contoined  in  it  is  on  integrol  port  of  the 
MICROSCOPE  pockoge  and  is  only  read  during  o  session.  Indeed,  it  is 
advisable  to  protect  the  help  file  ogainst  on  occidental  overwrite  if 
possible,  e.g.  by  uaing  a  READONLY  switch  in  on  OPEN  statement. 

Any  channel  con  be  rewound  by  the  RWIND  command.  This  facility 
con  be  used  e.g.  to  regenerate  a  specific  porameter  configuration  or  to 
overwrite  previously  generated  informotion  thot  is  no  longer  needed.  An 
infinite  loop  con  be  created  by  putting  a  rewind  command  ot  the  end  of  the 
input  file. 

The  PAUSE  command  can  be  used  to  view  the  output  of  o  MICROSCOPE 
program  at  leisure  during  run  time.  It  will  oak  lor  o  channel  number  and 
then  suspend  computation  until  o  dummy  character  is  received  from  thot 
chonnel.  The  channel  will  usually  be  the  terminol. 

Computation  by  a  program  con  be  suspended  by  changing  the  input 
device  to  the  terminol  (entering  interactive  mode)  ond  later  continued  by 
switching  the  input  channel  back  to  the  appropriate  device. 

Note  that  for  all  of  the  device  handling  commands  no  precaution 
has  been  token  to  ensure  that  only  logoi  channel  numbers  ore  given,  ond 
that  oppropriote  devices  are  assigned  to  the  channel  numbers. 

Comments  can  be  included  in  a  program  in  two  ways.  Firstly,  o 
’comment  line'  replaces  a  command  line  and  is  marked  by  o  "-f"  character  in 
column  1  or  2.  Secondly,  any  chorocters  contained  in  the  space  starting 
with  column  3  of  a  command  line  are  ignored  by  MICROSCOPE.  This  facility 
does  not  exist  for  lines  containing  data.  Indeed,  ony  spurious  chorocters 
may  lead  to  the  rejection  of  the  data  by  the  parser. 
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PLOT  enter  plotting  (<PL0T79>)  mode  -  -  - 

The  alphanumerical  dieplaye  ueuolly  generated  by  MICROSCOPE  ore 
adequate  to  obtain  the  desired  onewere  to  your  quest  lone,  but  they  moy  not 
be  suitable  e.g.  for  pubi icot ions,  reports,  etc.  Depending  upon  the 
availability  of  suitable  hard-  and  software  (<PL0T79>.  see  Beebe.  1979)  you 
may  be  able  to  generate  more  presentable  plots.  This  is  occomplished  by 
the  commend  PLOT  which  generates  a  drawing  of  the  current  alphanumerical 
disploy.  augmented  by  some  epeciol  options. 

Details  about  the  link  between  MICROSCOPE  and  <PL0T79>  ore  given 
in  section  S  and  your  system  manoger  may  have  put  out  suitoble  local 
information  in  this  respect.  For  the  remainder  of  this  subsection,  we  will 
assume  that  the  <PL0T79>  package  is  ovoiloble.  thot  you  hove  linked  it  with 
MICROSCOPE  and  your  own  program,  and  thot  you  hove  some  suitable  plotting 
device  available.  Appendix  tV  containes  copies  of  block  and  white  pen 
plots  corresponding  to  the  alphonumericol  plote  scattered  throughout  this 
monual . 

The  command  PLOT  puts  you  into  the  plot  submode.  The  precise 
kind  of  plot  you  will  eventuolly  generate  depends  upon  the  MICROSCOPE 
commands  given  so  far,  and  the  specifications  you  ore  going  to  moke  in  plot 
mode. 

We  first  discuss  the  PLOT  specifications.  After  entering  the 
PLOT  submode  you  will  see  a  list  of  current  options  on  the  screen.  You  con 
change  these  options  one  by  one.  After  each  change  on  updated  list  will  be 
displayed.  To  enter  an  option  change  enter  the  reference  number  of  the 
option  followed  by  the  appropriate  action  (os  prompted  by  the  program),  to 
activate  the  actual  plotting  process  enter  0.  ond  to  leave  the  plotting 
mode  and  return  to  command  mode  without  plotting  anything  enter  ony 
negative  number.  The  list  of  PLOT  options  is  left  unchanged  between  colls 
to  PLOT. 

There  are  12  PLOT  options: 

-1-  Title  (Default:  blank) 

A  title  printed  above  the  alphonumericol  display,  consisting  of 
up  to  72  characters  can  be  drawn. 

-2-  Legend  (Default:  blank) 

A  string  of  text  printed  below  the  alphanumerical  display, 
consisting  of  up  to  72  characters  can  olso  be  drown. 
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For  both  tha  titia  end  the  lagand,  o  dafoult  charactar  haight  ia 
providad  which  will  be  decraaaad  automatically  if  inaufficiant  apace  ia 
avai labia. 


-3-  Labeling  of  gropha  (Defouit:  ON) 

If  ON,  iabala  ara  drawn  giving  tha  degree  of  the  darivotiva  of 
each  graph.  On  a  crowded  diaplay.  thaaa  are  aomatimaa  awkwardly  plocad, 
and  marking  tha  grapha  (option  4)  may  be  preferable. 

-4-  Marking  of  Gropha  (Oafoult:  OFF) 

If  ON,  each  point  ot  which  a  darivotiva  hoa  bean  evaluated  ia 
marked  with  a  apacial  eymbol,  aimilorly  to  the  graphical  diaplay  on  tha 
tarmino  Tha  aymbola  for  tha  k-th  derivative  are 

k:  0  1  2  3  4  5  6 
Symbol;  dot  plue  oatariak  circle  x  aquora  triangle 

Marking  tha  grapha  inatead  of  labeling  them  hoa  the  advontaga 
that  it  ia  poaaibla  to  tall  from  any  amoll  port  of  tha  graph  to  which 
derivative  it  corraaponda,  making  it  unnacaaaory  to  trace  o  groph  from  tha 
point  of  intaraat  to  ita  label.  It  hoa  the  diaadvantoga  that  it  alowa  down 
aignif icontly  tha  apaad  with  which  tha  gropha  con  be  drown.  Alao, 
accentuating  a  graph  ia  only  poaaibla  with  tha  marking  on  (aaa  below).  Tha 
marking  aymbol  ia  repeated  in  tha  numerical  diaplay  if  tha  marking  option 
ia  on,  enabling  aaay  aaaociation  of  marka  ond  dagraaa  of  tongantlol 
darivotivaa. 


-5-  Drawing  tha  plot  in  color  (Default:  ON) 

Tha  graph  of  tha  tangential  darivotivaa  con  each  be  drawn  in  o 
diffarant  color,  if  aupportad  by  tha  exiating  hordwora.  In  the  moat 

axtrema  coaa,  when  all  darivotivaa  ora  to  be  plotted  in  different  colora, 
aaven  pane  for  their  grapha  plua  on  eighth  pen  for  tha  numaricol  diaplay, 
acalaa,  axea,  etc.  muat  be  euppliad. 

MICROSCOPE  aalacta  a  pen  by  calling  a  routine  SETPEN(I)  where 
I  B  1  corraaponda  to  tha  pan  drawing  tha  axea  etc.,  and  I  >  1  corraaponda 

to  tha  graph  of  tha  (I-2)th  tangential  derivative  (I  •  2,3 . B).  Thia 

range  of  I  haa  to  be  tronalotad  into  tha  oppropriota  <PL0T79>  ronga.  On 
aoffla  inatal lot iona  tha  pane  are  numbered  0  through  7,  on  othara  from  1  to 
8.  Tha  appropriate  SETPEN  routinaa  are  contained  in  tha  filaa  PEN07  ond 
PEN18,  reapactively  (aaa  aection  5).  If  only  a  amol ler  range  of  pena  ia 
ovailobla  you  may  want  to  write  your  own  pan  aalaction  routine,  tailoring 
it  parhopa  to  tha  apacific  ronga  of  darivotivaa  in  which  you  are 
intaraated.  Your  routine  muat  be  of  tha  following  form: 
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SUBROUTINE  SETPEN(I) 

INTEGER  I.J 

Cl  is  in  the  range  from  1  to  8 
J  •  sons  function  of  I 
C  J  is  in  the  ronge  of  pens  available 
C  on  the  local  plotting  device 
CALL  STPEN(J) 

RETURN 

END 

The  subroutine  STPEN  is  a  MICROSCOPE  routine  that  coiis  the 
<PLOT79  routine  SETPN(J)  if  the  volue  of  J  hos  changed  since  the  lost  coll 
to  STPEN.  SETPEN  Is  not  referenced  if  the  color  option  Is  turned  off. 

In  the  numerical  display,  the  ranges  of  the  tangential 
derivatives  ore  drawn  in  the  some  colors  os  their  graphs. 

Using  color  mokes  for  very  nice  drawings,  but  for  reproduction 
purposes  It  may  be  preferable  to  use  block  and  white  (like  the  <PL0T79> 
drawings  In  Appendix  IV). 

-8-  Drawing  a  page  number  (Default:  None) 

A  large  and  prominent  integer  con  be  drawn  below  the  legend.  We 
have  found  this  useful,  for  oxampla,  for  ossocioting  o  drawing  with  the 
rsforsnes  number  of  o  point  In  o  list  (s.g.  the  point  of  examination),  and 
associating  or  defining  o  figure  number  (like  in  this  manual). 

No  page  number  Is  drown  If  its  value  is  zero. 

-7-  Left  bottom  lobe  I  (Defoul t:  none) 

A  lobel  of  up  to  20  characters  of  smoll  type  can  be  drawn  in  the 
left  bottom  corner  of  the  plot. 

-8-  Right  bottom  label  (Default:  none) 

A  similar  label  con  be  drown  in  the  right  bottom  corner  of  the 

plot . 


-9-  Time  (Default:  off) 

The  time  at  which  the  plot  woe  generated  can  be  drown  in  the  top 
left  corner  of  the  plot. 


-10-  Dote  (Oefoult:  OFF) 

Similarly,  the  date  of  the  plot  generotion  can  be  drawn  in  the 
top  right  corner. 
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Both  time  and  dote  ore  printed  in  the  some  size  os  the  bottom 
labels.  They  ore  sometimes  useful  for  identification  and  record  keeping. 

-11-  Numericol  Display  (Oefoult;  ON) 

If  desired,  the  numericol  disploy  con  be  turned  off.  This  is 
sometimes  useful  to  obtain  quick  drawings. 

-12-  Prome  (Defoul t:  ON) 

If  desired,  the  frame  that  is  normolly  drown  oround  the  entire 
<PL0T79>  drawing  con  be  turned  off. 

We  now  turn  to  the  MICROSCOPE  options.  Essentially  everything 
that  oppeors  In  the  o I phonumer i co I  screen  display  is  also  contoined  in  the 
<PL0T79>  drawing,  although  it  moy  ossume  a  slightly  different  form. 

-1-  Horizontol  scale 

If  0  horizontal  scale  is  activated  (by  the  OSCALE  commond)  then  a 
standard  labeled  horizontal  scale  is  drown  on  the  bottom  of  the  grophical 
<PL0T79>  drawing.  The  labels  on  the  scale  give  the  distonce  from  the  point 
of  examination  in  units  of  the  direction  of  investigation.  Notice  thot 
this  may  be  the  normalized  version  of  the  direction  given  in  the  numerical 
display,  depending  on  whether  the  normalization  option  is  on.  If  it  is, 
then  the  scale  in  the  <PL0T79>  drawing  gives  the  Euclidean  distonce  from 
the  point  of  examination. 

-2-  X  Axis 

A  horizontal  line  is  drown  through  the  center  of  the  graphical 
display  if  the  x  axis  is  turned  on  (using  DXAXIS).  This  is  sometimes 
useful  for  illustrating  symmetry  properties. 

-V-  Marking  the  center  of  the  disploy 

If  this  option  is  turned  on  (using  DCENTER)  the  center  of  the 
graphical  display  is  marked  with  o  cross. 

-4-  Accentuoting  a  groph 

If  a  graph  is  to  be  occentuoted  its  marks  (os  discussed  above) 
are  drawn  twice  their  normal  size.  However,  if  color  is  available,  a 
better  emphasis  con  usually  be  obtoined  by  using  a  special  color. 
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4.1S  Ueer  Intervention 


USER  enter  ueer  node  -  -  - 

Upon  receiving  the  connond  USER.  MICROSCOPE  calls  a  subroutine 
SUBUSR  which  is  often  a  dummy  routine,  but  which  can  be  written  to  suit 
special  purposes.  SUBUSR  has  no  parometers,  but  it  can  communicate  with 
MICROSCOPE  routines  or  your  own  routines  via  COtMON  blocks.  In  utilizing 
this  facility  it  is  importont  to  keep  in  mind  that  ordinarily  MICROSCOPE 
has  no  means  of  determining  that  the  trial  function  hos  changed  and  that 
previaus  evaluations  of  it  ore  no  longer  valid.  Therefore  the  command 
FORCE  may  have  to  be  used  rother  thon  CO. 

Five  examples  of  possible  opplicotions  follow: 

-1-  Mixed  portiol  dsrivotives 

If  derivatives  other  thon  pure  tangential  derivatives  of  purs 

cross  derivatives  are  to  be  examined  this  con  bs  occomplished  by  dsflning 
the  trial  function  to  be  the  derivative  thot  is  to  bs  displayed,  or  to  be 
further  differentiated,  ond  moke  it  dependent  upon  parameters  that  ors 
modified  by  SUBUSR  as  needed.  In  this  context,  it  may  be  convenient  to  use 
the  routine 

DOUBLE  PRECISION  FUNCTION  OERIV(F.NDIM,P. IOROER,DIR,H) 

DOUBLE  PRECISION  F.P(1 ) .DIR(1 ) ,H 
INTEGER  NDIM.IORDCR 

RETURN 

END 

which  computes  ot  the  point  P  the  lORDERth  derivative  in  the  direction  DIR 
of  the  double  precision  function  F,  using  the  formulas  given  in  section  1.2 
with  the  discretization  parameter  H.  F  must  be  declared  EXTERNAL  in  the 
calling  program,  and  is  of  the  same  form  os  the  trial  function  described  in 
section  3.1.  NDIM  is  the  number  of  independent  voriobles  F  depends  on  ond 
must  not  exceed  3.  DERIV  is  used  within  MICROSCOPE  (to  compute  any  cross 
derivatives)  and  con  therefore  be  readily  called  by  your  own  programs. 

Obviously,  the  modified  trial  function  may  itself  be  a  tangential 
derivotive.  Then  the  range  of  derivatives  that  con  be  displayed  is  shiftsd 
from  the  normal  9th  through  6th  order  to  a  higher  one. 

-2-  Switching  Triol  Functions 

If  several  trial  functions  ore  to  be  investigated  in  a  single 
session  (or  with  a  single  linking  and  looding  process)  they  con  all  be 
written  into  one  DOUBLE  PRECISION  FUNCTION,  and  distinguished  by  o 
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paromster  that  is  passed  in  a  conenon  block,  and  that  may  be  altered  by 
SUBUSR.  The  examples  in  this  manual  or*  hondled  in  this  manner  (see 
appendix  III  for  details). 

-3-  Switching  Data  or  Parometers 

Often  the  trial  function  depend*  upon  dota  or  parameter*  that  con 
be  changed  by  SUBUSR  without  having  to  leave  MCRSCP. 

-4-  Fi I*  Handl ing 

File  handling  is  system  dependent.  You  moy  find  it  convenient  to 
handle  the  assignments  of  files  or  devices  to  chonnel  numbers  vio  USER. 


4.16  Pflffluli  apttingg 

In  this  subsection,  E  denotes  a  vector  of  1  to  3  component*  that 
are  all  equal  to  1.  The  default  settings  for  MICROSCOPE  are  a*  follows; 

Number  of  independent  voriobles:  2 

Point  of  examination:  8  (the  origin) 

Direction  of  Investigation;  £ 

Discretization  porometer;  ID-4 

Window  width:  12 

Cross  di rect ion:  E 

Order  of  Cross  Derivative;  0  (no  cross  derivative) 

Order  of  Tangential  Derivatives  to  be  displayed:  0 

No  scale  or  x  axis  is  drawn,  the  center  of  the  display  is  unmorked,  and 
no  derivative  is  accentuated. 

PLOT  options:  see  section  4.14 
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5.  Inttatlalion  Guidm 


5.1  Paierlatlon  of  tha  Pnehaoa 

In  this  ••etion,  which  la  diractad  mainly  to  tha  local  ayatam 
manogar,  wa  daacriba  tha  filaa  auppllad  with  tha  MICROSCOPE  packaga,  and 
how  thoy  con  ba  put  togathar.  Uauolly,  it  will  ba  tha  baat  proeadura  to 
aalact  or  writa  filaa  that  ganarota  tha  meat  powerful  ond  convaniant  local 
varalon  of  MICROSCOPE,  to  put  thorn  into  one  file  or  a  library  ond  to  moke 
tha  oppropriota  Information  available  to  your  local  uaara. 

MICROSCOPE  ia  0  portable  aat  of  FORTRAN  routinaa.  The  uaar  muat 
aupply  0  moln  program  and  one  or  two  aubroutinaa  that  ora  apacific  to  the 
application.  Thua  uaara  muat  hove  occaaa  to  at  laoat  tha  ralocotobla 
varaiona  of  a  aultobla  aubaat  of  MICROSCOPE,  if  not  tha  aourca  codaa. 

You  ahould  hove  obtolned  a  mognatic  tope  containing  21  filaa  In 
tha  following  aaquanca: 


1 .  HELPLC 

(902) 

2.  HELPUC 

(902) 

3.  LC 

(1889) 

4.  MAIN 

(48) 

5.  MANLC 

(5150) 

6.  MANUAL 

(7456) 

7.  MANUC 

(5150) 

6.  PENe7 

(9) 

9.  PENIS 

(8) 

10.  PENOMY 

(3) 

11.  PLTOMY 

(169) 

12.  PLTLC 

(978) 

13.  PLTUC 

(980) 

14.  POS 

(14) 

15.  SCREEN 

(89) 

16.  SCROLL 

(86) 

17.  SUPPRT 

(3894) 

16.  UC 

(1862) 

19.  USROMY 

(8) 

20.  USRLC 

(271) 

21.  USRUC 

(270) 

Numbara  in  paranthaaaa  give  tha  number  of  linaa  (of  up  to  60 
choroctara  aoch)  in  tha  filaa. 

Below,  tha  filaa  ora  divided  into  aavarol  groupa.  One  or  two 
filaa  hove  to  ba  picked  from  aoch  of  tha  groupa  1  through  6.  Within  a 
group,  filaa  ora  liatad  by  dacreoaing  convanianca,  and  incraoaing 
portability.  If  tha  bottom  filaa  ora  choaan  from  each  group  a  complataly 
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portable  version  of  MICROSCOPE  (except  for  the  use  of  a  few  non— FORTRAN 
characters,  namely  ?.<,>, I,:.:)  is  obtoined. 

-1-  Case  Dependent  Support  Routines 

LC:  Contains  those  subroutines  that  require  the  support  of 

lower  case  letters  in  order  to  work  properly.  This  feature  is  not 
portable,  but  widely  available,  and  greatly  increoses  the  convenience  in 
using  MICROSCOPE- 

UC:  To  be  used  in  case  of  LC.  if  lower  cose  letters  ars  not 

avoiloble.  UC  is  supplied  os  port  of  the  pockoge,  but  it  can  also  be 
generated  from  LC  by  replacing  all  lower  cose  letters  with  their  upper  cose 
equivalents,  and  by  replacing  the  subroutine  LCUC  by  o  dummy  routine. 
(Upper  cose  versions  of  all  other  files  con  be  generated  from  their  lower 
case  counter  ports  by  replacing  lower  cose  letters  by  uppercase  ones 
everywhere.  The  leading  comments  describing  the  file  will  of  course  no 
longer  apply.) 

If  UC  is  used  in  cose  of  LC  oil  output  will  be  in  upper  case 
letters,  and  lower  cose  input  will  not  be  recognized. 

LC  and  UC  contain  their  respective  versions  of  the  master  routine 
MCRSCP,  which  is  the  first  routines  in  eoch  of  the  two  fiies.  Ali  other 
routines  then  follow  alphabetically. 

-2-  Cose  Independent  Routines 

SUPPRT:  This  file  contains  the  buik  of  the  MICROSCOPE  pockoge. 

consisting  of  completely  portobte  routines. 

-3-  The  File  Access  Routine 

POS;  This  file  contains  o 
SUBROUTINE  POS(IDEV,N) 

which  positions  the  device  lOEV  at  the  N-th  record,  such  that  subsequent 
reading  from  IDEV  yields  the  N-th  record.  This  is  accomplished  by 
rewinding  IDEV  and  then  reading  a  dummy  chorocter  from  the  first  N-1  lines 
(records)  in  IDEV.  If  this  at  oil  time  consuming  it  moy  be  worthwhile  to 
replace  POS  by  a  faster  system  dependent  routine. 

-4-  The  Screen  Display  Pockoge 

SCREEN:  To  be  used  if  screen  editing  is  locaily  avoiloble,  ond 

highly  recommended  over  its  alternative,  SCROLL,  described  below.  To  toke 
odvontoge  of  screen  editing,  the  following  four  smoll  interface  routines 


Instol lotion  Guide 


5.1  Description  of  the  Package 


page  76 


MICROSCOPE: 
U«er‘3  Manual 


A  Software  Syetem  for  Multivariate  Analyeie  poge  77 

Vereion  of  April  1984 


must  be  written  for  each  possible  conputer/termi nal  conbination  (in  oil 
cases,  IDEVCE  is  the  channel  number  of  the  relevont  terminal): 

SUBROUTINE  CLSCRN( IDEVCE, I .J) 

INTEGER  IDEVCE, I. J 

C  Clear  the  screen  of  the  device  with  number  IDEVCE  below  the  point  with 
C  the  coordinates  (I,J}  where  I  is  the  horizontol  coordinote  meoeured 
C  from  the  left  of  the  screen  starting  ot  1,  and  J  is  the  verticoi 
C  coordinate  measured  from  the  top  of  the  screen  storting  at  1. 

RETURN 

END 

SUBROUTINE  PCURSR( IDEVCE, I ,J) 

INTEGER  IDEVCE, I, J 

C  Position  the  cursor  at  the  position  (I,J)  of  the  device  with  number 
C  IDEVCE,  where  l.J  are  as  for  CLRSOM 

RETURN 

END 

SUBROUTINE  PLCHRS( IDEVCE, I ,J ,N, CHARS) 

INTEGER  IDEVCE, I, J,N,CHARS(13S) 

C  Write  the  string  of  N  choroetere  contained  in  A1  format  in  the  integer 
C  array  CHARS,  starting  at  the  ioeation  (I.«l).  l.J,  and  IDEVCE  ore  os 
C  in  CLRSCRN. 

RETURN 

END 

SUBROUTINE  eLSCRN( IDEVCE) 

C  Blank  the  screen  of  the  device  with  number  IDEVCE 

RETURN 

END 

SCROLL:  SCREEN  con  be  reploced  with  SCROLL  if  the  SCREEN  package 

cannot  be  used.  This  will  have  the  consequence  that  the  entire  screen 
display  is  regenerated  after  each  change,  rather  than  just  modified. 

-5-  The  PLOT  pockage: 

PLTLC:  This  file  contains  the  interfoce  between  MICROSCOPE  and 

<PLOT79>.  The  routines  use  and  recognize  lower  case  letters. 

PLTUC:  The  equivalent  of  PLTLC,  but  without  the  lower  cose 

letters.  Note,  however,  that  <PLOT79>  supports  Hershey  Fonts  that  include 
lower  cose  letters  regardless  of  the  computing  system.  Hence  the  <PLOT7 
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drowings  generated  by  PLTUC  do  contain  lower  cose  letters.  The 
conmunicot ion  with  the  tsrninol  vio  the  coMtond  PLOT  is  completely  in  terms 
of  upper  cose  l^^tters,  however. 

With  both  the  PLTLC  or  PLTUC  package  o  pen  selection  routine  also 
hoe  to  be  supplied  (see  section  4.14).  The  MICROSCOPE  pockoge  contoins 
three  routines  in  separate  files:  PEN18  corresponding  to  8  pens  labeled  1 
to  8,  PEN07  corresponding  to  8  pens  lobeled  0  to  7,  and  PENOMY  to  bs  used 
if  color  is  not  available  (e.g.  on  dot  matrix  printer).  However,  you  may 
wish  to  leave  it  up  to  the  user  to  select  or  write  his  own  pen  selection 
rout ines. 


PLTOMY:  To  be  used  if  <PLOT79>  is  not  available.  In  this  case 

only  olphonumericol  screen  displays  liks  thoss  shown  throughout  this  manual 
can  be  obtained.  However,  the  dummy  routine  simulotes  the  interaction  with 
<PLOT79>  up  to  the  point  at  which  the  actual  drowing  takes  plocs.  This 
feature  is  useful  for  the  development  of  MICROSCOPE  programs,  and  mokss  It 
possible  to  run  programs  without  <PLOT79>  even  if  they  coll  the  pockoge. 

There  is  only  the  upper  cose  version  of  PLTOMY. 

-6-  The  Test  packags 

Ordinorily,  users  will  not  be  interested  in  the  specific  routine 
SUBUSR  employed  in  the  generation  of  this  monuol.  They  will  wont  to  write 
their  own.  However,  the  manual  version  of  SUBUSR.  ond  the  routines  called 
by  it,  ore  contained  in  the  files: 

USRLC:  contains  the  lower  cose  version. 

USRUC:  contains  the  upper  cose  version,  to  be  used  if  lower 

cose  letters  are  not  available. 

USRDMY:  contains  a  dummy  SUBROUTINE  SUBUSR,  which  disables  the 

USER  conmond.  If  the  MICROSCOPE  routines  ore  held  in  a  library  it  moy  be 
possible  to  keep  this  dummy  routine  in  it  os  a  stondby  if  the  user  dees  not 
supply  his  own  version. 


-7-  The  <PLOT79>  package 

Obviously  the  <PL0T79>  package  must  bs  linked  with  MICROSCOPE  if 
it  is  to  be  utilized.  This  is  a  lorge  graphics  pockoge  that  con  be  eosily 
interfaced  with  almost  any  grophics  device.  It  con  be  obtained  for  o 
nominal  charge  from  Or.  Nelson  H.F.  Beebe,  Department  of  Physics, 
University  of  Utah,  Salt  Lake  City.  Utah  84112,  Tel.:  801-S81-5254. 

-8-  Doto  files 
Sample  Driver  Program 
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Yh*  fil*  MAIN  contoin*  th«  drivar  program  aa  it  wo*  utad  in  th* 
gonaration  of  thia  manual.  It  i»  liatod  oa  o  doto  fila  bacouaa  tha  program 
naada  to  ba  modifiad  to  auit  locol  ci rcumatoncaa.  In  particular,  tha  OPEN 
atotamant  ia  not  portobla. 


InUtacLbn  iltLa-fiJM 


Tha  poekaga  Ineludaa  two  halp  fllaa  that  contain  ralavont  halp 
information  (but  no  coda).  Again,  thara  ora  tha  lowar  and  uppar  coaa 
varaiona: 


HELPLC 


HCLPL  : 


Ona  of  thaaa  fiiaa  (prafarably  HELPLC)  muat  ba  raodobla  by  tha 
uaar  programa,  ond  tha  uaara  muat  ba  informad  how  to  occaaa  tham.  It  la 
poaaibla  to  daaign  a  aaguanca  of  MICROSCOPE  commonda  that  ovarwritaa  tha 
halp  flla,  and.  If  poaaibla,  tha  uaar  ahould  ba  daniad  writa  occaaa  to  it. 
You  may  olao  wont  to  kaap  aoaily  accaaalbla  bockup  coplaa  of  thaaa  fllaa. 

Tha  halp  fiiaa  aoch  contain  tha  following  Itama: 

1.  On  tha  firat  lina  thara  ora  thraa  intagara  I,M,N  In  313 
format.  Thaaa  hava  tha  fol lowing  maoninga:  I  ia  tha  numbar  of  llnaa 
oceuplad  by  tha  tobla  of  commando  that  con  ba  occaaaad  by  tha  commond 
HSUMRY.  M  ia  tha  total  numbar  of  linaa  occupiad  by  tha  data! lad  information 
occaaaad  by  HELP.  N  ia  tha  numbar  of  oval  labia  commando. 


2.  Tha  naxt  four  linaa  ora  occupiad  by  N  numbara  writtan  in  2812 
format.  Tha  k-th  ouch  numbar  givaa  tha  langth  (In  linaa)  of  tha  datollad 
HELP  information  on  tha  k-th  commond.  Thaaa  N  numbara  muat  odd  up  to  M. 


3.  Tha  naxt  I  linaa  contain  tha  HSUMRY  Information 

4.  Tha  naxt  M  linaa  contoln  tha  dotal  lad  HELP  information. 


S. 

contoina  tha 


Tha  naxt  lina  contoina  on  intagar 
numbar  of  linaa  occupiad  by  tha  nawa. 


J  in  13  format  that 


6.  Tha  loot  J  linaa  contain  tha  nawa. 


Monugl  fllti 

MANUAL;  contoinaa  a  varbotim  machine  raodobla  copy  of  thia 
manual.  It  ia  writtan  with  FORTRAN  corrioga  control  ond  raody  for  printing 
(ot  60  linaa  par  page  and  up  to  75  choroctara  per  lina).  Underlining  ia 
oecompllahad  by  having  tha  undarlinaa  in  o  aaporota  lina  with  a  “4" 
character  in  tha  firat  column,  indicating  that  tha  lina  ia  to  ovarwrita  tha 
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prectding  one.  There  ie  no  upper  ease  vereion  of  thie  file.  If  lower  caee 
lettere  are  not  availoble  a  xeroxed  copy  of  the  manual  rather  than  a 
printed  vereion  should  be  distributed. 

This  installation  guide  is  chapter  five  of  the  monuol. 

MANLC:  Thie  file  contains  o  version  of  the  manuol  that  is 
suitable  for  on-line  reading  and  editing.  Poge  breaks  and  underlines  have 
been  eliminated,  otherwise  the  text  is  identical  to  thot  in  MANUAL. 

MANUC:  The  upper  case  equivalent  of  MANLC. 


S.7  A  Key  to  Selecting  the  Appropriote  Files 

The  following  is  a  quick  key  to  selecting  the  oppropriote  subset 
of  the  above  files.  Proceed  by  answering  the  queetions  ond  following 
instructions  until  you  reach  a  line  giving  the  relevant  file  nomes  or 
descriptions.  The  pen  selection  routine  is  described  in  section  4.14  (and 
can  be  supplied  by  the  user  if  necessary),  oil  other  ingredients  ore 
described  in  this  Installation  guide. 


1 . 

Are  lower  cose  letters  available  ^ 

on  your 

mach i ne? 

Yes 

GO 

TO 

2 

No 

GO 

TO 

3 

2. 

Have  you  supplied  an  Interface 

to 

ocroon 

edi t ing? 

Yes 

GO 

TO 

4 

No 

GO 

TO 

5 

3. 

Have  you  supplied  an  interface 

to 

tcroon 

editing? 

Yes 

GO 

TO 

6 

No 

GO 

TO 

7 

4. 

Is  <PL0T79>  available? 

Yes 

GO 

TO 

8 

No 

GO 

TO 

9 

5. 

Is  <PL0T79>  available? 

Yes 

GO 

TO 

18 

No 

GO 

TO 

11 

6. 

Is  <PLOT79>  available? 

Yes 

GO 

TO 

12 

No 

GO 

TO 

13 

7. 

Is  <PLOT79>  ovoilable? 

Yes 

GO 

TO 

14 

. 

No 

GO 

TO 

15 

8.  LC,  SUPPRT,  POS  or  local  equi volant,  SCREEN,  PLTLC,  pen  selection 
routine,  possibly  USRLC,  screen  interface,  HELPLC,  MANUAL,  MANLC 

9.  LC,  SUPPRT,  POS  or  local  equivalent,  SCREEN,  PLTOMY,  possibly  USRLC, 
screen  interfoce,  HELPLC.  MANUAL.  MANLC 
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10.  LC,  SUPPRT,  POS  or  locol  aquivolant,  SCROLL,  PLTLC,  pen  ••laction 

routina,  potsibly  USRLC,  HELPLC.  MANUAL,  MANLC 

11.  LC,  SUPPRT,  POS  or  local  oquivolant.  SCROLL,  PLTDMY,  poatibly  USRLC, 
HELPLC,  MANUAL,  MANLC 

12.  UC,  SUPPRT,  POS  or  local  aquivolant,  SCREEN,  PLTUC,  pan  aalaction 

routina,  possibly  USRUC,  scrssn  intsrfoca,  HELPUC,  MANUC 

13.  UC,  SUPPRT,  POS  or  loco!  aquivolant,  SCREEN,  PLTDMY.  possibly  USRUC, 
scrssn  intsrfoca.  HELPUC,  MANUC 

14.  UC.  SUPPRT,  POS  or  local  aquivolant,  SCROLL,  PLTUC.  pan  salaction 

routina.  possibly  USRUC.  HELPUC,  MANUC 

15.  UC.  SUPPRT.  POS  or  local  aquivolant.  SCROLL,  PLTDMY,  possibly  USRUC, 
HELPUC.  MANUC 
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Aoeendix  1:  Numerical  Pi f ferent i at i on 

In  this  appendix,  we  will  discuss  in  some  detoil  the  properties 
of  the  formulas  listed  in  section  1.2.  Without  loss  of  generality  we 
assume  the  point  of  exominotion  is  the  origin.  We  also  restrict  our 
attention  to  purely  tangential  derivatives.  Thus  we  may  assume  there  Is 
only  one  independent  vorioble.  Of  porticulor  interest  ore  step 
discontinuities  in  certain  derivatives  (for  other  types  of  reduced 
smoothness  see  the  examples  at  the  end  of  Appendix  III).  Thus  we  consider 
the  trial  function 


F(x)  “  PHI(x)  +  S(m,eto.x) 

where  PHI(x)  is  o  smooth  function  (i.e.  all  its  relevant  derivotives  are 
continuous)  and  S  is  determined  by  the  requirements  thot  its  m-th 
derivative  be  piecewise  constant  and  have  a  jump  discontinuity  of  magnitude 
eta  at  x  >  8.  Explicitly,  S  is  defined  by 


m 

S(m,eta,x)  >0  if  x  <  0  and  eta*x  /ml  otherwise 

where  the  exclamation  mark  denotes  the  foctoriol. 

To  simplify  matter*  we  wi 1 1  assume  S  is  evoluotsd  sxoctly,  ond 
all  round-off  errors  occur  in  PHI.  The  properties  of  the  numerical 
differentiation  formulas  depend  very  critically  on  the  relative  sizes  of 
the  derivatives  of  PHI,  about  which  little  con  be  ossumed  in  generol.  Our 
approach  will  be  to  derive  general  farmulas  incorporating  the  derivatives, 
but  to  illustrate  the  relevant  phenomena  by  using 

X 

PHI(x)  -  e 

since  the  derivatives  of  the  exponentiol  ore  all  identical  (assuming  the 
normalization  option  is  on). 

We  wilt  denote  the  degree  of  o  deri votive  by  porenthesized 
superscripts. 

A  simple  calculation  shows  that,  for  k  m,  the  approximation 

(k)  (k) 

0  of  f  at  X  •  0  is  given  by  the  fundamental  equation 
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(k)  (k) 

0  -  PHI  trua  valua 


2  (k+2) 

h  *gaiima  •PHI 
k 

-k 

4  tau*h  •(1-t-aba(PHI))*atpha 

k 


truncation  arror 

(1) 

round-off  tarn 


la— k 

•f  ato'dalto  'h  dalto  tarn 

m,k 

If  k  <  III  tha  dalto  tarn  la  obaant .  Wa  dlacuaa  tha  individual 
tarna  in  turn: 


(k) 

PHI 


-  Tha  trua  valua 
Thia  la  tha  axpraaalon 
la  avaluatad  at  zaro. 


wh I ch  wa 


ora 


trying  to  epproxinota. 


-  Tha  truncotion  arror. 


Thia  tarn  can  ba  obtainad  by  axpanding  tha  arror  into  o  Taylor 
aariaa  about  h  ■  0,  oa  la  dona  in  any  taxtbook  on  nunarieol  anolyaia.  Tha 
(k+2) 

function  PHI  la  to  ba  thought  of  oa  having  baan  avaluatad  at  a 

auitobla  (unknown)  point  in  tha  donoin  of  tha  diaploy.  Tha  gonno  ora 

k 


conatonta  givan  in  tabla  1. 


-  Tha  round-off  tarn 


Wa  aaauna  rounding  ia  taking  ploca  aa  daacribad  in  aaction  0, 
i.a.  if  X  ia  tha  quontity  to  ba  roundad  to  D  digita,  aoy,  than  x  ia 
raplocad  by 


—0  — 0 

z  :■  (1+apa  *10  )»x  +  apa  alO 

1  2 


(2) 


whara  apa  and  apa  ora  random  numbara  batwaan  -1  ond  +1 . 
1  2 


in  (1)  ia  o  auitobla  ovaroga  of  tha  apa  and  apa  tarna  in 

1  2 


Tha  quantity  tou 
(2).  Wa  will 
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sat  it  to  19  The  function  PHI  ia  ogoin  evaluated  at  a  auitoble  point  in 
the  diaploy,  ond  the  alpha  are  the  aufli  of  the  abaolute  valuea  of  the 

k 

coefficienta  in  the  k-th  di f ferent lot  ion  formula.  Their  valuea  are  liated 
in  table  1 . 

-  The  del ta  term. 

We  uae  the  expreaaion  "delta  term"  becauae  it  ariaea  from 
differentiating  the  atep  function  a(m.eta,x),  and  the  firat  derivative  of  a 
atep  function  ia  o  Dirac  delta  function.  The  delto  ore  conatanta  that 

m,k 

can  be  read  off  the  MICROSCOPE  plota  in  aection  4.2.  They  are  olao  given 


in  toble  1 . 

k  alpha 

k 

gamma 

k 

m  •  9 

1 

delta 

2 

m.k 

3 

4 

5 

6 

1 

2 

1/6 

1/2 

1 

- 

- 

- 

- 

- 

2 

4 

1/12 

1 

1 

1 

- 

- 

- 

- 

3 

24 

1/16 

4 

2 

1 

1 

- 

- 

- 

4 

256 

1/24 

48 

16 

3 

1 

1 

- 

- 

5 

2439 

1/27 

365 

61 

14 

3 

1 

1 

- 

6 

46656 

1/36 

7299 

1458 

146 

27 

4 

2 

1 

Table  1:  Conatanta  for  Differentiation  Formulae 


We  now  oddreaa  three  queationa: 

1.  When  do  round-off  effecta  become  viaible? 

2.  When  do  round-off  effecta  exceed  the  truncotion  error? 

3.  How  email  a  diacont inuity  can  be  detected  with  MICROSCOPE? 
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t.1  Th>  Round-Off  Thr««hold 

W*  are  concerned  with  the  queetion  of  when  round-off  effecte 
begin  to  obecure  the  graph  of  o  derivotive.  More  precieely,  we  oek  when 
the  round-off  tern  beconea  lorger  then  o  certain  froction  of  the  rite  in 

M 

the  graph  of  the  derivative  PHI  over  on  interval  of  length  e.  An 
equation  for  the  critical  atep  cize.  which  conetitutee  a  limit  on  h  beiow 
which  no  cieor  gropha  can  be  obtained,  ie  generated  by  ignoring  the  delta 

(k) 

term  and  the  truncation  error,  and  approximating  the  riae  of  PHI  by  the 
mean  value  theorem.  Thia  laoda  to: 

-k  (k+1) 

pal  eh  •(1-faba(PHI))*tou*alpho  >  2*h*abt(PHI  )/w 
1  k 

(k+1) 

where  w  *  2h/a  ie  the  window  width,  PHI  ogoin  ie  evoiuoted  at  a 

euitoble  point,  and  pai  ie  a  eofety  factor.  The  iorgar  pai  ,  the  clearer 
1  1 
(W) 

will  be  the  graph  of  PHI 

Solving  for  h  yielda: 


(k+1)  ( 

h  ■  [pal  (1+aba(PHI))»tau*olpha  •m  /(2*aba(PHI  ))]  1/(k+ 
1  k 


Table  2  containa 
number  of  digite  carried)  and 

exponential.  All  derivotivea 


threahold 
pai  -  5. 
1 

hove  been 


voluea  for  vorloua  voluea  of  D  (the 
for  the  coae  that  PHI  la  the 

evoiuoted  at  x  ■  B. 
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Digits  (D) 

Degree  of  Derivot i 

ive  (k) 

12  3  4 

5 

6 

1 

3.46D^e 

2.8804^8 

3.460488 

4.340488 

4.940488 

6.880+88 

2 

1 .180408 

1  .3404-88 

1 . 950488 

2.740489 

3.370+89 

4.320488 

3 

3.46D-81 

6.210-81 

1.180488 

1 . 730488 

2.290+88 

3.110+88 

4 

1.180-81 

2.880-81 

6.160-81 

1 .890488 

1 . 560+88 

2.240+88 

5 

3.460-82 

1 .340-81 

3.460-81 

6.880-81 

1  .860488 

1 .610+88 

6 

1.180-82 

6.210-82 

1 .950-81 

4.340-81 

7.250-81 

1.160+88 

7 

3.460-83 

2.880-82 

1.180-81 

2.740-81 

4.940-91 

8.340-81 

8 

1.180-83 

1 . 340-82 

6.160-82 

1.730-91 

3.370-81 

6.880-81 

9 

3.460-84 

6.210-83 

3.460-82 

1 .890-81 

2.290-81 

4.320-81 

18 

1.180-84 

2.680-83 

1 . 950-82 

6.880-82 

1 .560-81 

3.110-81 

11 

3.460-85 

1 . 340-83 

1.180-82 

4.340-82 

1  .860-81 

2.240-81 

12 

1 . 180-85 

6.210-84 

6.160-83 

2.740-82 

7.250-82 

1.610-81 

13 

3.460-86 

2.880-84 

3 . 460-83 

1 .730-82 

4.940-82 

1.160-81 

14 

1.180-86 

1 . 340-84 

1.950-83 

1.890-82 

3.370-82 

8.340-82 

15 

3.460-87 

6.210-85 

1.180-83 

6.880-83 

2.290-82 

6.880-82 

16 

1.180-87 

2.880-85 

6.160-84 

4.340-83 

1 . 560-82 

4.320-82 

17 

3.460-88 

1 .340-85 

3.460-84 

2.740-83 

1 .860-82 

3.110-82 

18 

1.180-88 

6.210-86 

1 . 950-84 

1.730-83 

7.250-83 

2.240-82 

19 

3.460-89 

2.68D-86 

1.180-84 

1 .890-83 

4.940-83 

1.610-82 

28 

1.180-89 

1.340-86 

6. 160-85 

6.880-84 

3.370-83 

1.160-82 

21 

3.460-18 

6.210-87 

3.460-85 

4.340-84 

2 . 290-83 

8.340-83 

22 

1.180-18 

2.880-87 

1.950-85 

2.740-84 

1 . 560-83 

6.880-83 

23 

3.460-11 

1 .340-87 

1.180-85 

1 .730-84 

1.860-83 

4.320-83 

24 

1.180-11 

6.210-88 

6.160-86 

1 .890-84 

7 . 250-84 

3.110-83 

25 

3.460-12 

2.880-88 

3.460-86 

6.880-85 

4.940-84 

2.240-83 

26 

1.180-12 

1 .340-88 

1.950-86 

4.340-85 

3.370-84 

1 .610-83 

27 

3.460-13 

6. 210-89 

1.190-86 

2.740-85 

2 . 290-84 

1.160-83 

28 

1.180-13 

2.860-89 

6.160-87 

1 .730-85 

1 .560-84 

8.340-84 

29 

3.460-14 

1 . 340-89 

3.460-87 

1.890-85 

1 .860-84 

6.880-84 

38 

1.180-14 

6.210-18 

1 .950-87 

6.880-86 

7.250-85 

4.320-84 

31 

3.460-15 

2.680-18 

1.180-87 

4.340-86 

4.940-85 

3.110-84 

32 

1.180-15 

1 .340-18 

6.160-88 

2.740-86 

3.370-95 

2.240-84 

33 

3.460-16 

6.210-11 

3.460-88 

1 .730-86 

2 . 290-85 

1  .610-84 

34 

1.180-16 

2 . 860-1 1 

1 .950-88 

1 .890-86 

1 .560-85 

1.160-84 

35 

3.460-17 

1 .340-11 

1.190-88 

6.880-87 

1 . 860-85 

8.340-85 

36 

1.180-17 

6.210-12 

6.160-89 

4.340-87 

7 , 250-96 

6.880-85 

37 

3.460-18 

2.880-12 

3.460-89 

2.740-97 

4 . 940-86 

4.320-85 

38 

1.180-18 

1  .340-12 

1 .950-89 

1 .730-87 

3 . 370-86 

3.110-85 

39 

3.460-19 

6.210-13 

1.190-89 

1 .890-87 

2 . 290-86 

2.240-85 

48 

1.180-19 

2.880-13 

6.160-18 

6 . 880-88 

1 . 560-86 

1 .610-85 

Table  2:  Guide  to  selecting  a  value  of  the  discretization  poronete 
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1.2  Round-off  vraua  Truneotlon  Errori 

In  working  with  MICROSCOPE,  ono  will  raroly  bo  intoraatad  in 
obtaining  aecurata  numarlcal  valuaa  of  darlvotlvaa.  Rothor,  quolltatlva 
faaturaa  will  ba  of  primary  concarn.  Howavar,  if  accuracy  la  paramount, 
than  tha  diacratiaation  paromatar  ahould  ba  picked  oa  amoll  aa  poaaibla 
without  tha  round-off  tarm  awomping  tha  truncation  arror.  Tha  optimal  atop 
aiza  ia  obtained  by  equating  tha  round-off  tarm  with  tha  truncation  arror 
and  aalving  for  h.  Thia  yialda: 


(ka2)  1/(k-t-2) 

h  ■  C(iaoba(PHI))*tau*aipha  /  (gamma  •PHI  )] 

k  k 


Table  3  I  late  valuaa  of  h  for  tha  exponent iai. 
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Digits  (D) 

1 


1 

1 . 690488 

2 

7.830-81 

3 

3.630-81 

4 

1.690-81 

5 

7.830-82 

6 

3.630-82 

7 

1 .690-82 

8 

7.830-83 

9 

3.630-83 

18 

1 . 690-83 

11 

7.830-64 

12 

3.630-64 

13 

1 . 690-84 

14 

7.830-85 

15 

3.630-85 

16 

1 . 690-85 

17 

7.830-86 

18 

3.630-86 

19 

1 . 690-86 

28 

7,830-87 

21 

3.630-87 

22 

1.690-87 

23 

7.830-88 

24 

3.630-88 

25 

1 .690-88 

26 

7.830-89 

27 

3.630-89 

28 

1 .690-89 

29 

7.830-16 

36 

3.630-16 

31 

1.690-18 

32 

7.830-11 

33 

3.630-11 

34 

1.690-11 

35 

7.830-12 

36 

3.630-12 

37 

1 .690-12 

38 

7.830-13 

39 

3.630-13 

48 

1.690-13 

Table  3: 


Degree 


2 

3 

2.490488 

4.580488 

1 . 480488 

2.840488 

7.870-81 

1 . 790488 

4.430-81 

1 . 130488 

2.490-81 

7.130-61 

1 .480-81 

4 . 580-81 

7.870-82 

2.840-81 

4.430-82 

1 . 790-81 

2.490-62 

1.130-81 

1 . 480-82 

7.130-82 

7.870-83 

4.580-82 

4.430-63 

2.840-82 

2.490-83 

1 .790-82 

1 . 480-83 

1 .130-82 

7.870-84 

7.130-83 

4.430-64 

4.580-83 

2.490-84 

2.840-63 

1 .480-84 

1 . 790-83 

7.870-85 

1. 130-83 

4.430-85 

7.130-84 

2 . 490-65 

4 . 580-84 

1 . 480-85 

2.840-84 

7.870-86 

1 . 790-84 

4.430-86 

1.130-84 

2.490-86 

7.130-85 

1 . 480-86 

4.560-85 

7.870-87 

2.840-85 

4.430-87 

1.790-85 

2 . 490-87 

1.130-85 

1 .480-87 

7.130-66 

7 . 870-88 

4.580-86 

4.430-88 

2.840-86 

2.490-68 

1 .790-86 

1 . 480-88 

1.130-86 

7.870-89 

7.130-87 

4.430-89 

4.580-67 

2.490-69 

2 . 840-87 

1 . 480-89 

1 . 790-87 

7.870-16 

1.130-87 

4.430-16 

7.130-88 

Values  of  h  for 


Derivative  (k) 


4 

5 

8.250488 

1.180481 

5.620488 

8.490488 

3.630488 

6.110488 

2.610488 

4.480488 

1.780488 

3.170488 

1 .210488 

2.280488 

6.250-81 

1.640488 

5.620-81 

1.180488 

3.830-61 

8.490-81 

2.610-61 

6.110-81 

1.780-81 

4.480-81 

1 .210-81 

3.170-81 

8.250-82 

2.280-81 

5.620-82 

1.640-61 

3.830-82 

1.180-81 

2.610-62 

8.490-82 

1.780-62 

6.110-82 

1 .210-82 

4.480-82 

8.250-83 

3.170-82 

5.620-83 

2.280-82 

3.830-83 

1 . 640-82 

2.610-83 

1.180-82 

1.780-63 

8.490-83 

1 .210-83 

6.110-83 

8.250-84 

4.460-83 

5.620-84 

3.170-83 

3.830-84 

2.260-83 

2.610-64 

1 .640-63 

1 .780-64 

1.180-83 

1.210-64 

8.490-84 

8.250-85 

6.110-84 

5.620-85 

4.460-84 

3.830-85 

3.170-84 

2.610-85 

2.280-84 

1 . 780-85 

1 .640-84 

1 .210-65 

1.180-84 

8.250-86 

8.490-85 

5.620-86 

6.110-85 

3.630-86 

4 . 460-65 

2.610-86 

3.170-65 

moximum  occurocy 

Appendix  I:  Numerical  Differentiation 

1.2  Round-off  versus  Troncotion  Errors 


6 


1 .880481 
1 .410481 
1.860481 
7.930488 
5.9S0488 
4.46D488 
3.340488 
2. 810488 
1 .880488 
1 .410488 
1.860488 
7.93D-81 

5. 950- 81 

4. 460- 81 

3.340- 81 

2.510- 81 
1.880-81 

1.410- 81 
1.860-81 

7.930- 82 

5.950- 82 

4.460- 82 

3.340- 82 

2.510- 82 
1 .880-82 

1.410- 82 
1 .660-82 

7.930- 83 

5.950- 83 

4.460- 83 

3.340- 83 

2.510- 83 
1 .680-83 

1.410- 83 
1 .660-83 

7.930- 84 

5.950- 84 

4.460- 84 

3.340- 64 

2.510- 84 
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l. 3  Dorlvativi  of  Stop  Funetioni 

It  turn*  out  that  vary  amall  jump*  (of  magnitud*  ata,  aoy)  in 
aoffl*  darlvotlv*  can  b*  datactad  only  by  axomining  o  highar  ordar  darivativa 
rathar  than  tha  diaeontinuou*  darlvatlv*  Itaalf.  Wa  datarmina  optimal 
volua*  of  h  and  minimum  voluaa  of  ata  by  Impoaing  two  ragui ramanta.  Firat, 
th*  dalta  tarm  ahould  ba  aignif Icantly  lorgar  than  round-off  arrora,  and, 
aacond.  tha  dalta  tarm  ahould  ba  atgnl f leant ly  largar  than  tha  riaa  of 

(k) 

PHI  over  half  tha  window.  Thia  laoda  to  th*  aquation* 

ar-k  -k 

atoadalto  *h  >  pal  ah  *(14ab*(phi))*tau*olpha 

m. k  2  k 

and 

m-k  (k+1) 

atoadalta  *h  *  pal  ahaPHI 

m.k  3 


wharo  pal  and  pal  arc  aofoty  factor*. 

2  3 

Wa  eonaidar  two  caaaa.  If  m  «  6  (l.a.  thar*  I*  a  atap  in  tha 

trial  function)  aquation  (4)  ia  indapandant  of  h.  Solving  for  ata  yiald* 
th*  minimum  valu* 

ato  >  pai  •(l4aba(PHI))*tau*olpho  /  dalta 

2  k  e.k 

Th*  maximum  valu*  of  h  (abov*  which  th*  dalta  tarm  will  b* 
(k) 

obacurad  by  th*  grawth  In  PHI  )  ia  obtained  by  aolving  (5)  for  h: 

(k+1)  1/(k+1) 

h  ■  [*ta«d*lta  /  (pai  aPHI  )J 
e.k  3 

If  m  >  6  (i.a.  th*  diacontinuity  occur*  in  o  darivativa)  th* 
minimum  volua  of  ata  can  b*  computed  by  aliminating  h  between  (4)  and  (5), 
ond  aolving  for  ato.  Thia  give* 
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(k+1)  m 

•to  -  [(psi  •abs(PHl  )) 

3 


k+1-m  1/(k+1) 

•  (psi  •(1-faba(PHI))*tau*alpha  )  )] 

2  k 


/  delta  m, 
( 


The  corresponding  optimal  value  of  h  (above  which  the  delta  term 

(k) 

is  too  small  to  be  visible  in  the  general  growth  of  PHI  ond  below  which 
it  is  obscured  by  round-off  errors)  con  be  obtained  by  eliminating  eta 
between  (4)  and  (5)  and  solving  for  h,  giving 


1/(k+1) 

h  •  (psi  •(1-»-abs(PHI))*alpho  etou  /  psi  ) 

2  k  3 


Notice  that  this  expression  is  independent  of  m,  the  degree  of  the 
derivative  in  which  the  jump  occurs. 


Table  4  lists  minimum  values  of  eto  (in  the  lost  seven  columns) 
and  the  corresponding  optimal  volue  of  h  (in  the  second  column)  for  the 
case  of  the  exponential,  and  for  a  few  volues  of  tou.  Notice  how  the  size 
of  the  jump  that  con  be  detected  is  sometimes  decreased  by  several  orders 
of  magnitude  If  higher  order  derivatives  ore  exomined. 

too  •  1 .60-16 

Degree  of  discontinuous  derivative 
k  optimal  hO  1  2  3  4  5  6 


0 

1 

2 

3 

4 

5 

6 


5.00-10 

1 .50-09 

4.50-05 

1.20-08 

1.60-03 

1.20-08 

1.20-02 

1 .80-08 

4.80-02 

1.60-08 

1.20-01 

2.00-08 

2.40-01 

1.90-08 

1.30-04 

7.60-06  4.80-03 
2.90-06  4  60-04 
1.00-06  1.10-04 
7.80-07  3.90-05 
4.00-07  1.70-05 


3.70- 02 

6.90-03  1.40-01 
1.60-03  4.60-02 

3.70- 04  1.00-62 


3.50-01 

8.70-02  7.20-01 
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tau  ■  1 .00-17 

Degree  of  dieeont inuous  derivative 
k  optimal  h  0  1  2  3  4  5  6 

0  5.00-17  1.50-16 

1  1.40-08  1.20-154.20-08 

2  7.40-06  1.20-151.60-102.20-05 

3  2.20-04  1.60-15  1.60-11  1.50-07  6.60-04 

4  1.90-03  1.60-15  2.50-12  7.00-09  1.10-05  5.70-03 

5  7.90-03  2.00-15  1.10-12  8.30-10  4.90-07  1.90-04  2.40-02 

6  2.40-02  1.90-154.00-131.70-103.70-081.00-058.70-047.20-02 

tou  -  1 .00-24 

Degree  af  discantinuous  derivative 
k  optimal  h  0  1  2  3  4  5  6 

0  5.00-24  1.50-23 

1  4.50-12  1.20-22  1.30-11 

2  3.40-08  1.20-22  3.50-15  1.00-07 

3  3.90-06  1.80-229.10-174.60-11  1.20-05 

4  7.60-05  1.60-22  6.30-18  4.40-13  1.70-08  2.30-04 

5  5.40-04  2.00-22  1.70-18  1.80-14  1.60-10  8.70-07  1.60-03 

6  2.40-03  1.90-22  4.00-19  1.70-15  3.70-12  1.00-08  8.70-06  7.20-03 

tau  -  1 .00-31 

Degree  of  discontinuoue  derivotive 
k  optimal  h  0  1  2  3  4  5  6 

0  5.00-31  1.50-30 

1  1.40-15  1.20-29  4.20-15 

2  1.60-10  1.20-297.60-204.80-10 

3  7.00-08  1.80-29  5.10-22  1.50-14  2.10-07 

4  3.00-06  1.60-291.60-232.80-172.80-110.10-06 

5  3.70-05  2.00-292.50-243.90-194.90-144.00-091.10-04 

6  2.40-04  1.90-29  4.00-25  1.70-20  3.70-16  1.00-11  8.70-08  7.20-04 


Table  4;  The  euperiority  af  testing  Delta  Functians  and  their  o  vatives 
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Appendix  11: _ Orooniiotion  of  MICROSCOPE 

This  appendix  contains  (1)  a  brief  discussion  of  the 
organ i zot i ono I  framework  upon  which  MICROSCOPE  is  built,  (2)  o  detailed 
explanation  of  the  sampling  scheme  and  logic  which  MICROSCOPE  uses  to 
calculate  the  tangential  derivatives,  and  (3)  a  brief  description  of  the 
common  blocks  used. 

MICROSCOPE  hos  evolved  over  o  period  of  several  years, 
progressing  through  several  preliminary  versions.  The  core  of  the  present 
version  was  designed  and  written  by  Bill  Harris  during  the  summer  of  1983. 
MICROSCOPE  then  grew,  under  the  hands  of  Peter  Alfeld,  to  o  substont iai ly 
larger  pockoge  which  is  the  current  version.  This  explains  lack  of 
uniformity  in  the  commenting  of  the  source  code.  It  also  accounts  for  the 
inconsistency  in  the  notation  used  within  the  program  versus  thot  employed 
in  this  manual.  However,  only  very  serious  users  interested  in 
understanding  and  possibly  manipulating  the  code  (for  whome  this  oppendix 
is  written)  should  suffer  from  such  idiosyncrasies  of  the  pockoge. 

11.1  Qraonizotion  of  the  Progrom  Poekooe 

MICROSCOPE  follows  o  common  and  well  established  format  for 
interactive  programs  in  FORTRAN.  The  user  must  provids  o  short  main 
program,  which  assigns  FORTRAN  Input/Output  device  numbers  end  screen 
parameters,  and  then  colls  a  moster  routine  (MCRSCP) .  MCRSCP  is  o  driver 
program  from  which  oil  other  routines  ore  colled  ond  to  which  execution 
returns  after  each  command  has  been  completed.  This  is  accomplished  by  ( i ) 
prompting  for  a  new  command,  (ii)  reading  and  recognizing  the  new  command, 
(iii)  using  o  computed  GO  TO  to  bronch  to  the  appropriate  routine,  (iv)  a 
bit  of  logic  and  computation  (i.e.  sompling  the  function  and  calculating 
the  derivatives),  (v)  displaying  the  results  in  some  graphic  form,  and  (vi) 
updating  an  options  table  and  returning  to  step  (i). 

To  help  understanding  the  following  flowchart  several  ideos  and 
their  representative  variables  need  to  be  introduced.  The  logical  variable 
LGO  indicates  whether  the  screen  is  active  (LCO  •  .TRUE)  or  inactive  (LGO  « 
.FALSE.).  It  is  set  and  reset  at  o  number  of  ploces  in  the  pockoge.  The 
only  place  where  the  value  of  LGO  hos  an  effect  on  the  progrom  flow  is  in 
the  routine  MCRSCP,  between  statements  7990  and  8300. 

There  are  4  basic  types  of  commands.  These  correspond  to  the 
'nteger  variable  ITYPE  which  is  set  oppropr i ote I y  by  each  command. 
Commands  of  ITYPE  ••  1  imply  that  the  function  must  be  sompled  from  scrotch. 
Since  other  changes  may  be  desired  before  computation  begins,  these 
commands  deoctivate  the  screen  ond  this  is  accomplished  by  setting 
LGO  •  .FALSE..  Commands  which  ore  of  ITYPE  •  2  require  some  computation. 
Two  special  ITYPE  •  2  commands  are  FORCE  and  GO.  These  cause  LGO  to  be  set 
to  .TRUE,  and  ore  the  only  commands  that  activate  an  inactive  screen.  The 
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ramaining  coflMionda  of  ITYPE  *  2  only  moka  modifi cat  Iona  of  tha  currant 
aatup  and  uaually  only  raquira  amoll  omounta  of  odditionoi  computation. 
Whan  LGO  •  .TRUE,  any  ITYPE  ■>  2  commond  ia  axocutad  bafora  tha  naxt  command 
ia  prompted  for.  ITYPE  -  3  commonda  do  not  raquira  computation  but  couaa 
tha  alphonumarical  diaploy  to  ba  altarad.  Pinal ly,  commanda  of  ITYPE  *  4 
do  not  altar  tha  diaplay  and  do  not  raquira  computotion. 

Tha  other  major  logical  control  vorlobla  la  LSCRN.  Ita  main  uaa 
ia  in  acraan-adi t ing  varaiona  of  MICROSCOPE.  LSCRN  ■  .FALSE,  maana  that 
tha  alphonumarical  diaplay  haa  bean  ovararittan  by  text  and  muet  bo 
raganaratad.  LSCRN  ■  .TRUE,  maona  that  tha  diaplay  la  intact  from  tha 
pravioua  command  and  only  aalaetad  pixala  need  to  ba  updotad.  LSCRN  ia 
uaad  by  tha  diaplay  routine  SCRAPH  in  tha  file  SCREEN  . 

Computation  time  for  evaluating  diaplay  functiona  con  ba  aizobla. 
Tha  program  logic  and  aampling  achama  wara  daaignad  to  utillza  pravioua 
eomputotiona  and  allminata  unnacaaaory  racomputot Iona  to  tha  axtant 
poaalbla.  For  axampla,  ahifta  and  mogni f icot ion  chongaa  oftan  raquira  only 
a  few  additional  avqluationa. 

Tha  flowchart  that  followa  givaa  an  overview  of  tha  logicol  flow 
in  tha  MCRSCP  driver  routine.  It  concantrotaa  on  tha  portion  batwaan 
atatamant  1S06  and  8366  whara  tha  logic  io  aomawhot  involvad. 
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Epitome  Elowchort 


from  Main  Progrom 


isee 


7909 


8100 


•  output  messoge  identifying  MICROSCOPE 

•  check  screen  porometers 

*  copy  channel  numbers,  screen  porometere,  and  set  flogs 

*  set  LGO  -  .FALSE. 

*  compute  round-off  characteristics 

•  set  defoult  values  for  some  of  the  variables  contained 

the  various  common  blocks 

input  some  of  the  help  documentation  os  well  os  news 


(NUMOIG) 


(INHELP) 

(DIALOG) 


profflpt/input/porse  commend 
reset  error  flogs 

branch  to  the  appropriate  subroutine  (using  a  computed  GO  TO) 

This  corresponds  to  statements  1600  through  7800.  Many 
of  those  commands  alter  variables,  with  subscript  2,  in 
the  common  block  OPTION.  These  will  then  be  used  by 
later  routines. 

if  command  is  of  ITYPE  •  1  (i.e.  requires  complete  recomputation) 
then  set;  LGO  -  .FALSE, 
if  ITYPE  •  2  and  LGO  -  .TRUE. 

then  colculote  volues  for  NUM  and  DEN  corresponding  to 
sampling  changes.  (See  documentot ion  in  port  2 
of  this  appendix)  (SGAtMA) 

check  for  errors  and  output  appropriate  messages  if  any  occurred 
i f  error  ■  .TRUE. 

- then  go  to  1500 

if  LGO  ■  .FALSE,  and  the  Refresh  Screen  (RSCREEN)  command  is  given 
then  display  the  old  graphical  informotion  (SGRAPH) 

if  LGO  -  .FALSE. 

- then  go  to  1500 

branch  to  the  oppropriote  section  of  the  progrom  depending  on 
the  value  of  ITYPE  (using  a  computed  GO  TO) 

(if  ITYPE  -  2) 

•  check  to  see  which  ospects  of  the  sampling  have  been  oitered 

by  the  command  and  set  the  oppropriote  logical  flogs  (CHKCMP) 
•coll  the  routines  flogged  by  the  CHKCMP  routine. 

(the  octuol  sampling  ond  computation  of  derivatives  occurs 
only  in  the  routine  SCMUPD) 

(if  ITYPE  -  3) 

•  set  up  the  graphics  dato  (i.e.  update  plot  array)  (SGRUPD) 

•  graph  the  doto  (SGRAPH) 

if  logging  is  required 

then  output  a  copy  of  the  grophicol  display  in  scroll  mode 
copy  current  set  of  porometer  values  from  new  to  old.  (OPCOPY) 
go  to  1500 
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11.2  Samplino  Scheme  ond  Logie 


In  the  remoinder  of  this  appendix  an  explonation  of  the  sompling 
scheme  and  computational  logic  is  given.  The  varioble  h,  os  used  in  the 
source  code,  is  not  the  some  os  the  variable  lobeled  os  h  on  output.  To 
reduce  confusion  we  will  uee  the  following  notation  In  this  appendix  (but 
nowhere  else  in  this  monuol): 

Variable  in  Variable  in  Variable  in  Output 

Source  Code  Appendix  and  elsewhere  in 

the  monuol 


h  ss  (step  size)  s 

sw  tw  (stencil  width)  h 


The  output  names  correspond  to  the  common  usage  whi le  the  names 
used  in  the  program  and  appendix  hove  the  meanings  defined  below. 

The  basic  Idea  of  the  sampling  scheme  consists  of  representing 
the  relevant  port  of  the  line  of  investigation  in  two  linear  arrays,  one 
corresponding  to  points  in  the  domain,  and  one  containing  corresponding 
values  of  the  display  function  (If  it  has  been  evaluoted).  Any  changes  in 
the  stencil  width  and  step  size,  os  well  as  small  shifts  of  the  point  of 
e;amination  are  then  implemented  by  shuffling  the  informotion  in  these  two 
arrays.  The  size  of  the  orroys  has  been  chosen  to  accomodate  foiriy  iarge 
changes  in  parameter  values  with  tittle  computation.  The  key  voriobtes  and 
arrays  used  in  the  program  are  as  follows: 
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XS(IXSIZE);  tha  sampling  orray  for  tha  domain  space,  giving  the 
distance  from  the  point  of  examination  in  units  of  tha  direction 
of  invastigotion. 

FS(IXSIZE):  tha  sompiing  orroy  contoining  tha  value  of  the  display 

function  f  at  tha  corresponding  point  in  XS() . 

IXSIZE:  tha  dimension  of  the  XS()  and  FS()  arrays,  (tha  currant 

value  of  IXSIZE  ia  5377} 

ICENTR;  tha  index  of  tha  canter  alamant  of  tha  XS()  and  FS()  arrays. 

(the  currant  value  of  ICENTR  is  2689) 

ROPSTS(2):  this  is  the  "step  sisa"  corresponding  to  8  points  in  tha 

XS()  orray.  (sea  tha  illustration  below)  It  oppaors  in  tha  common 
block  OPTION.  Tha  raoaon  for  introducing  it  hare  is  that  for  most 
computations  it  will  remain  fixed  while  tha  other  poromatars, 
such  os  sw  and  ss,  will  be  doubled,  halved,  etc. 
sa:  (step  siie)  tha  spacing  between  tha  points  plotted  olong  tha 

line  of  investigation.  This  is  h  in  tha  source  coda, 
sw:  (stencil  width)  tha  width  of  tha  interval  over  which  tha 

f ini ta-di f faranca  formulos  ora  calculated.  This  is  h  on  output, 
nh:  ths  spacing  between  plotted  pointa  in  tha  index  space.  It  must 

be  a  positive  integer  multiple  of  4. 
nw:  half  tha  stencil  width  in  index  spoca.  It  must  olso  bs  a 

positive  integer  multiple  of  4.  Tha  number  of  points  in  index 
space  corresponding  to  the  stencil  width  is  1  -f  2»nw. 

The  sampling  paints  are  spaced  unevenly  olong  the  line  of 
investigation,  but  evenly  in  index  spoca.  In  the  following  figures, 
osterisks  denote  the  points  at  which  the  display  function  nseds  to  be 
evaluated.  A  hyphen  marks  a  point  for  which  spocs  for  a  function  volus  hos 
bsen  provided  in  the  arrays  XS  and  FS,  but  where  no  evaluation  is 
necessary.  Consider  a  point  p  at  which  we  wish  to  evaluate  tha  display 
function  ond  its  derivatives  up  to  order  6.  The  stencil  is  centered  obout 
P- 
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(p) 

(1) 

(p) 

(2) 

(P) 

(3) 
(P) 

(4) 
(P) 

(5) 
(P) 

(«) 

(P) 


3  4  8 


8  4  3 


(h  on  output) 


•w 

2 

-> 


Tho  figuro  doocriboo  tho  octuol  oompllng  itrueturo  olong  tho  lino 
of  fnvoot (got ion.  Spoeo  hoo  boon  providod  only  for  pointa  at  which  tho 
diaploy  function  may  hova  to  bo  ovoluotod.  Providing  apoco  for  aompling  ot 
ineromanta  of  ow/12  would  bo  aim'^or  but  woatoful  ainco  ot  moat  only  9 
pointa  (not  13)  aro  noodod  in  ooch  atoncil. 

Noxt,  conaidor  tho  probiom  of  oxtonding  tho  aompling  ovar  on 
ontira  lino  aogmont.  Tho  porticulor  oxtonaion  uaod  in  yiCROSCOPE  woa 
chooan  ao  that  changing  tho  atop  aizo  and  atancil  width  by  powora  ot  two 
would  roquiro  minimal  computation.  Othor  foctora  ora  ollowad  but  roquiro 
grootor  computation. 

Conaldoring  a  portion  of  tho  lino  aogmont  containing  a  aoquonco 
of  throo  plotted  pointa:  p1.p2,  ond  p3.  tho  apocing  (for  tho  ratio  ew/aa  - 
2/1,  nw  >  4,  and  aaauming  all  tongontlol  dorivativoa  aro  to  bo  plotted) 
looka  Ilka; 
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.  • . •  .  . 

pi 

P2 

P3 

• . 

(stencil  for 

pi)- 

-4  -3  -2  -1 

e 

1  2  3 

4 

•  - 

•  .  .  . 

.  .• 

.  .  .  .  • 

(stenei 1 

for  p2). 

-4 

-3  -2  -1 

e 

1 

2  3 

4 

.(stenei 1  for  p3)  *. . 

.  • 

.  • . • 

--4 

-3 

-2  -1 

e 

1  2 

3  4 

< - 

-  9W  - 

- > 

<— 

-  ss  - 

— >  <— 

“  *8  - 

> 

Note  how  the 

stenei  Is 

ove  r 1  op , 

allowing 

the  use  of  one 

value  of 

the  display  funct 

on  for 

more  than 

one  (here  two  or 

three)  values  of  o 

derivative.  If 

the  stencil  width  is  doubled  while  the 

step  size 

is  kept 

constant  (corresponding 

to  an  increase  in 

the 

w i ndow 

width 

by  a  factor  2) 

wo  obtain: 

.  • . *  ,  . 

P1 

P2 

p3 

.  ♦ 

.  .  . 

.  .  .  »  ♦ 

(stenei  1 

for  pi). 

-4  -3 

e 

3 

4 

5 

6 

• 

0 . 

*  .0.  * 

. • 

-8  -5 

-4 

-3 

6 

3 

4 

5 

6 

.(stenei 1  for  p3)  *. . . 

........ 

.  • 

— .  . 

.  • . •  .  . 

-5 

-4 

-3 

e 

3  4 

<_ -  - > 

< -  •»  - >  < -  ss  - > 


Tha  obova  oparatlon  roquiraa  additional  avaluotion  of  the  dieplay 
function  only  at  the  boundary  of  the  interval  of  inveetigation.  The 
reverse  process  of  halving  the  window  width  would  require  no  additional 
evaluations  at  all. 

So  that  all  derivatives,  up  to  and  including  the  sixth,  can  be 
evaluated,  the  minimal  number  of  points  betwsen  plotted  points  (such  as 
p1,p2,p3)  is  4.  Furthermore,  the  distance  is  alwoys  a  multiple  of  4. 
Similarly,  the  minimal  number  of  points  nseded  in  o  stencil  is  9  (i.e. 

-4 . 0 . 4),  and  it  is  always  of  the  form  l^-Can  integer  multiple  of  6), 

although  evaluation  takes  piece  ot  at  most  9  points.  These  features  are 
illustrated  in  the  following  progression  of  increasing  stencils: 
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...... 

.  *  . 

P2 

•  .  . 

.  .  .  • 

.  • 

.  •  . 

.  .  • 

-4 

-3 

-2 

-1 

9 

1 

2 

3 

4 

• .  . 

.  .  .  • 

...... 

.  .  • 

-8 

-5 

-4 

-3 

9 

3 

4 

5 

8 

...... 

-8 

-6 

-4 

9 

4 

6 

8 

Thooo  3  ttoncilo  span  tndox  aoto  of  9.17.  ond  25  points, 
roapoetivoly  (although  tha  third  aat  would  not  fit  on  thia  poga).  Tha 
varlobla  nw  for  thasa  la  4.8.  ond  12.  raapaetivaly. 

Tha  ralotion  batwaan  tha  guontltlaa  aa  and  R0PSTS(2)  ia  givan  by 
nh 

aa  -  —  •  ROPSTS(2) 

8 

Tha  program  atrivas  to  maintain  a  valua  of  nh  ■  8.  For  a  mora 
datoilad  diacuaaion  aaa  tha  daaeription  of  tha  routino  SCAtitA  baiow. 

Tha  raiation  batwaan  tha  point  of  axominotion.  P.  o  point  of 
avoluotion.  PS.  tha  distanea  t  of  pa  from  P.  tha  diraetion  of  axomlnatlon. 
d.  and  tha  quantitiaa  in  MICROSCOPE,  ia  givan  by; 


XS(ICENTR-2)  • 
XS(ICENTR-I)  • 
XS(ICENTR)  - 
XS(ICENTR+1)  - 
XS(ICENTRi'2)  • 
XS(ICENTR+3)  - 
XS(ICENTR+4)  - 
XS(ICENTR+5)  - 
XSdCENTR-^)  - 
XS(ICENTR+7)  - 
XS(ICENTR+8)  - 


-ROPSTS(2)/4 

-ROPSTS(2)/6 

9 

ROPSTS(2)/6 

R0PSTS(2)/4 

ROPSTS(2)/3 

R0PSTS(2)/2 

2*ROPSTS(2)/3 

3»R0PSTS(2)/4 

5*R0PSTS(2)/e 

R0PSTS(2) 


and 
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P  -  R0PPNT(»,2) 
d  -  R0PUDI(*.2) 
t  -  XS(i) 

PS  -  P  +  t  •  d 
FS(i)  -  f(PS) 

LOEF(I)  -  .TRUE.  <— >  FS(1)  hoa  baan  avoluatad 

To  axpiain  tha  calculation  of  tha  darivotivaa,  conaidar  a  point  p 
at  which  wa  wiah  to  avaluota  tha  function  and  ita  darivotivaa.  To  avoluota 
all  darivotivaa  up  to  tha  aixth  wa  naad  9  pointa.  Wa  nunbar  thaia  oa 
fol Iowa: 

1  234  5  678  9 

• . . * . . • 

p  ... 

aw  aw  aw  aw  0  aw  aw  aw  aw 


2346  643  2 

Inataod  of  uaing  tha  formulos  givan  in  aaction  1.2  diractly,  tha 
following  aqui volant  but  mora  afficiant  formuloa  ora  uaad: 


gO  ■ 

f5 

g1 

f1  + 

f9 

g2  • 

- 

( 

f1  - 

19 

g3  B 

4 

• 

( 

f3  + 

17 

g4  ■ 

4 

0 

( 

f3  - 

17 

g5  m 

4 

0 

( 

f2  - 

f8 

g6  <• 

-  6 

0 

( 

f2  4 

18 

) 

) 

) 

)  +  5  •  (  f6  -  f4  ) 

)  +  15  •  (  f4  +  f6  ) 


(«) 

0  “  gO 


(1) 

0  -  (1/aw)  •  g2 
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(2)  2 

Q  «  (2/ow)  •  (  g1  -  2»g0  ) 

i 


(3)  3 

0  -  (4/aw)  /4  •  (  2*g2  +  g4  ) 

i 


(4)  4 

0  «  (4/aw)  •  (  g1  -  g3  +  6*g0  ) 

i 


(5)  5 

0  -  (6/aw)  /2  .  (  g2  4  g5  ) 

I 

(6)  6 

0  ■  (8/aw)  •  (  g1  +  g6  -  20*g0  ) 


Tha  aaguanea  In  which  MCRSCP  col  la  tho  laoin  routinaa  ia: 

SGAMMA  -  coleulotaa  gamma  (ao#  following  diaeuaaion) 

CHKCMP  -  ehaeka  to  aaa  whot  chongaa  in  eomputotiena  ora  naadad 
SSMPUD  -  ahift  and  copy 
SCMUPD  -  compute  now  voluoa 

Whan  0  atop  aiza/atoncil  width  chonging  command  ia  given,  tha 
program  will  ottompt  to  change  nh  and  nw  without  having  to  chonge  R0PSTS(2) 
(ond  hence  copy  or  throw  out  data).  Otherwiee,  MCRSCP  ottempte  to  chooee  a 
change  in  ecoling  (R0PSTS(2))  which  eovee  oe  much  dato  ae  poeeible.  Thie 
ia  the  purpoae  of  the  routine  SGAInNA  which  we  now  diecuae  in  eome  detail. 
Firet  we  introduce  additional  nototion. 


width:  the  number  of  columne  on  the  terminal  ecreen  ueed  for  plotting 
m:  ecoling  factor  for  etep  eize  (either  a  poeitive  integer  or  1/(poaitive 

integer)  ) 

v:  ecoling  factor  for  etencil  width  (or  window),  (either  poeitive  integer 

or  1/(poeitive  integer)  ) 

rho:  ioweet  common  multiple  of  the  denominotore  of  m  and  v. 
e:  ai  R0PSTS(2).  (boee  etep  aize  when  nh  «  8),  (thie  dif fere  from  the  a 

thot  denotea  the  interval  elaewhere  in  thie  monuol) 

Primee  will  denote  the  new  voluee  for  nh,nw,ae,and  e.  For 
example:  h'  denetee  the  etep  eize  after  chongee  hove  been  made. 

Define:  gommo  •  e/e'  *  num/den  >  0 


Appendix  II:  Orgonizotion  of  MICROSCOPE 

1 1. 2  Sompling  Scheme  ond  Logic 


page  101 


MICROSCOPE: 
Usar's  Manual 


A  Softwara  Syatem  for  Multivoriota  Analyaia  poga  102 

Varaion  of  April  1984 


Conatrainta:  Tha  conMonda  aet  the  voluaa  of  at  and  v.  Thera  ora  S 
conatrainta  which  are  then  used  to  chooaa  tha  volua  of  gamiaa.  They  ora  not 
all  independent  of  each  other. 


(1) 

4  <■  nh ‘ 

(2) 

4  <«  nw  * 

(3) 

nh‘/4  must  be 

(4) 

nw‘/4  muat  be 

(5) 

effective  aompling  width 

in  the  index  apace 

We  hove: 

as  -  (nh/8)  *  a  i 

Al  ao: 

aa' 

Hence : 

aa‘  «  (nh’/8)»a' 

— > 

(nh'/O)  •  •* 

— > 

nh*  ■  m  •  nh  •  a/a' 

Simi lor ly : 

— > 

aw'  -  (nw'/4)  •  a' 

— » 

ixsi^t/wiatn 


aw 


gomaa 


In  aumaory; 


IXSIZE 


nh'  •  la  •  nh  •  gomao 
nw'  «  V  •  nw  •  gonvao 

Referring  to  the  conditiona  numbered  above,  divide  (1)  by  m  •  nh  and  (2)  by 
V  •  nw  to  get: 


la  •  nh 


<•  gamma  <• 


IXSIZE 

width 


1 


nh 


ond 
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V  •  nw 


<m  goiMO  <- 


Lot 


4  4 

a  ■  niax(  -  ,  -  ) 

m  •  nh  V  •  nw 


IXSIZE  1 

■  •  - "  ' 

2  V  •  nw 


and 

IXSIZE  1  IXSIZE  1 

b  »  min(  -  •  -  ,  -  •  -  ) 

width  Ik  •  nh  2  v  •  nw 

ProA  obova,  gonmo  muat  aotiafy: 

a  <<•  gomko  <>  b 

Furtharikoro; 

nh'  ■  m  •  nh  *  gomno  ond  mod(nh*,4)  ■■  0 

■-•>  (•)  mod(  A  •  nh  •  gonMO  ,  4  )  •  0 

nw’  >  V  •  nw  *  gamma  ond  Acd(  v  •  nw  •  gamma  ,  4  )  •  0 

■«>  (••)  mod(  V  •  nw  •  gommo  ,  4  )  •  0 

Th«  rout  in*  SGAMilA  u«e«  tho  obova  two  conditiena  (•)  ond  (••)  to 
chooaa  occaptobla  voluaa  for  gommo  (ond  hanca  num  and  dan  ainca  goniao  » 
num/dan) . 
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Outline  of  SGAMmIA 

•  check  to  eee  if  the  changee  m  and  v  con  be  occomodated  by  changing 

nh  and  nw,  but  not  ae  (h  in  code). 

•  if  no  changes  in  as  are  needed 

<  -  then  num  ■  den  ■  1  •“>  gamma  •  1  and  return 

•  calculate  o  and  b 

•  i  f  I*  "  V 

<  -  then  set  gamma  -  m,  R0PSTS(2)  ■  m  •  R0PSTS(2),  ond  return 

•  calculate  rho 

m'  w  rho  •  m  (m*  ond  v*  ora  poaitiva  integers) 
v’  ■  rho  •  w 

•  determine  "all"  possible  numerators  and  denominators  for  gommo 

(up  tp  32)  which  satisfy: 

a  <m  num/den  <••  b 
mod(  ffl  •  nh  •  gommo  ,  4  )  •  0 
mod(  V  •  nw  •  gamma  ,  4  )  >  0 

«  rank  the  condidate  values  of  gommo  according  to  4  criteria 

(1)  nh*  •  width  +  2  •  nw'  <“  IXSIZE 

(2)  8  •  width  <m  nh'  •  width  +  2  •  nw'  <■  ICENTR 
(this  allows  for  halving  and  doubling  without  rescoling) 

(3)  gamma  a  power  of  2  (positive  or  negative) 

(4)  nh'  >  4  and  nw'  >  4 

•  choose  the  combination  of  num  and  den  which  hoe  the  highest  cor¬ 

responding  rank  and  set  the  variables  occordingly. 

< — •  return 

11.3  Cnmninn  Blocks 

Here  we  list  the  COkMON  blocks  used  in  MICROSCOPE.  This 
information  may  be  useful  in  modifying  the  code  ond  in  occesaing  the  dota 
structures  of  the  package  via  the  command  USER  (see  section  4.15).  For 
each  common  block  we  list  in  parentheses  the  number  of  words  (i.e.  2  words 

for  each  DOUBLE  precision  variable  and  1  word  for  eoch  REAL.  INTEGER,  or 
LOGICAL  variable)  occupied  by  the  block.  The  blocks  ore  listed  in  the  smne 
sequence  os  they  appear  in  the  routine  MCRSCP. 


CB  (16)  controls  cross  differentiation. 


CBOSW 

CB0(3) 

CB0U(3) 

ICBO 

LCBD 


width  of  differentiation  stencil 
normoliied  cross  direction 
unnormolized  cross  direction 
order  of  cross  differentiation 
•TRUE,  if  cross  differentiation 


i  s  on . 


otherwise 


.FALSE. 
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IQ  (5)  contain*  the 
RECORD 

GRAPHO 

HELPD 

RSTRTD 

INPUTO 


input/output  d*vic*  number*. 

Record  device  number.  (Output  data  file  for  recording 
ecreen  imogea  and  comment*) 

Graphic*  device  number,  (terminal  or  graphic*  di*play) 
Help  file  device  number. 

Restart  file  device  number,  (unformatted  output  data 
file  for  ua*  in  restarting  at  o  loter  dote) 

Input  device  number,  (usually  the  terminal) 


iOfi  (1)  cant  role  logging 

LCHN  Logging  device  number  (*  0  (ogling  Is  off 


LQGCQM  (5391)  contain* 
L0F(7) 

LPLT(7) 


LDEF(5377)  : 


logical  arrays  used  primarily  by  SCMUPO 
to  l(**p  trock  of  what  hoe  alreody  been  computed. 

Denotes  whether  the  (i-1)‘th  derivative  ha*  been  calcu¬ 
lated.  LDF(i)  m  .false,  means  the  (i-1)*th  derivative 
has  not  been  calculoted.  LOF(t)  >  .TRUE,  mean*  it  has. 

Denotes  whether  the  (i-l)'th  derivotiv*  is  currently 
being  displayed.  LPLT(i)  ■  .FALSE,  mean*  the  (i-l)’th 
derivative  I*  not  being  disployed  and  LPLT(i)  ■  .TRUE, 
mean*  it  is. 

Denotes  whether  or  not  the  interpolot ion  function  ho* 
been  eolcuieted  at  the  coordinates  corresponding  to 
XS(I).  (i.e.  whether  FS(i)  has  been  ossigned) 

LOEF(i)  >  .FALSE,  mean*  the  i'th  value  has  not  been 
sampled  and  LOEF(i)  >  .TRUE,  meon*  it  ho*. 


FUNCQM  (23,426)  conto 
XS(5377) 

FS(5377) 


DF(135,7) 


DFMMX(2.7): 


n*  array*  pertaining  to  the  values  of  the 

interpolation  function  ot  it*  sompied  point*. 

The  dietonc*  array,  olwoys  centered  about  0.  It*  defini¬ 
tion  and  usage  i*  given  in  the  second  port  of  this 
appendix. 

The  array  containing  the  value*  of  the  interpolation 
function  corresponding  to  the  points  specified  by  the 
XS()  array.  See  the  second  port  of  this  oppendix  for 
discussion. 

The  finite  differences.  DF(l,j)  corresponds  to  the  i'th 
plotted  point  and  (j-l)’th  derivotiv*.  These  ore  the 
p1,p2,p3's  of  the  discussion  in  the  second  port  of 
this  appendix. 

These  are  the  minimums/maximum*  of  the  finite  difference* 
In  the  0F()  orroy. 
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ERRCOM  (2)  uaad  to  hondlt  error  checking  and  reporting. 

When  errors  ore  detected  in  any  of  the  sub  routines 
or  cofliputot ion  routines,  the  errors  ore  ftogged  ond 
then  reported  by  the  CHKERR  routine. 

ERRCOO  This  is  on  integer  code  number  for  the  type  of  error 

detected. 

ERR  This  is  the  logicol  flog  which  denotes  whether  or  not 

on  error  hoe  occurred.  ERR  •  .FALSE,  means  no  error 
hos  occurred  ond  ERR  •  .TRUE,  means  on  error  has 
occurred. 


PLTCQM  (16342)  contoins  display  versions  of  the  derivatives. 

IPLOT(i,j)  is  the  "screen  coordinate”  (integer)  of  the 
(j-1)‘th  derivotive  in  the  i*th  column. 

The  previous  screen  image.  ISCRNI(i.J)  contains  the 
Hollerith  representation  of  the  chorocter  to  be 
displayed  ot  the  (i,j)’th  pixel  on  the  screen. 

(where  indexing  starts  at  the  upper  left  corner 
of  the  screen),  i  is  the  horizontal  component  and 
j  is  the  vertical  component. 

The  current  screen  image.  This  has  the  some  description 
as  ISCRN1(,)  obove. 

The  scale  factors  used  to  individually  scale,  up  or  down, 
the  finite  difference  "curves"  (in  the  0F(,)  array) 
so  thot  the  Curve  fills  up  the  entire  screen  in  the 
sense  thot  the  maximum  point  lies  on  the  top  row  ond 
minimum  point  lies  on  the  bottom  row  of  the  screen. 


IPL0T(135,7): 
ISCRN1 (135,57); 


ISCRN2(135.57); 

SCALE(7) 


LOADIO  (3)  contains  channel  numbers  of  the  loading  devices. 


PLOAO 

DLOAO 

CLOAO 


point  loading  device 
direction  loading  device 
cross  direction  looding  device 
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SCREEN  (7)  containa 

OUTPUT 

LINES 

WIDTH 

ILP 

IPRMPT 

lOSPLA 

LSCRN 


tha  valuta  of  poromatara  daacribing  tha 

davica  on  which  moat  intaroction  occura  (i.t.  tha 
acraan  of  tha  tarminol  In  moat  coaaa) 

Tha  davica  numbar  of  tha  grophica  davica  (uauolly  tha 
tarminol ) 

Tha  numbar  of  uaabla  llnaa  on  tha  acraan.  Nota  that 
thta  moy  hova  to  ba  ona  laaa  than  tha  numbar  of 
diapiayabla  linaa. 

Tha  numbar  of  columna  to  ba  diaployad.  It  la  boat  to 
chooaa  on  odd  numbar  for  thia  ao  that  tha  cantar 
can  ba  plottod.  It  muat  ba  laaa  thon  or  aqual  to 
135.  75  la  0  common  valua. 

Tha  numbar  of  linaa  in  tha  curve  plotting  portion  of  tha 
diaploy.  It  muat  ba  laaa  thon  or  equal  to  57. 

Tha  numbar  of  linaa  obova  tha  acraan  bottom  at  which 
tha  prompting  commonda  ora  givan.  It  muat  ba  at 
laaat  2. 

Tha  number  of  linaa  above  tha  acraan  bottom  ot  which  tha 
numaricol  doto  region  bagina.  It  muat  ba  greater  than 
IPRMPT  by  6  or  mora. 

Tha  logical  voriobla  denoting  whether  tha  curvaa  ora 
currently  being  diaployad  or  they  hova  bean  over¬ 
written  by  text.  Thia  ia  to  ollow  tha  program  to 
toka  odvontoga  of  tarminola  which  allow  acraan 
editing.  LSCRN  *  .FALSE,  maona  that  tha  acraan  hoa 
bean  ovarwrittan  ainca  tha  pravioue  command  and 
LSCRN*  .FALSE,  maana  thot  tha  acraan  ia  intoct  from 
tha  pravioua  commond. 


OPTION  (122)  contoina  moat  of  the  vorioblaa,  other  than  tha 

aompling  orraya  in  FUNCOM,  which  moy  chonga  due  to 
commonda.  Thera  ora  2  copiaa  of  aoch  variable. 

(1)  ia  the  pravioua  volua  ond  (2)  ia  tha  currant 
volua.  If  miatokea  ora  mode  in  the  inactive  (daoctiv) 
mode  (i.a.  LGO  *  .FALSE.)  than  tha  UNDO  command  will 
aroaa  tha  chongea  ond  copy  tha  old  aot  (1)  to  (2). 

Tha  vector  opacifying  tha  direction  of  i nvaat i got i on . 

Nota  that  for  oil  vector  valued  quontitiaa,  the 
dimanaion  ia  apacifiad  by  IOPDM()  below. 

Tha  vector  opacifying  tha  canter  of  tha  line  of 

i nvaat i got  ion.  If  the  ahift  ia  6  (i.a.  IOPSH(2)  *  6) 
then  thia  ia  tha  point  which  ia  ot  tha  center  of  tha 
acraan. 

Currently  unuaad  vector. 

Currently  unuaad  vector. 

Tha  nominal  apocing  between  the  pointa  at  which  tha 
finite  diffarancea  are  calculated.  Thaaa  will  ba 
the  plotted  pointa.  '’Nominol"  refara  to  the  baaa 
atap  aiza  when  nh(2)  *  6,  aa  diacuaaad  in  tha  aacond 


ROPDI(3.2) 


ROPPNT(3.2): 


R0PDR1(3,2) 

R0P0R2(3.2) 

R0PSTS(2) 
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KPSVM(2) 

ROPUOI(3.2) : 

NH(2) 

NW(2) 

REFT(2) 

IR1GHT(2) 

I0P0M(2) 

I0PSH(2) 

10PSS{2) 

10PWN(2) 

OPDC(2) 

0PDS(2) 

OPOX(2) 

0PCMP(2) 

0PSPL(2) 
OPOF1(7.2)  : 

0PDF2(7.2)  ; 


port  of  thie  oppendix.  Ae  noted  eorlier,  this  it 
not  the  h  which  oppeore  on  output. 

The  etencil  width  (ew).  Thie  correcponde  to  the  h  which 
appears  on  output.  It  is  discussed  in  the  second 
port  of  this  oppendix. 

The  noritialized  (unit)  vector  specifying  the  direction 
of  invest igot ion.  This  is  the  normalized  version 
of  the  ROPDlO  vector. 

The  step  size  in  index  space.  See  the  discussion  in  the 
second  part  of  this  appendix. 

Half  of  the  stencil  width  in  index  spoce.  See  the 
discussion  in  the  second  port  of  this  oppendix. 

The  index  of  the  left-most  element  of  the  FS()  array 
which  has  been  assigned  o  value. 

The  index  of  the  right-most  element  of  the  FS()  array 
which  has  been  assigned  a  value. 

The  dimension  of  the  domoin  space.  It  con  be  1,2  or  3. 

The  shift,  in  integer  increments,  of  the  plotted  points 
from  the  center  defined  by  ROPSTS(*,2). 

The  scaling  factor  for  the  step  size.  It  is  sst  by 
the  commands. 

The  scaling  factor  for  the  stencil  width.  This 

multiplies/divides  the  current  window,  os  it  appears 
on  the  screen.  It  is  set  by  the  commonds. 

The  logical  flog  determining  whether  or  not  the  center 
of  the  screen  is  to  be  marked. 

The  logieol  flog  determining  whether  or  not  o  scale  is 
to  be  displayed  along  the  horizontal  oxis. 

The  logical  flog  determining  whether  or  not  o  line  of 
dashes  will  be  disployed  to  represent  the  x-axis. 

The  logical  flog  denoting  whether  or  not  any  computotions 
hove  yet  been  done  for  the  current  set  of  ROPPNT  and 
ROPDI . 

Currently  not  in  use. 

The  logical  array  specifying  which  derivotives  ore  to 
be  displayed.  It  is  set  by  the  DGRAPH  and  ECRAPH 
commonds . 

The  logical  array  specifying  which  derivotives  ore  to 
to  be  accentuated  os  they  ore  displayed.  It  is  set 
by  the  ACCENT  commond. 


NRML2E  (1)  controls  normalization  of  derivotives 

NORMAL  ■  .TRUE.  <"“>  normalization  is  on 
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ROOM  (2)  controls  whathar  ona  line  of  the  grophicol 

diaploy  is  tronsferred  to  the  numerical  display. 
This  is  only  necessary  if  all  six  tongential 
derivatives  of  o  cross  derivative  are  plotted. 
ILPUSR  User  specified  space  for  the  graphical  display 

IDSUSR  User  specified  spoce  for  the  numerical  disp  lay 

The  correaponding  variobles  ILP  ond  IDSPLA  in 
the  common  block  SCREEN  ore  modified  as 
necessary. 

HELPER  (7330)  contains  help  information 


HELP 

JHELP1 

JHELP2 

JHELP3 

JHELP(09.2) 


IHELP(72,99)-.  (•,!)  prompt  text 

USER  (5)  see  appendix  III 
PLOMN  (296)  contains  currant  PLOT  options 


Help  device  chonnel  number 
Length  of  the  HSUMRY  information  (in  lines) 
Length  of  the  HELP  information  (in  lines) 
Number  of  available  commonda 

(1.1)  I  ins  in  which  datoiled  help  for  I-th 
command  begins 

(1.2)  length  of  detailed  help 
(*,1)  prompt  text  for  1-th  commond 


1TITLE(72) 

1B0TTM(72) 

NUMBR 

BL(20) 

BR(2e) 

Title 

Legend 
poge  number 
left  bottom  label 
right  bottom  label 

FRAME 

•  .TRUE. 

<— ■> 

frome  is  drown 

COLOR 

-  .TRUE. 

<— > 

plot  is  in  color 

NUMRCL 

-  .TRUE. 

<— > 

numericol  disploy  is 

LMARK 

-  .TRUE. 

<— > 

gropha  ore  marked 

LBLS 

-  .TRUE. 

<■—> 

grophs  ore  labeled 

DATE 

-  .TRUE. 

<•— > 

dote  is  drawn 

TIME 

-  .TRUE. 

<— > 

time  is  drown 

(1)  controls  revaluation 

after  the 

command  FORCE  has 

LFO 

been  issued. 

.  .TRUE.  <— > 

reevol uats 
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AooandiK  111:  Tha  Taat  Poekooa 

Tha  aaomplaa  in  thia  nonuol  con  ba  raproducad,  and  axparianca 
laith  MICROSCOPE  can  ba  gainad,  by  uaing  tha  taat  routinaa  (a  apacial  choica 
of  tha  trial  function  F  and  tha  uaar  intarvantion  routina  SUBUSR,  plua  a 
faw  aupport  routinaa)  auppliad  with  tha  pockoga.  Tha  taat  coda  ia 
contoinad  in  tha  filaa  USRLC  and  USRUC  (aaa  choptar  5). 

Savaral  diffarant  trial  functiona  ora  aval  labia  for  axominot ion, 
and  rounding  to  any  numbar  of  digita  up  to  tha  numbar  aupportad  on  your 
local  inatal lotion  con  ba  aimulotad. 

All  optiona  ora  aalactad  uaing  tha  coomond  USER,  which  puta  you 
into  uaar  noda.  You  con  changa  optiona  aiiailorly  aa  in  tha  PLOT  moda  until 
you  ora  aotiafiad,  and  than  raturn  to  tha  coiMond  moda. 

Optiona  ora  choaan  by  poaaing  numbara  to  tha  progron  according  to 
tha  prompta  givan.  A  currant  lict  of  optiona  ia  mointoinad  on  tha  acraan. 

Tha  trial  function  ia  aalactad  by  giving  ita  rafaranca  numbar. 
Tha  poromatar  ata  ia  1  by  dafault  and  con  ba  chongad  by  typing  -3.  By 
dafault  rounding  to  10  digita  ia  aimulotad.  Rounding  con  ba  chongad  by 
typing  -1  ond  giving  tha  daairad  numbar  of  aignificont  digita.  It  can  ba 
turnad  off  by  giving  0  aa  tha  nixitbar  of  digita.  Tha  following  trtol 
functiona  ora  aval  labia: 


X 


1: 

F(x)  - 

a(0.ata,x) 

or 

F(x)  - 

a(0,ato,x) 

•f  e 

2: 

FCk)  - 

a(1 ,ata,x) 

or 

F(x)  - 

s(1 .ato.x) 

X 

+  0 

3: 

F(x)  - 

a(2,eta,x) 

or 

F(x)  - 

8(2.ato,x) 

X 

+  0 

4: 

F(x)  - 

a(3,eta,x) 

or 

F(x)  - 

a(3,ata, x) 

X 

+  9 

5: 

F(x)  - 

a(4,ata,x) 

or 

F(x)  - 

a(4, ato.x) 

X 

+  9 

6: 

F(x)  - 

a(5,ata,x) 

or 

F(x)  - 

a(S, ato.x) 

X 

+  9 

Tha  addition  of  tha  exponential  in  tha  above  functiona  ia  off  by 
dafault.  and  ia  toggled  on  or  off  by  typing  -2. 

ota*x 

7:  F(x)  -  a 
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2 


6: 

r(x,y)  ■  oto*oba(x)*x  ♦  (1-ata)*obs(y)*y 

9: 

2 

f((».y)  ■  *  •obs(x)*y»aba(y) 

19: 

F(x)  *  tha  apllna  Intarpolant 

X 

to  a  used  in  section  2 

11: 

X 

F(x)  •  a  -  function  10 

(i.a.  F(k)  Is  tha  error  in 

the  spline  19) 

12: 

F(x)  -  9 

13: 

F(x)  m  ata*x 

14: 

ata 

F(x)  -  X 

IS: 

ata 

F(x)  -  X 

Tha  uaar  laopa  la  laft  by  typinp  laro  (or  Juat  RETURN). 

Tha  eoiMionda  STORE  and  RESTART  ara  affactiva  on  tha  taat 
paramatara  which  ara  atarad  In  tha  COMMON  block  USER.  Tha  poraiaatara  hova 
tha  following  maanlnga: 

ETA:  DOUBLE  PRECISION  variabla  containing  tha  valua  of  ata. 

IROUNO;  Numbar  of  digita  to  which  function  voluaa  ora  roundad  (IROUND  -  6 
no  rounding  la  taking  ploca). 

N:  rafaranca  nuiabar  of  tha  trial  function. 

ADO:  LOGICAL  variabla  that  Indicotaa  whathar  on  axponantlol  tarn  la  to  ba 
oddad  to  tha  trial  functiona  dafinad  by  N>  1,...,6. 

Tha  following  tabia  givaa  tha  USER  dafoulta  and  tha  poronatar 
aattinga  nacaaaory  to  ganarata  tha  figuraa  in  thia  manual.  N  is  tha 
rafaranca  nunbar  of  tha  trial  function,  0  Is  tha  numbar  of  digits  roundad 

to.  Tha  antry  " - "  naans  "doss  not  apply”.  Tha  dafoult  sattlnga  ara 

affactiva  on  tha  first  USER  and  aftar  giving  tha  command  SETOP.  Other 
poramaiars  nscassory  for  reproducing  the  axonplas  in  this  manual  (point, 
direction,  dimension.  Interval  ate.)  con  ba  read  off  the  a  I phonuma r I co I 
display  In  each  figure. 
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Figure 

N 

D 

ETA 

exponent i a  1 

Defoul t 

1 

10 

1 

NO 

1 

9 

10 

— 

— 

2-6 

11 

10 

— 

— 

7-8 

10 

10 

— 

— 

9 

12 

10 

— 

— 

10 

13 

10 

2D-9 

— 

11-12 

10 

10 

— 

— 

15-18 

3 

10 

-5D-3 

YES 

19 

1 

10 

NO 

20 

2 

10 

-t 

NO 

21 

3 

10 

-1 

NO 

22 

4 

10 

-1 

NO 

23 

5 

10 

-1 

NO 

24 

6 

10 

-1 

NO 

25-26 

7 

10 

1 

— 

27 

9 

10 

— 

— 

28 

14 

10 

1.5 

— 

29 

15 

10 

1.5 

— 

Table  of  test  parameters 

Because  of  the  random  nature  of  round-off  errors,  somewhot 
different  displays  may  be  obtained  even  with  the  above  parometer  settings. 

Some  of  the  above  test  functions  so  for  hove  not  been  used  in 
this  manual.  They  illustrote  more  subtle  phenomena  thot  will  not  normolly 
occur  in  multivoriote  interpolation  end  opproximot ion.  In  a  typical 
applicotion,  the  trial  function  will  hove  a  finite  degree  of 
differentiability  across  the  front,  and  be  infinitely  often  differentiable 
in  any  direction  contained  in  the  front.  Function  8  i I  lust  rotes  a 
situation  where  the  degree  of  differentiability  is  finite  in  ony  direction, 
but  is  dependent  upon  the  direction  of  i nvest i got i on .  The  function  is  once 
differentiable  in  the  y  direction  ond  twice  di f f erent i obi e  in  the  x 
direction.  The  relative  sizes  of  the  jumps  in  the  oppropriote  derivotives 
can  be  controlled  by  the  parameter  eta. 

The  functions  14  and  15  con  be  used  to  illustrate  discontinuities 
other  thon  step  discontinuities.  The  functions  differ  in  the  woy  they  are 
defined  for  negotive  values  of  x,  for  which  the  volue  of  function  14  is 
negative,  and  thot  of  function  15  positive,  if  the  integer  port  of  abs(eta} 
is  even,  and  of  the  opposite  sign  if  it  is  odd. 

For  both  functions  14  and  15,  the  value  at  x  -  0  is  0  if  eta  is 
negative.  This  prevents  a  floating  point  overflow  and  a  meaningless 
display. 
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Figure  29  Another  3/2  times  differentiable  function 


Figure  29  shows  a  cusp  in  the  first  derivative  whereas  in  figure 
28  the  first  derivative  passes  verticolly  through  the  origin.  Both 
derivatives  ore  continuous  at  the  point  of  exominotion  ond  hove  a  vertical 
tangent  there.  The  foct  that  the  second  dsrivotive  has  a  pole  can  be 
verified  by  observing  that  the  extreme  function  value  increases  os  the 
discretization  porameter  is  decreased.  For  a  cusp,  the  extreme  function 
volue  approaches  a  limit  os  h  decreoses.  This  phenomenon  con  be  utilized 
to  distinguish  between  a  cusp  and  a  pole. 

It  is  of  course  st roightforword  to  program  other  or  odditionol 
examples  for  a  trial  function  by  modifying  the  routines  in  the  teet  package 
oppropr lately. 
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Point  =0.0  s  =  5.00E-03 

Direction  =  1.000  h  =  0.0300 

FO  •  (-0.03332,  0.06764)  Number  of  Evaluations:  249 
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Point  =  0.0  s  =  2.50E-03 

Direction  =  1.000  h  =  0.0150 

FO  •  (-0.02027,  0.02884)  Number  of  Evoluotions:  287 
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Point  =  2.00  s  =  5.00E-03 

Direction  =  1.000  h  =  0.0300 

FO  •  (2.638E-11,  0.01415)  Number  of  Evaluations:  362 
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Point  =  2.00 

Direction  =  1 .000 

FO  •  (2. 638E-11.  0.01415) 

Number  of  Evoluotions;  374 


s  =  5.00E-03 
h  =  0.0300 

FI  +  (-0.061 11.  0.1768) 


This  plot  corresponding  to  figure 
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Point  =  2.00  s  =  2.50E-03 

Direction  =  1.000  h  =  0.0150 

FO  •  (2.638E-11,  3.029E-03)  FI  +(-0.04209.0.07089) 

Number  of  Evaluations:  418 


This  plot  corresponding  to  figure 
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Point  =  0.0 
Direction  =  1 .000 
FO  •(0.8828,  1.117) 

F2  ♦  (-0.2526,  0.2526) 


s  =  2.50E-03 
h  =  0.0150 

FI  +  (1.263,  1.274) 

Number  of  Evoluotions:  505 


This  plot  corresponding  to  figure 
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Point  =  1.000 
Direction  =  1.000 
FO  •  (2.486.  2.973) 

F2  *  (2.478.  3.106) 
Number  of  Evaluations:  592 


s  =  2.50E-03 
h  =  0.0150 

FI  +  (2.387.  2.898) 

F3  0  (2.730.  4.058) 


This  plot  corresponding  to  figure 
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Point  =  0.234  s  =  2.50E-03 

Direction  =  1.000  h  =  0.0150 

FO  •  (-9.993E-11,  9.780E-11)  Number  of  Evoluations:  667 
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Point  =  0.234  s  =  2.50E-03 

Direction  =  1.000  h  =  0.0150 

FO  •  (2.162E-10,  7.314E-10)  Number  of  Evoluotions:  742 


This  plot  corresponding  to  figure 
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Point  =  0.234 
Direction  =  1 .000 
FO  •(1.062,  1.563) 

F2  *(0.1344,  1.145) 
Number  of  Evaluations:  885 


s  =  5.00E-03 
h  =  0.0300 

FI  +(1.266,  1.503) 

F3  0(2.731,  2.731) 


This  plot  corresponding  to  figure 
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Point  =  1.000  s  =  0.0540 

Direction  =  4.00  h  =  0.324 

FO  •  (0.3678,  20.06)  Number  of  Evaluations:  960 
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Point  =  1.000 
Direction  =  4,00 
FO  •  (0.3678,  20.06) 
Number  of  Evaluations:  972 


F6  A  (-60.68,  62.69) 


s  =  0.0540 
h  =  0.324 


This  plot  corresponding  to  figure 
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F6  ^  (-60.68,  62.69) 
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to  figure 
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Point  =  1.000 
Direction  =  4.00 
F2  *(0.1480,8.078) 
Number  of  Evaluations:  989 


F6  A  (-60.68,  62.69) 


s  =  0.0540 
h  =  0.324 
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:-17  s  =  2.11E-04 

lO  h  =  1.26E-03 

1.002)  Number  of  Evaluations;  1353 
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Point  =  0.0 

Direction  =  1 .000 

FO  •  (-1.0000,  9.977E-1 1) 

F2  *(-1.0000,  1.0000) 

F4  X  (-48.00,  48.00) 

F6  ^  (-7.290E  +  03,  7.290E  +  03) 


s  =  0.0333 
h  =  1.000 

FI  +(-0.5000,  7.683E-11) 

F3  O  (-4.000,  4.000) 

F5  □  (-2.430E  +  02,  3.645E  +  02) 
Number  of  Evaluations;  1662 


This  plot  corresponding  to  figure 
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Point  =  0.0 
Direction  =  1 .000 


FI 

+  (-1.0000,  4.284E-1 1) 

F2 

F3 

O  (-2.000,  2.000) 

F4 

F5 

Q  (-81.00,  81.00) 

F6 

Number  of  Evaluations:  1797 

s  =  0.0333 
h  =  1.000 

♦  (-1.0000,  2.851E-10) 
>‘(-8.000,  16.00) 

(-1.458E  +  03,  9.720E+02) 


This  plot  corresponding  to  figure 
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Point  =  0.0 

Direction  =  1 .000 

F2  *(-1.0000,  1.987E-10) 

F4  X  (-2.666,  2.666) 

F6  A  (-1.458E  +  02,  1.458E  +  02) 


s  =  0.0333 
h  =  1.000 

F3  o(- 1.0000,  1.736E-09) 

F5  o  (-8.977.  13.50) 

Number  of  Evaluations:  1932 


This  plot  corresponding  to  figure 
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Point  =  0.0 

Direction  =  1 .000 

F3  O  (-1.0000,  1.600E-09) 

F5  o  (-3.000,  3.000) 

Number  of  Evaluations:  2067 


s  =  0.0333 
h  =  1.000 

F4  X  (-1.333.  1.139E-08) 

F6  A  (-12.84,  27.00) 


This  plot  corresponding  to  figure 
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Point  =  0.0 

Direction  =  1 .000 

F4  X  (-1.0000,  9. 600E-09) 

F6  ^  (-4.128,  4.128) 


s  =  0.0333 
h  =  1.000 

F5  D  (-1.375,  6.609E-08) 

Number  of  Evaluations:  2202 


This  plot  corresponding  to  figure 
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Point  =0.0  s  =  0.0333 

Direction  =  1.000  h  =  1.000 

F5  a  (-1.000,  5.0g3E-08)  F6  a  (- i  .650,  2.666E-06) 

Number  of  Evoluotions:  2337 
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Point  =  0.0 
Direction  =  -1.000 
FO  •  (0.5397,  1.852) 
Number  of  Evaluations:  2424 


FI 
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h 

+  (-1.855,  -0.5406) 


0.0166 

0.1000 


This  plot  corresponding  to  figure 
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CD:  1  Dir  =  (1.000,  -1.000) 
Point  =  (0.0,  1.000) 

Direction  =  (1.000,  0.0) 

FO  •  (-7.491E-04,  7.284E-04) 
F2  ♦  (-4.401,  4.183) 

Number  of  Evaluations:  2598 


ch  =  6.00E-03 
s  =  5.00E-04 
h  =  3.00E-03 
FI  +(0.01050,  0.08002) 

F3  o  (-13.42.  1.002E  +  03) 


This  plot  corresponding  to  figure 
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Point  =  0.0 

Direction  =  1.000 

PO  •  (2. 141 E- 11.  1.800E-03) 

F2  *  (6.166,  1.414E  +  02) 


s  =  4.00E-04 
h  =  2.00E-04 
FI  +  (-0.1824,  0.1824) 

Number  of  Evaluations:  2901 


This  plot  corresponding  to  figure 
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Point  =  0.0 

Direction  =  1 .000 

FO  •  (-1.800E-03,  1.800E-03) 

F2  ♦  (-38.13,  38.13) 


s  =  4.00E-04 
h  =  2.00E-04 
FI  +  (0.01414,  0.1824) 

Number  of  Evaluations:  3052 


This  plot  corresponding  to  figure 

29 


P  L  0  T  7  9 


University  of  Utah 


MICROSCOPE: 
U««r‘a  Manual 


A  Softwor*  Syttom  for  Multivarlata  Analyila  pogo  US 

Vara  Ion  of  April  1984 


Raforaneaa 


1.  Abromowliz.  M. ,  and  Stogun.  1.,  1968,  Handbook  of  Mothanatlcol 

Functions,  Dovor  Publication 


2.  Alfald,  Pator,  A  Trivarlata  Clough-Tochar  Schama  for  Tatrohadrol  Doto 
to  appaar  In  Computar  Aided  Caomotric  Design,  An  Internatlenal 
Journal,  North  Holland 


3.  Besbo,  Nslson  H.F.,  1988,  PQCmCNT.  A  Portable  Text  Formatting 

Pragrom,  Department  of  Physics,  University  of  Utoh,  Solt  Loks  City, 
Utah  84112 


4.  Beebe,  Nelson  H.F.,  1979,  A  User's  Guide  to  <PL0T79>,  Department  of 
Physics,  University  of  Utah,  Salt  Lake  City,  Utoh  84112 

5.  Barnhill,  R.E.,  October  1983,  A  Survey  of  the  Representation  ond 
Design  of  Surfoess,  IEEE  Computer  Graphics  and  Apl 1 1  cot  I  one,  v.  3, 
Number  7,  pp.  9-18 

6.  ds  Boor,  C. ,  1978,  A  Preeticel  Guide  to  Splines,  Springer  Verlog 

7.  Courant,  R. ,  end  Hilbert.  D..  1962,  Methods  of  Mathematical  Physics, 
V.  2,  Wl  ley  8  Sens 

8.  Ryder.  B.G..  The  PFORT  Verifier,  In  Software  Practice  end  Experience 
4.  pp.  389-377,  Wiley  8  Sons 

9.  Schumaker,  L.,  1981,  Spline  Functions:  Basic  Theory,  Wiley  8  Sons 


References 


poge  118 


MICROSCOPE; 
U««r'«  Manual 


A  Softwar*  Syst«Ri  for  Multivariate  Analysis  page  117 

Version  of  April  1984 


Acknowledgements 

MICROSCOPE  was  developed  over  severol  yeors  ond  has  evolved 
through  severol  preliminary  versions.  This  process  was  shaped  by  the  work 
of  the  Computer  Aided  Geometric  Design  Group  at  the  Mathematics  Department 
of  the  University  of  Utah  and  has  benefited  from  the  interoctions  with  Paul 
Arner,  R.C.  8arnhill,  Nelson  Beebe.  Gerold  Eorin,  Tom  Jensen,  Chip 
Petersen,  Bruce  Piper,  Tracy  Whelan.  Andrew  Worsey,  ond  others. 

This  work  wos  supported  portly  by  the  Deportment  of  Energy  under 
Contract  No.  DE-AC02-62ER 12046  ABOB  to  the  University  of  Utah,  by  the 
United  Stotes  Army  under  Contract  No.  DAAC29-8BC-0B41 ,  by  two  grants  from 
the  University  of  Utah  Research  Committee,  and  by  o  sobboticol  leave  of  the 
first  author  from  the  University  of  Utoh.  A  large  port  of  the  program 
development,  and  this  documentation,  were  carried  out  at  the  Mothemotics 
Research  Center  at  the  University  of  Wisconsin.  Fred  Souer  of  MRC  hos  been 
most  helpful  with  all  computer  related  problems. 


I 


Acknowledgements 


page  1 1 7 


MICROSCOPE; 
User's  Menus  I 


A  Software  Systsn  for  Multivoriots  Analysis  pOQO  tfO 

Version  of  April  1984 


Abromowitz,  13,  116 
ACCENT.  51.  54-55. 
168 

Alfeld,  1.  7.  36. 

92.  116 

a  I pho ,  83-85 .  87 , 
89-96 

Bernhi 11,8.  116-17 
Beebs,  6.  8,  69.  78. 
116-17 

bivariate,  9,  16,  36 
BLSCRN.  77 
de  Boor ,  21 ,  116 
C1CROSS.  51.  65-66 
C2CROSS.  51,  65-66 
C3CROSS.  51.  65-66 
CBO.  164 
CBOSW,  164 
CBOU.  164 
CCHANNL.  51,  63. 
67-68 

CDIRCTN,  51.  65-66 
CHVALUE,  51,  65-66 
Clough-Tocher,  36. 
116 

CLSCRN,  77 
coflOTond,  3,  5-6, 

17- 18.  22,  33. 

35-38,  41. 
44-74,  78-79, 
92-94,  161-2. 
164,  167-11 

CORDER,  51.  65-66 
Courant,  45,  116 
critical .  9-16,  27, 
85;  see  olso 
set 

cross,  3.  14-16, 

18- 19.  24,  36. 
44,  51-53,  58. 
63-66,  68, 
72-74,  164, 

166,  169;  see 
also 

derivot ive, 
di rect ion 

CRT,  8.  16,  13.  15 
cusp,  9,  49.  114 


INDEX 

CWINOOW,  51.  60-62 
D,  7.  14.  22-32, 
38-43,  46-49, 
51-53.  56-57, 
59.  64.  67.  83. 
85-86,  88. 
111-14.  116 
d.  10-11,  16.  18. 
51.  54-55. 
59-66.  99-166 
01CH0IR.  51,  59 
02CH01R.  51.  59 
03CH0IR,  51.  59 
OCENTER.  51.  54-55. 
72 

Del  to,  16  .  46  ,  45. 
49-56,  83-85, 
89-91 
DEN,  94 

derivative,  3,  5, 
9-11,  14-16. 

19.  21.  24-27. 

36- 32.  37-66. 
64-66.  76. 
73-74.  82. 
84-86.  88-91, 
98.  165-6,  169, 
113-14 

DF,  165-6 
DGRAPH,  39.  51. 

54- 55.  108 
Dirac,  16,  45,  84 
direction,  3,  16-11, 

14-16.  18. 
22-32.  38-43. 
46-49.  51-53, 

55- 66.  63-68. 
72-74,  96.  99, 
164,  166-8, 
111-14;  see 
also 

i nvest igot ion 
display,  3-6,  8-11. 
13-26,  22-35, 

37- 42.  44,  46. 
49-51,  53-61. 
63,  65-66, 
69-72,  74, 


76-78,  83-84, 
93-98,  165-7. 
169.  111-12, 
115;  see  also 
f unct ion 

DIVIDE,  51.  66-62 
DMNSN,  38.  51.  54. 
58.  65 

DOCUMENT.  6.  116 
DOUBLE.  34,  51, 

66-62.  73.  164, 
111 

DSCALE,  39,  51, 
54-55,  72 
DXAXIS,  51.  54-55. 

72 

e.  8-11,  15-19, 
21-22.  27-28. 
34.  37-39.  52. 
54-55,  58-61. 
63-64,  68-69, 
71,  78.  82-83. 
89,  92,  94,  98, 
164-5,  167, 
116-11.  113 
ECRAPH,  46,  51. 

54-55,  168 
eps,  7,  83 
sto.  36.  42-43, 

82-84.  89-96, 
116-13 

sxaminotion,  3,  16, 
14-16.  19,  23. 
28,  31.  42. 

56-53.  55, 
58-59,  63-64, 
66,  68,  71-72, 
74.  82.  95-96, 
99.  lie,  114 
EXIT,  51.  62-63.  68 
External,  16 
F,  16-13.  34.  73. 
76.  82.  116-11. 
116 

FLIP.  51.  54-56 
Font,  77 

FORCE.  33.  51-52. 

54,  61-64,  73. 


Index 


page  118 


MICROSCOPE:  A  Software  Syaten  for  Multivariate  Analyeie  page  119 

liter's  Monual  Vereion  of  April  1984 


92.  109 

INPEX 

51-53.  55. 

66.  108 

front.  7.  10.  66. 

58-60,  63-64, 

NUM,  94 

112 

66,  68.  72,  74, 

NUMRCL,  34-35,  109 

FS.  96.  100,  105, 

95-98.  107-8, 

nw,  96-97,  99, 

108 

112 

101-4 

function.  2-5,  8-10. 

IPOINT,  51,  58 

OUTPUT.  34-35.  52, 

1>-17.  19-22. 

ITYPE,  92-94 

65.  107 

24.  26-31. 

k.  11.  13.  15-16, 

P.  10-13,  16.  45. 

34-36.  38,  40. 

45.  50-55,  57. 

51-53,  58.  73, 

42-50.  55-56. 

65-66,  70. 

99-100 

61-64.  66.  71. 

82-91 

P1CHPNT.  52.  58 

73-74.  82-84. 

LC.  75-76,  80-81 

P2CHPNT.  52.  58 

89-93.  95-98. 

LCBD.  104 

P3CHPNT,  52,  58 

100.  105, 

LCROSS,  51,  65-66 

package,  3-4,  6-9, 

110-14,  116 

LDEF,  100,  105 

13,  17,  20.  35. 

Gateaux,  18 

LORCTN,  51.  59 

38.  64,  68-69, 

GO.  18,  33,  38, 

LGO,  92-94,  107 

75-79,  92-94, 

51-52,  54,  58. 

LINES.  34-35,  107 

104,  110-14 

61-64,  73.  80. 

LIST.  51.  64-65 

PAUSE,  52,  67-68 

92,  94 

LOG.  51,  65.  105 

PCURSR,  77 

GRAPHC,  34 

LPOINT,  51.  58-59 

PEN07.  70.  75,  78 

h.  5.  11-12.  14. 

LSCRN,  93,  107 

PEN18.  70,  75,  78 

16,  22-32, 

MANLC,  75.  80-81 

PFORT,  8,  116;  see 

38-43.  46-52, 

MANUAL.  75,  79-81 

also  Verifier, 

56-57.  60-61, 

MANUC,  75,  80-81 

Portoble, 

63-64,  66-67, 

MCRSCP,  3,  35-35, 

Portobi 1  i  ty 

83.  85.  87-91, 

58.  62-63,  68. 

PLCHRS,  77 

95-97,  102, 

74,  76,  92-93, 

PLOT.  3.  34-35,  37, 

104,  108, 

101,  104 

52,  69-72,  74, 

113-14 

MICROSCOPE.  1-4, 

77-78,  109-10 

HALVE,  41,  51.  60-62 

6-8,  10-11, 

PLOT79.  3-4,  8.  35, 

Harris.  1.  92 

13-14,  17-23. 

52.  54,  69-72. 

HELP,  34,  51.  58. 

25  .  32-45  ,  49. 

77-78,  80. 

64-65.  79.  109 

52-54.  79. 

115-16 

HELPLC.  75.  79-81 

84.  87.  a2-110. 

PLTDMY.  75.  78, 

HELPUC,  75.  79.  81 

117 

80-81 

Hershey,  77;  see 

Mode,  3.  35-37, 

PLTLC,  75.  77-78, 

also  Font 

52-53.  57.  65, 

80-81 

Hilbert,  45,  116 

68-69,  73,  94, 

PLTUC,  75.  77-78,  81 

HSUMRY.  51.  64-65, 

107,  110-11 

point.  3.  6,  10-11, 

79,  109 

multivariate.  1-3. 

13-16,  18-19, 

ICBD.  104 

8-9,  18-19,  44, 

22-32,  35. 

IDIRCTN,  51.  59 

112 

37-43,  45-49, 

IHVALUE,  51,  68-62 

MULTPLY,  52.  60-62 

51-59,  62-64, 

IINTVL,  38,  51, 

NEWS,  52,  64-65 

66-68,  70-74, 

58-59 

nh,  96,  99,  101-4, 

77-78,  82-85, 

investigation,  3,  7, 

107 

95-100,  105-8, 

10,  14-19,  35, 

NORMAL.  52.  59-60, 

111-14;  see 

Index 


page  119 


MICROSCOPE: 
User's  Manual 


A  Software  System  for  Muffivoriote  Anolyaia  P09*  ^20 

Version  of  April  1984 


a  I  so 

exam  I  not i on 
pole.  9.  113-14 
Portobi I i ty ,  6.  33, 
75 

Portable,  2,  8,  33, 
75-76,  79.  116 
POS,  75-76,  80-81 
precision,  3,  6, 

28-29,  36.  73. 
104 

pr imi live,  8,  10, 

21- 22,  27.  44; 
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